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HOW TO KICK OFF A ROCKET 
in the right direction 


Load! Aim! Launch! And off goes a rocket from shipboard . . . but not 
without precise devices that start the rocket in the right direction. 
Devices ... such as Ford Instrument has manufactured for the 
Armed Forces since 1915. For from the vact engineering and pro- 
duction facilities of the Ford Instrumer Company, come the 
mechanical, hydraulic, electro-mechanica’ » ‘netic and electronic 
instruments that bring us our “tomorrows” tuday. Control problems 
of both Industry and the Military are Ford specialties. 
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FORD INSTRUMENT COMPANY 


DIVISION OF THE SPERRY CORPORATION 
31-10 Thomson Avenue, Long Island City 1, N. Y. 


You can see why a job with Ford Instrument offers young 
engineers a challenge. If you can qualify, there may be 
a spot for you in automatic control development at Ford. 
Write for brochure about products or job opportunities. 
State your preference. 4 
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Summers new PAR Gyro ends guidance headaches. Complicated 
automatic pilots are no longer necessary. Heretofore, position 

and rate signals were obtained by a position gyro and a separate 
rate unit, respectively, then algebraically added by a third unit 
Summers integrating PAR Gyro provides position plus rate on a 
single pickoff. Possessing but one degree of freedom it avoids 
the complexity and limitations of multi-gimbal gyros. Unlimited 
maneuvering may be called for by applying an appropriate com- 
mand voltage directly to the Summers single-gimbal PAR Gyro. 
Drift rates are in the order of 0.1 degree per minute. 


Both weight and cost of the PAR Gyro are about one-fourth 
of the weight and cost of the apparatus it obsoletes. 


Write for details on the PAR Gyro. or for information on Summers’ 
facilities for developing and producing components and systems. 


YROSCOPE HLA. Webb, 34 Mann Steet, Fairborn, Ohio 
Cc. ° M P A N Y 


W. A. Laukaitis, Suite 724, Cafritz Building, 
1625 Eye Street N.W., Washington, D.C 


George E. Harris & Co., Inc., 
2328 BROADWAY * SANTA MONICA, CALIFORNIA 1734 No. Hillside, Wichita, Kansas 


SUMMERS DESIGNS AND PRODUCES ALL WEATHER AUTOPILOTS + RATE GYROS * FREE GYROS * VERTICAL GYROS * TORQUE-TYPE ACTUATORS * POSITIONING-TYPE ACTUATORS * GYRO SERVO 
ACTUATORS * INTEGRATING MOTORS * ALTITUDE CONTROLS * MAGNETIC AMPLIFIERS * FLIGHT TEST TABLES * INVERTERS * PENDULUM POTENTIOMETERS * RATE INTEGRATING DETECTORS * FREQUENCY 
DOUBLERS * MAGNETIC FRICTION CLUTCHES * CONTROL AMPLIFIERS * CONTROL SYSTEM TEST EQUIPMENT * SPECIAL CONTROL SYSTEM COMPONENTS * COMPLETE CONTROL SYSTEMS * RADIO RECEIVERS * 
RADIO TRANSMITTERS * FLIGHT COMPUTERS * CONTROL SYSTEM ANALYSIS SERVICE * ANALOG SIMULATION AND SERVICE * SUMMERS DESIGNS, DEVELOPS AND PRODUCES INSTRUMENTS IN ANY QUANTITY 


71 


POSITION AND 
YRO 
5 
\ 
5 
14 
ly 


GOODYEAR 
AIRCRAFT CORPORATION | 
takes its 
ENVIRONMENTAL TEST PROBLEMS 


Trouble-free performance of electronic circuits, ... smooth operation of 
servo-mechanisms, auto pilots, and complete guidance systems — upon 
these things depend the accurate flight of guided missiles and the accom- 
plishment of their mission. Goodyear is one of the many members of the 
great industrial team behind the U. S. Air Force Guided Missile Program. 
To meet rigid Air Force requirements, many electronic components — as 
well as cross-sections and complete missiles — undergo complete environ- 
mental testing under accurately simulated operating conditions ... in 
this Tenney Altitude Chamber. 


SPECIFICATIONS 
For all types of testing — development, Altitude: 80,000 ft., rated 
research, environment, specification, and Humidity: 20% to 95% 
production—a Tenney-engineeredchamber Temperature: —100°F. to +200°F. 
vill insure dependability and precisely — pyj|.down: to —85°F. with 1000-lb. mass load in 90 minutes 
sontrolled test data for your requirements. MEGS . 
Ser full information write Dissipation: 3 kw at —100°F. and 50,000 ft. altitude 
Tenney Engineering, Inc., Dept. JE, Dimensions: (inside) 7’w x 9’d x 7’h (outside) 10'2’w x 18’d x 10%'h 
1090 Springfield Road, Union, N. J. Temperature and humidity program control instrumentation for this 
chamber produces moisture resistance cycles according to JAN-C-25 or 
Plants: Union, : J and Baltimore, Mil-T-27, and other conditions to meet other government specification 


Ron. cycling tests. Air-tight door slides on tramrail and is) sealed by 

pneumatic closing device. Accessories include frost-free viewing 

te 2 Pree windows, power shafts, terminal pads, provision for inside 

monorail for heavy equipment, and externally operated 
conditioning dampers. 
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MAGINE a motor powerful enough 
to propel a missile at speeds ex- 
ceeding 3000 mph...so pow- 

erful that its developed thrust can 
amount to tens of thousands of 
pounds. One of Ryan’s most chal- 
lenging current assignments is the 

complete production of such a 

motor for an Army Ordnance sur- 

face-to-surface missile. 

There was a many-sided problem 
of fabrication, welding and ma- 
chining that Ryan had to solve in 
connection with rocket motor man- 
ufacture. Because a rocket motor is 
capable of burning as much as a ton 
of fuel a minute at temperatures up 
to 5000 degrees F., terrific internal 
pressures are created that must be 
contained in a very compact pack- 
age of complex construction and 
exacting dimensions. 

The solution was found in per- 
fecting new techniques for forming, 
welding and machining the special 


YAN AERONAUTICAL COMPANY 


ctory and Home Offices: Lindbergh Field, San Diego 12, California 


pment ETHER OFFICES: WASHINGTON, D.C.; DAYTON, OHIO; SEATTLE, WASH.; NEW YORK CITY 


alloy materials. Ryan devised new 


methods of controlling work to very 
close tolerances; ingenious electric 
resistance and arc welding processes 
and a better furnace brazing system 
... plus intricate machine operations 
that had to be jewel-like in precision. 

Ryan’s proved ability in the pro- 
duction of complete rocket motors 
is due in large measure to its long 
experience in building the “hot end” 
of jet and piston engines. Its versa- 
tility in many specialized fields is an 
important advantage in each new 
assignment, for it enables every di- 
vision to draw on 31 years of first 
hand experience in the most ad- 
vanced phases of aviation engineer- 
ing and production. 

Thus, each year more unique 
technical engineering and produc- 
tion projects are awarded Ryan... 
an integrated company with superior 


abilities in meeting the challenges of jf 
today’s high-speed air age. /; ee 
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First carrier-based airplane 


to 7 a world’ s speed record 


Streaking each mile in less than 5 sec- 
onds, during four passes at a 3-kilometer 
course, a Douglas F4D Skyray returned 
he official world speed record to the 
United States. 

Two weeks later, Skyray blazed into a 
100-kilometer course, cracked all records 


Be a Naval Aviator— 
write Nav. Cad. 

Washington 25, D.C. 


by an even wider margin. Now in pro- 
duction for the Navy, this delta wing jet 
interceptor adds terrific climb and fire- 
power to its speed—plus the ease of 
handling needed in carrier landings. The 
Douglas F4D Skyray has now passed its 
initial carrier tests, in service will guard 


our fleets against the fastest of modern 
jet bombers. 

Performance of the Navy’s F4D Sky- 
ray is another example of Douglas leader- 
ship in aviation. Faster and farther with 
a bigger payload is always the basic rule 
of Douglas design. 
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SPINNING, STAMPING, PUNCHING, DEEP 
DRAWING, HYDROFORMING, ANNEALING, 
SPOT WELDING, ASSEMBLING, TOOL MAKING 


Kaupp can supply your instrument cases to exact 
specification quickly and economically. Special 
shapes and odd sizes are a specialty at Kaupp and, 
in most cases, can be turned out on reasonably short 
notice. Kaupp has the experience and the metal 
working facilities for precision forming of intricate 
shapes to close tolerances. Gauges .002 to % stock 
in stainless steel, Inconel, aluminum, cold rolled 
steel, brass and other alloys. Check with Kaupp on 
your metal parts needs, now! 
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CAN HAS BUILT MORE AIRPLANES THAN ANY OTHER COMPANY IN THE WORLD 


an example of North American 
Rocket Research and Development 


Developed by the Propulsion Section, North American Aero- 
physics Department for the U.S. Air Force... this versatile, 
liquid fuel rocket engine produces a thrust of 50,000 pounds. 

Installed on a research sled built by Cook Research Laboratories, 
} the engine accelerates the 5000 pound vehicle to over 
5 1500 MPH in 4.5 seconds to test parachutes and instruments. 
BRE. was =F Many long range programs for developing high thrust rocket 
eae or engines offer you a real career if you are experienced in: 
on : ae Preliminary Analysis & Design Valves — Regulators Turbines — Pumps Controls 
My Stress & Dynamic Analysis Combustion Devices Instrumentation 


You are invited to inquire regarding these professional opportunities. Send a resume fo: 
ENGINEERING PERSONNEL, 12214 Lakewood Bivd., Downey, California 


organization, facilities and experience keep 


North American AVIation Ad. 


in aircraft... atomic energy ... electronics... guided missiles AWA 
research and developement. 
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Eyes in the night . . . certain 
as a pendulum .. . to pierce 
the veil of time and distance. 

Thus we have the moving eye— 
symbol of knowledge — on the — = a 

ever predictable pendulum, : eye sS in the n § ght 


<i Aly + The victory over time and darkness is certain with Kollsman instruments. 
. Certain because of our quarter century dedication to accuracy in controls — 
‘out Tal Today our encompass four fields: 
RADIO COMMUNICATIONS AND NAVIGATION EQUIPMENT 
A 


ko | | S m q n INSTRUMENT CORP. 


ELMHURST, NEW YORK « GLENDALE, CALIFORNIA ¢ SUBSIDIARY OF Standard COIL PRODUCTS CO. INC. 
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Guided missiles launched from submarines promise 
to be major offensive weapons in case of war. A mis- 
sile of this type travels to its distant destination 
under unerring electronic orders. The brain center for 
such missiles will be typical of the electronic systems 
developed and manufactured by Arma Corporation. 

In close collaboration with the Armed Forces 
since 1918, and more recently with the Atomic 
Energy Commission, Arma has contributed much 


in basic research, design, development and manu- 
facture to the advancement of electronic and electro- 
mechanical weapon control, navigation, and other 
precision remote control systems. There is every 
reason to believe that engineering background and 
techniques—first used successfully in these devices 
—will see widespread industrial applications. Arma 
Corporation, Brooklyn, N. Y.; Mineola, N. Y. 
Subsidiary of American Bosch _— 
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By using a few basic principles of compressible flows, it 
is shown that even simple ducts will show interesting prop- 
erlies at supersonic speeds. This knowledge is applied 
to the diffuser problem and a picture is given of the pres- 
sure recoveries and drags to be expected and of the hys- 
that 


given of the behavior of the flow in a ramjet, including 


teresis phenomena will occur. discussion is 
some of the problems of aerodynamic control. It is 
that 
Mach numbers is poor, and ways of using more complex 


A final 


section is devoted to a consideration of some of the flow 


shown the recovery of simple diffusers at high 


diffusers to improve the efficiencies are discussed. 


instabilities that can occur in ramjets. 


Behavior of a Simple Duct at Supersonic 
Speeds 


N DISCUSSING the problem of introducing and expelling 
air from missiles and aircraft at supersonic speeds, I 
believe it is instructive to start with the simplest kind of a 
duct and investigate its characteristics, seeing if it offers flow 
phenomena of interest. We start with a section of cylindrical 
pipe oriented with its axis parallel to a supersonic airstream. 
We imagine that the pipe is so constructed that its internal 
area can be reduced, as shown by the dotted lines in Fig. 1 


FIG. 1 


SCHEMATIC DIAGRAM 


OF A SIMPLE DUCT 


For the moment we keep the picture as simple as possible by 
neglecting skin friction and assuming that the internal flow 
is uniform across any normal cross section. 

For the unrestricted pipe, the flow goes through at the 
speed of the free stream. If the internal area is reduced, the 
velocity of the flow, in contrast to the subsonic case, goes down 
as is known from elementary analysis in compressible fluids. 
As the area reduction is carried farther, the Mach number at 
the throat decreases until the limiting value of Mf = 1 is 


Nore: This paper is based on a lecture given origi- 
nally at the Symposium on Supersonic Aerodynamics held at 
The Johns Hopkins University, December 6, 1945. It. was se- 
lected by the editors for publication at this time in response to 
suggestions that have been received from ARS members to the 
effect that a general introduction to the ramjet diffuser, nonmathe- 
matical yet advanced in viewpoint, would be a welcome addi- 
tion to the literature. 

? Chairman, Department of Aeronautics. 
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reached. Throughout this process, the mass flow has been 
constant and equal to product of the free stream density and 
velocity times the entranee area. However, if the internal 
area is reduced to less than its value when JJ = Lat the throat, 
this condition of constant mass flow is no longer possible. 
The minimum area acts as a “choke” or throttle, reducing 
the flow through the duet. This abhi is communicated 
to the oncoming flow by a compression wave originating at the 
throat and traveling upstream in the duet. It is clear that 
this compression wave cannot remain stationary in the duct. 
It quickly develops into a shock wave and the resulting en- 
tropy increase of the fluid passing through it lowers the pres- 
sure at the minimum section (but does not change the corre- 
sponding velocity or temperature). The resulting decrease 


FIG. 2. 


FLOW PATTERN 


WHEN DUCT IS CHOKED 


Shock 


As a 
consequence, a second compression wave is formed which 
overtakes and coalesces with the first wave, moving upstream 


in density causes a further reduction in mass flow. 


in the duct to a position of higher Mach number. This 
process will be repeated until the shock wave is expelled at 
the entrance and forced far enough upstream so that a portion 
of the oncoming flow can spill around the lips to take care of - 
the reduced flow through the duct (see Fig. 2). 

If the internal area continues to be reduced, the shock wave. 
does not inerease in intensity but moves forward enough to- 
allow a greater amount of spillage around the lips. 

Now let us imagine the procedure to be reversed, i.e., the 
throat area to be increased. The question arises: When | 
the throat area exceeds the value at which the choking first: 
occurs (and the shock appears), is the shock swallowed in the | 
duct from whence it came? It is apparent that it does not 
because the reduced density fluid passing through the throat. 
cannot permit the shock to retreat to the entrance down which 
it must be swallowed. In faet, the internal area must be in= 
ereased by the inverse ratio of the lowered throat density” 
before the shock will be swallowed and the regime of constant — 
mass flow re-established. 

Let us re-examine the above by considering the pressure — 
distribution of the fluid as it passes through the duct. When 
there is no reduction in the internal area, the pressures are 
equal to the free stream value throughout the duct, as shown 
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FOR VARIOUS FLOW 
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FIG. PRESSURE 


by line A in Fig. 3. When the area reduction is moderate, 
the pressure rises to a maximum and falls as shown by B. 
When the area is reduced so that MM = 1 at the throat, the 
pressure is given by C. With further reduction in area, the 
pressures no longer follow this pattern. Instead, a shock 
forms in front and the pressures become nonsymmetrical as 
shown by D. It will be noticed that there is an abrupt rise 
in pressure through the shock, a further subsonic compression 
to the entrance of the duct, then an expansion through the 
velocity of sound at the minimum section to supersonic ve- 
locities in the rear. The pressure at the minimum section, 
corresponding to M = 1, is lower than it was before the shock 
appeared because of the entropy rise through the shock. 
This causes the flow to reach the exit at less than free-stream 


pressure and compression must take place in the exit jet.. 


Continued reduction of the internal area simply exaggerates 
this pressure distribution without changing its pattern as ex- 
emplified by E. 

If, now, the area is increased beyond its value corresponding 
to D, the pressure pattern will not immediately return to the 
symmetrical case typified by C and B but will continue on 
the nonsymmetrical sequence F and G. At G, the area has 
been increased to such an extent that the shock wave has 
moved back to the mouth of the entrance. Any further 
increase in area allows the wave to move back into the duct 
where it becomes unstable. The procedure of originally 
forming the shock will be reversed and the shock swallowed 
by the duct. However, the shock disappearance occurs at a 
different area than that at which it appeared. 

An interesting phenomenon takes place at the exit during 
the sequence DFG. As mentioned previously, when the 
shock first appears D, the flow arrives at the exit with a pres- 
sure lower than free stream and compression must take place 
in the exit jet. However, in the process of increasing the 
throat area, it turns out that a condition F is reached, before 
the shock disappears, for which the expansion from throat to 
exit is so small that the fluid arrives at the exit with free-stream 
pressure. In the limiting case G, the exit pressure is always 
above free stream and expansion will take place in the exit 
jet. Here we have the interesting phenomenon of flow experi- 
encing an entropy rise and yet rejoining the isentropic ex- 
ternal flow, at an increased pressure. 

When we examine the pressure distributions, we see that 
under our assumptions the duct has no drag for the cases 
A, B, and C because the pressures are symmetrical and have 
no resultant rearward component. However, when the shock 
appears in cases D, E, I, and G, higher pressures act on the 
forward facing areas than on the rearward facing areas, re- 
sulting in a large drag. One frequently sees statements in 
connection with tests of ducted models that a shock wave in 
front causes high external drag. I believe that the observed 
drag increase is really atrributable primarily to the internal 
flow and not the effect of the shock on the external flow. . 
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: Hysteresis of Entry Flow 


In our consideration of flow in a simple duct, we have no- 
ticed that there is a hysteresis in the flow pattern as thie 
throat area is changed. This is probably best illustrated by 
the graphs of Fig. 4 showing mass flow and drag as a function 
of the throat area. 

When the throat area is equal to the entrance area, we start 

at the point A, and as the throat is reduced we travel to B. 

During this time the drag has been zero and the mass flow 

constant. At B, we have M = 1 in the throat and the shock 

appears ahead of the duct. The mass flow drops and the 
drag rises to C. A further reduction in area takes us to the 
point D. Upon reversing the procedure and increasing the 
area, the conditions do not follow the line from C to B, but 
the area must be increased to E before the shock retreats to 
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FIG. 4 HYSTERESIS OF DRAG AND MASS FLOW IN A SIMPLE DUCT 


the mouth of the duct. The shock then disappears and we 
go to F, back on the original curve from A to B. It is clear 
that BCIF is hysteresis loop and that we must know the 
previous history as well as the present geometry to ascertain 
the flow conditions existing in the duct. 

Under our previous assumptions of no friction and uniform 

“flow across any duet section, it is a straightforward procedure 
to compute the mass flows and drags for several Mach num- 
bers, as well as the area ratios and pressure ratios as a function 
of Maeh number. This has been done and the results pre- 
sented on the four accompanying graphs of Fig. 5. 

The first graph shows internal drag coefficient versus the 
ratio of throat area to entrance area for the Mach numbers 
1, 2, 3, 4. It will be noticed that the drags are quite high, 
being of the same order as the external drag for projectiles. 

The second graph presents mass flow coefficient versus 
area ratio for the same set of Mach numbers. It shows 
clearly the large reductions in mass flow that occur when the 
shock appears at high Mach numbers. Both the drag curves 
and the mass flow curves illustrate how rapidly the region of 
hysteresis broadens out with increasing Mach number. This 
indieates that at Mach numbers only a little above 1, a small 
amount of unsteadiness in the airstream causes the flow pit- 
tern to shift from one type to the other, in some cases even 
obscuring the existence of a hysteresis phenomenon. How- 
ever, at higher Mach numbers each flow pattern has wide 
ranges of stable existence; so much so that major changes 
in duet geometry are required to go from one pattern to the 
other. 

The third graph presents the ratio of throat area to en- 
trance area for which the shock appears and disappears as 3 
function of Mach number. It shows how great the constric- 
tion must be to produce a shock at high Mach numbers and 
how slight the constriction permissible to insure that no 
shock is present. 

The fourth graph shows how much the pressure is reduced 
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at the throat by the presence of the shock wave. This curve 
‘an also be interpreted as the ratio of the maximum pressure 
recovery possible behind the shock to the recovery if the com- 
pression were reversibly adiabatic. It shows quite clearly 
the enormous losses incurred if the shock is present at high 
Mach numbers. 


Behavior of Ducts at Different Speeds 


In all of our previous discussion we have varied the throat 
area of the duct in a constant Mach number stream. Let 
us now observe the behavior of a fixed duct as the Mach num- 
ber .s changed. This is the case of interest in actual missile 
or aireraft work. We begin our observation at the instant 
when the missile is accelerating through the Mach number 
one. We see from the third graph that if the throat is con- 
stricted at all, the only possible flow pattern is with the shock 
wave in front. This is highly significant, because it means 
that a duet system, with even the slightest constriction, 
enters the supersonic regime with a shock wave in front. 
Consequently, no fixed geometry duct of this type could be 
expected to operate with shock free flow using a throat area 
designed to reduce the Mach number to nearly one at operat- 
ing speed unless the missile is first overspeeded until the shock 
is swallowed and then allowed to return to the operating 
speed (a maneuver not frequently possible in practice). 

It is interesting to compare the results predicted by the sim- 
plified theory given above with some results actually observed 
in practice. The first experiments made on such a simple 
duct were those by Ferri (1). He was investigating a two- 
dimensional biplane proposed by Busemann, who arrived 
at a conclusion similar to ours (by different reasoning) that 
an internal flow system (in his case, two-dimensionally dis- 
posed to form wings) could have very small drags because of 
the nearly symmetrical fore and aft pressure distribution. 
However, when first proposed, it was not realized that the 
nonsymmetrical pressure distribution with shock in front was 
also possible and that this latter case is the one present when 
the supersonic regime is first entered. Ferri’s model, tested 

at a Mach number of approximately 2, had an area ratio of 
0.705, which is in the region of hysteresis. Consequently, 


* Numbers in parentheses refer to the References on page 94. 
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when the tunnel was brought up to speed, the observed flow 
pattern was always the one with the shock present and not 
the one predicted by Busemann, which proved somewhat 
disconcerting. Later it was found that by spreading the two 
wings apart (increasing the ratio A2/A;) the shock could be 
made to disappear suddenly and the Busemann-type flow 
established. 

The drag of the biplane was measured for each of the types 
of flow possible and it was found that Cp = 0.298 when the 
shock was present ahead of the model. Our simple theory, 
neglecting friction, predicts a value of Cp = 0.30. A rea- 
sonable estimate for skin friction (C; = 0.003) increases the 
predicted value of Cp to 0.329 (the ratio of wetted area to 
frontal area is 9.6:1). This agreement is considered encourag- 
ing. When no shock was present, the drag dropped to about 
1/, of the value with shock. Since this low drag ean be largely 
accounted for by skin friction, this result also indicates no 
significant disagreement with theory. 

An interesting discrepancy between our conclusions, based 
on the assumption of uniform flow across any section and 
the actual two-dimensional flow, can be seen in the schlieren 
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FIG. 6 BEHAVIOR OF SHOCK WAVE AT DUCT ENTRANCE 
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photographs of Ferri, drawings of which are reproduced in 
Fig. 6. We predicted that when the shock had retreated to 
the mouth of the entrance it would become unstable and be 
swallowed. Actually, the shock wave retreats a considerable 
distance inside the entrance before it is swallowed. The 
normal shock curves into diagonal shoeks at the ends and 
the entire system is swallowed when very little of the normal 
shock is left. These extremal diagonal shocks correspond 
to the “‘second” solutions, 6 and 6b’ of the shock polar. The 
interior portions of the shock wave correspond to the portion 
of the shock polar between 6 and 6’. After this system is 
swallowed, the completely supersonic flow predicted by Buse- 
mann is established with diagoni il shocks corresponding to 
the “‘first”’ solutions, a and a’ of the shock polar. 


Diffusers for Supersonic Ramjets 


At this point, we attempt to apply some of the results of © 
our previous discussion to the problem of obtaining efficient — 
pressure recovery in an air intake for a supersonic ram jet. 
A glance at Fig. 7 shows that the maximum possible recoveries : 
at high Mach numbers are remarkably large, being of the 
order of 20 or more times atmospheric pressure. This is in 
contrast to the more familiar subsonic pressure rises which are 
only a fraction of an atmosphere. The figure also shows how 
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FIG. 8 DIAGRAM OF RAMJET WITH CONVERGENT-DIVERGENT 
‘DIFFUSER 


much the maximum pressure recoveries are lowered by the 
presence of a normal shock wave ahead of the intake at high 
Mach numbers. Consequently, we are tempted to consider 
the arrangement shown in Fig. 8 

The air is decelerated in a simple converging entrance to a 
Mach number a little above 1 at the throat. The transition 
to subsonie Mach numbers is by means of a normal shock 
maintained in the diverging portion as shown (if the shock is 
in the converging portion, it is unstable as we have seen 
previously). Since the normal shock takes place at lower 
Mach numbers where the losses are small, this arrangement 
‘an provide acceptable recoveries if the subsonic diffusion 
losses are kept low. However, this arrangement is not with- 
out its problems. First, the normal shock must be main- 
tained in its position behind the throat. This can be accom- 
plished either by varying the exit throat area (the final exit 
area cannot be used for this control because in its usable range 
it has no effect on the combustion chamber pressure) or by 
varying the density of the exhaust gases by the heat released 
in the combustion chamber. If the combustion chamber 
pressure is decreased by these controls, the normal shock will 
move back in the diffuser to a point of higher Mach number 
and its Joss will rapidly increase. If the chamber pressure 
is inereased, the shock will move into the throat and, if over- 
controlled, will pass through the throat and out the converg- 
ing section to assume a position ahead of the entrance, nulli- 
fying the advantage of this arrangement. If the control 
is reversed and the chamber pressure lowered, the shock will 
not be swallowed, as we have seen from our previous discussion. 
What oecurs is that the air will accelerate in the converging 
section to M = 1 at the throat and a second shock will form 
in the diffuser. The only way to remove the forward shock 
is either to speed up the vehicle or provide some means of 
varying the ratio of throat area to entrance area. Neither 
expedient can be said to be convenient. 

This forward shock condition is encountered not only 
when the throat shock is inadvertently regurgitated, but is 
sure to be present after the supersonic regime is first entered 
by the vehicle. It is again encountered when the question 
of the stability of the flow is considered. In order to maintain 
low loss, the shock must be positioned near the throat where 
its stability is not great. The inevitable huffing and puffing 
of the combustion process will repeatedly force the shock out 
the intake. After each expulsion the process of swallowing 
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FIG.9 DIVERGENT DIFFUSER FOR LOW SUPERSONIC MACH NUMBER 
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FIG. 10 CONVERGENT-DIVERGENT DIFFUSER FOR INTERMEDIATE 
SUPERSONIC MACH NUMBERS 


the shock must be repeated. All things considered, this 
arrangement does not appear to have too bright a future. 
The foregoing suggests that, in practice, intakes may be 
classed in three Mach number ranges. 
1 In the range of 1<4/<1.5, simple diverging entrances 
such as the one shown in Fig. 9 may be used. In this Mach 
number range, the pressure loss through the shock will be 
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about 10 per cent or less. This simple diverging entrance 
avoids all hysteresis phenomena and eliminates the necessity 
of changing entrance geometry to swallow the shock. The 
greatest recoveries possible with this entrance are given by 
curve 3 of Fig. 7. 

2 In the range 1.5<)/<1.7, entrances may be used in 
which a converging-diverging nozzle is used (Fig. 10), but 
the constriction is so slight that at operating speed, hysteresis 
is always avoided. Reference to the third graph of Fig. 5 
shows that the constriction will be about 10 per cent. This 
modest area reduction proves sufficiently beneficial so that 
the maximum recoveries in this higher Mach number range 
are as great as those of the simple entrance in the lower range. 
When used at speeds lower than the design speed, this en- 
trance will probably have a shock in front most of the time, 
but at these lower speeds, the shock losses are small enough 
to cause no trouble. This type of entrance has been tested 
by Kantrowitz and Donaldson (2) and the preliminary results 
show that it ean be successfully used. The greatest recoveries 
possible with this entrance are given by curve 2 of Fig. 7. 

3 Inthe range M>1.7 it is desirable to use entrances hav- 
ing more complex supersonic flow patterns which are designed 
to avoid the troubles we encountered in our high efficiency 
converging~liverging entrance. We shall discuss such en- 
trances at greater length in the following paragraphs. 


Diagonal Shock Diffusers 


In our original discussion of flow in a simple duct, we postu- 
lated a simple one-dimensional flow. This postulate has 
been retained throughout the above discussion and this reten- 
tion has foreed us to overlook some very useful means of 
decelerating a supersonic stream efficiently. The unidimen- 
sional postulate permitted only normal shock waves which are 
unaeceptably inefficient at high Mach numbers. However, 
if we consider the possibility of using diagonal shock waves, 
this proves to be a fruitful field of investigation. We find 
that a single diagonal shock wave can reduce the Mach num- 
ber from 2.5 to 2.0 with only a 5 per cent loss in recoverable 
pressure. If we place a slightly constricted entrance in the 
flow behind a diagonal shock, as shown in Fig. 11, the limiting 
recoveries are about 81 per cent. If a diverging entrance 
alone is used the limiting recoveries are about 50 per cent. 
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FIG. 12. DIFFUSER USING MULTIPLE DIAGONAL SHOCKS 


A slightly constricted entrance by itself gives limiting recov- 
eries of about 62%. The advantage of using the diagonal 
shock is quite apparent. 

With such gains possible, it appears worth while to consider 
using multiple diagonal shocks, possibly followed by a simple 
diverging entrance as shown in Fig. 12. At M = 2.5, this 
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MACH NUMBER. 
FIG. 13° PRESSURE RECOVERY OF DIAGONAL SHOCK DIFFUSERS 
will not improve the recovery greatly; however, at higher 
Maeh numbers this procedure can be used quite effectively 
to maintain high levels of efficieney, as shown by Fig. 13. 
This was taken from a paper by Oswatitsch (3). It shows 
the limiting recoveries for multiple diagonal shocks and illus- 
trates clearly the advantage of their use at high Mach num- 
bers. 

This type of entrance was first discussed by Ferrari, Oswat- 
itsch, and Busemann. During World War IT, tests were made 
in Germany at the Kaiser Wilhelm Institut, G6éttingen, The 
Hermann Goring Institut, Braunschweig, and the Wasser- 
bau Versuchsanstalt, Kochelsee, and the results indicated 
that superior efficiencies could be attained in practice. A 
few of Oswatitseh’s results are indicated in Fig. 13. In all of 
the afore-mentioned tests, the entrance was constructed with 
circular symmetry by revolving the entrances shown above 
about the straight side parallel to the airstream. 

The great disadvantage that this type of entrance appears 
to possess is the difficulty of keeping the drag within reasona- 
ble limits. High recoveries and low drags are somewhat in- 
compatible. It will be remembered that the drag of a given 
surface element increases approximately as the square of its 
angle of inclination. On the other hand, high recoveries 
dictate large angles of inclination of the entrance. Some im- 
provement ean be had by reducing the projected area of the 
lower lip with the shocks focusing at the tip of the upper lip, 
as shown in Fig. 14. It is also helpful to curve the upper lip 
as shown. However, this must be done with care to prevent 
difficulties in the subsonic diffuser. In spite of these refine- 
ments, the drags are still large. 

In searching for a means of reducing the drag, the question 
If multiple diagonal shocks are used, do they all 
need to come from the same surface, or can they come from 
opposite surfaces, as shown in Fig. 15? The reduction in 
frontal area and drag is immediately apparent. it will be 


FIG. 14 DIAGONAL SHOCK ENTRY WITH CURVED DIFFUSER 
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offered little gain over multiple diagonal shocks, and the for- 
mer has the construction difficulty of a very fine point. 


Unstable Flow in Ramjets 


_¢ At this point, we turn our attention to the effects of the 

7: 15 entrance on the internal characteristics of a ramjet, and in 
1G. 15 DIFFUSER USING REFLECTED DIAGONAL SHOCK 
olen ‘ particular the role the entrance can play in promoting oscil- 


~ yemembered that the important parameter in determining lations in such a system. There are many ways in which a 
hysteresis is the ratio of entrance area to throat area (when ramjet ean have dynamic coupling that will produce. oscilla- 
this ratio is one, as in the simple diverging entrance, no hys- tions, particularly as its control and fuel metering systems 
teresis is experienced). This type of entrance is evidently become more complex. Only a few simple examples will 
somewhat inferior in this respect to the higher drag entrances be treated here to illustrate the possible difficulties that ean 
discussed above. be encountered, 

In examining the above entrance, it might be suspected Our first example considers a ramjet operating with a : 
that the interaction of the second diagonal shock with the simple diverging entrance. The ramjet is to be ignited in js 
sizable boundary Jayer of the lower lip might cause trouble. flight, but when tried this cannot be successfully accomplished. ti 
However, this is not the case because the tendency is to di- The test data have shown that the mixture strength is too |, 
minish the thickness of this boundary layer rather than cause lean at the instant of ignition. To remedy this, the fuel s 
separation. Such cases were encountered by Ferri (1). flow has been considerably increased.- The burner now ig- " 
He found that when the shock does not intersect the bound- nites properly but combustion temperatures increase the | (, 
ary layer at the kink in the surface, the boundary layer (lis- back pressure of the exit nozzle so much that air flow is re- | y 
torts so that at its outer surface the kink is obligingly moved duced to about one third of its value before ignition. The ” 
to the shock intersection point, as shown in Fig. 16. This mixture strength now becomes so very rich the fire is extin- | 
helpful role of the boundary layer prevents the shock wave guished and the cycle begins over again for as long as the ig- | |, 
from being reflected by the solid surface and followed or nition system is left on. eo 
preceded by an expansion fan when the shock does not exactly In our second example, we suppose that the ramjet is } ,, 
intersect the kink as designed. equipped with a slightly constricted entrance in which the | o¢ 

throat is large enough to avoid hysteresis. The fuel for the | gj, 


Shock Exponsion ramjet is supplied by constant flow nozzles which, by design or | gi, 


a Waves Woves accident, are delivering an amount of fuel in excess of stoichio- | a, 
_ Shock behind : metric proportions. The exit throat, normally adjusted so | jg, 
that a normal shoek would be just behind the entrance throat 
. (as shown in Fig. 18), is actually constructed slightly smaller. | gy, 
4 ¢ ‘ . This forces the shock wave through the throat and out in | g¢y 
layer front. As soon as this Scouts, the air flow diminishes ei 
Shock cheod of the greater-losses through the shock in front, causing a] my, 
ot Kink decrease in the air to fuel ratio and the temperature of com: | tio 
ae bustion. This in turn decreases the back pressure of the exit | pe. 

; nozzle, allowing the shock to be swallowed. However, in | js 
; the beginning we saw that this condition was unstable and | te, 
FIG. 16 SCHEMATIC DIAGRAM OF SHOCK-WAVE BOUNDARY-LAYER started the eycle of oscillation. Fo 

INTERACTION 

dis 
Having gone this far in the use of diagonal shock waves to the 
improve the pressure recovery characteristics of supersonic stre 
intakes, it is only natural to ask about the possibility of using anc 
diagonal compression waves, as shown in Fig. 17. In this the 
case the extended lip is smoothly curved and isentropic com- alw 
pression takes place without shock losses. A glance at Fig. whe 
13 shows that such an entrance would have a useful applica- 7 cer! 


phy 
FIG. 18 TYPICAL RAMJET CONFIGURATION dev 
an 
For our last example we consider a ramjet missile having | tro 
7 an entrance with a large constriction. It is operated in the 
range of hysteresis. If we suppose that the burner used 
operates poorly on lean mixtures and that the drag is less 
than anticipated in the design calculations, then we have the] 4 . 


following possible sequence of events. After acceleration to] ¢, - 


FIG. 17 OSWATITSCH DIFFUSER USING COMPRESSION WAVES 
. design speed, a shock wave is present in front of the en-]| D - 


tion at very high Mach numbers where the number of diagonal trance. However, the thrust still exceeds the drag and the] f= 
shocks required to obtain good recoveries would be so great missile continues to accelerate with the drag, increasing only | ?> 
that they would be almost indistinguishable from a uniform ~ slightly more rapidly than the thrust. The missile reaches] % = 
fan of compression waves. The low losses possible with this == =a Mach number at which the shock is swallowed. This re 
configuration makes its use attractive at lower Mach numbers sults in a large air flow because of both the increased re} —— 
as well. covery and the high Mach number. The mixture ratio 
Oswatitsch made early tests on the entrance (see Fig. 17) becomes quite lean and the burner operates poorly or not at], 7 
which makes use of a fan of compression waves rather than all. The thrust decreases abruptly and the missile deceler-} 14 
discrete diagonal shocks. No unexpected difficulties were ates toa Mach number where the forward shock reappears, 
encountered. However, at the Mach numbers of his tests starting the cycle over again. 
(about 2.9), Oswatitsch felt that the compression waves (Continued on page 94) 
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The analytic investigation of laminar combustion 
processes which are essentially two- or three-dimensional 
present some mathematical difficulties. There are, how- 
ever, several examples of two-dimensional flame propaga- 
tion: which involve transverse velocities that are small in 
con: parison with that in the principal direction of flow. 
Such examples occur in the problem of flame quenching 
by « cool surface, flame stabilization on a heated flat plate, 
combustion in laminar mixing zones, etc. In these cases 
the problem may be simplified by employing what is 
known in fluid mechanics as the boundary-layer approxi- 
ma'ion, since it was applied first by Prandtl in his treat- 
ment of the viscous flow over a flat plate. Physically it 
consists in recognizing that if the transverse velocity is 
small, the variations of flow properties along the direction 
of main flow are small in comparison with those in a 
direction normal to the main flow. The analytic descrip- 
tion of the problem simplifies accordingly. The present 
analysis considers the ignition and combustion in the 
laminar mixing zone between two parallel moving gas 
streams. One stream consists of a cool combustible mix- 
ture, the second is hot combustion products. The two 
streams come into contact at a given point and a laminar 
mixing process follows in which the velocity distribution is 
modified by viscosity, and the temperature and composi- 
tion distributions by conduction, diffusion, and chemical 
reaction. The decomposition of the combustible stream 
is assumed to follow first-order reaction kinetics with 
temperature dependence according to the Arrhenius law. 
For a given initial velocity, composition, and temperature 
distribution, the questions to be answered are: (1) Does 
the combustible material ignite; and (2) how far down- 
stream of the initial contact point does the flame appear 
and what is the detailed process of development. Since 
the hot stream is of infinite extent, it is found that ignition 
always takes place at some point of the stream. However, 
when the temperature of the hot stream drops below a 
certain value, the distance required for ignition increases 
so enormously that it essentially does not occur in a 
physical apparatus of finite dimension. The complete 
development of the laminar flame front is computed using 
an approximation similar to the integral technique in- 
troduced by von Karman into boundary layer theory. 


Nomenclature mas: a 


A = activation energy, cal/mole/°K lente 
= specific heat at constant pressure, cal/gm/°K 
D = binary diffusion coefficient, em?/sec ; : 
f = dimensionless stream function 
P, = Prandtl number = cpu/d 
@ = rate at which heat is added per unit volume, cal/sec/unit 


Presented at the 8th ARS National Convention, New York, 
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cial a of the Ordnance Corps, U. 8. Army. 
2 Associate Professor, the Daniel and Florence Guggenheim 
Jet Propulsion Center. 
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gq = average heat release per unit mass of fuel gas, cal/gm , 


R = universal gas constant, cal/mole/°K 7 

S. = Schmidt number = y/pD 
T = local stream temperature, °K | 
u = local axial stream velocity, cm/sec w 
v = local transverse stream velocity, cm/sec ; 
8 = boundary of thermal zones 

A = boundary of diffusion zones 

6 = boundary of velocity boundary layer 
= dimensionless temperature = 7'/T» 
= dimensionless activation energy = A/RT> 


kK local relative mass concentration of combustible gas a 
A = velocity ratio = w2/u at 
\ = coefficient of thermal conductivity, cal/em/sec/°K > 
u = dynamic viscosity, gm/cm/sec 
vy = kinematic viscosity, cm/sec = u/p 
= density, gm/cm? 
7 = characteristic time of the chemical reaction 

Introduction 


NALYTIC studies of combustion in flowing gaseous 

mixtures (1, 2, 3)* are confined, for the most part, to the 
plane one-dimensional flame front. The reason for this is not 
entirely because of the relative simplicity of this problem but 
also because, due to the very small thickness of the laminar 
combustion zone at ordinary pressures, many actual situations 
are well approximated by the plane flame front in spite of the 
fact that the actual flame front is curved. 

There exist some problems, however, in which the physical 
situation may in no way be approximated by a plane flame 
front. For example, the thermal quenching of a gaseous 
mixture near a cool wall (4) or the thermal ignition near a 
heated wall essentially involve two-dimensional fields. Like- 
wise, combustion processes in free jets and combustion under 
conditions of mixing between two gaseous streams require 
considerations beyond the one-dimensional flame theory. In 
many problems associated with thermal jet-propulsion sys- 
tems and rocket motors, the circumstance arises where a 
flowing stream of combustible mixture is ignited by contact 
and mixing with another stream of hot gas or products of 
combustion. Such a process is certainly of importance near 
the injector of a rocket motor, in the flow and combustion of 
gas through the small ports of turbojet cans, and plays a 
definite and vital role in the flame stabilization on bluff bod- 
ies. This, in common with the other two-dimensional prob- 
lems which have been mentioned, has the property that, for 
large stream velocities the variations of temperature, com- 
position, and velocity are much larger in the direction normal 
to the stream than they are parallel to the direction of flow. 
In fluid mechanics, problems exhibiting this characteristic are 
treated by the so-called boundary layer approximation which 
simplifies the description of the problem by deleting certain 
variations along the direction of flow in comparison with those 
across the flow. An extension of this idea to flow with com- 
bustion allows an analytic treatment of this entire class of 
problems. The present treatment is restricted to laminar 
flow processes, in particular that of the ignition and develop- 
ment of a flame front in the laminar mixing zone between 
parallel streams of combustible gas and hot products of com- 
bustion. 


_ 4 Numbers in parentheses refer to the References on page 94. 
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After discussing the physical and chemical relations which 
apply to the problem and simplifying the equations through 
extension of the boundary layer concept, the initial portion of 
the mixing zone is investigated where the heat evolved through 
’ chemical reaction is yet small and the problem may be solved 
through a perturbation to the solution for mixing without 
combustion. ‘Then using the integral technique introduced 
by von K4érman into the study of the boundary layer, the 
Af development of combustion is traced approximately from 
J ignition through development of the plane laminar flame 
front. 
‘The authors have benefited from their many opportunities 
of discussing this work with Professor H. S. Tsien, through 
whose advice the problem has become much clearer and 
greatly simplified. They are most grateful for his patient 
and continued interest in the analysis. 


ah. 
be Formulation of the Combustion Problem _ 


The problem of combustion in a laminar mixing zone may 


be formulated by considering two semi-infinite streams of gas 
flowing steadily parallel to the z-axis, one of a combustible 
mixture at temperature 7; flowing with a velocity ™, the 
other consisting of combustion products at a temperature 7», 
less than the adiabatic flame temperature for the combustible 
mixture, moving with a velocity ws. Suppose the two 
streams begin to mix at « = 0. Then if no reaction takes 
place the process is simply one where the velocity profile 
develops in a well-known manner (Fig. 1) and the temperature 
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VELOCITY AND TEMPERATURE PROFILES FOR A MIXING 
REGION WITHOUT CHEMICAL REACTION 


profile grows in a somewhat similar fashion. In addition, 
diffusion takes place between the two streams of different 
gaseous composition, the combustible gas from the upper 
stream diffuses into the hot combustion products of the lower 
stream. Likewise products of combustion diffuse from the 
lower stream to dilute the combustible upper stream. 

Now when the combustible stream reacts chemically, 
producing heat in the process, the change in the resulting flow 
is most easily observed from the development of the tem- 
perature profiles. As indicated in Fig. 2, the first noticeable 
modification will be the small increase in temperature caused 
by reaction of the gas which has diffused into the hot stream. 
The temperature increase takes place predominantly in the 
hotter part of the stream because of the extreme importance 
of the ambient temperature in promoting reaction. At a 
later stage the temperature profile is sufficiently distorted to 
produce a ‘‘bulge’” where the local temperature exceeds that 
of the hot gas stream. From this point forward, the strongest 
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- reaction takes place above the maximum temperature p»int, 


_ development of the temperature field through conduction be- 
- tween the two streams and heat generated by chemical action; 


| 


PURE MIXING 
COMBUSTION 


FIG. 2 TEMPERATURE PROFILES AND DEVELOPMENT OF LAM!NAR 
FLAME IN MIXING REGION WITH CHEMICAL REACTION 


due to the fact that the greater concentration of combus: ible 
gas is exposed to the higher temperatures, and consequently 
the point of maximum temperature moves toward the upper 
stream. Eventually this zone of reaction develops inio a 
thin region resembling the laminar flame front and propagates 
at the appropriate speed. Since the mixing zone grows 
roughly as VJ. 2, it is clear that the reaction zone will move out 
of the laminar mixing region into the undisturbed combustible 
mixture and hence develops into a true plane laminar flame 
front. If the velocities of the two gas streams are such that 
laminar flame makes a very small angle with the z-axis, that 
is, if wu; >> laminar flame speed, then the boundary layer 
approximation is applicable throughout the region of flame 
development. 

Three physical processes are involved in this problem: (A) 
The fluid mechanical development of the laminar mixing zone 
between the two streams through viscous action; (B) the 


(C) molecular diffusion of each gaseous component into the 
other stream together with annihilation of combustible gas 
and generation of combustion products through chemical 
action. The first of these processes is described by the equa- 
tions of motion for a viscous compressible fluid. Under the 
boundary layer assumptions only the dynamic equilibrium in 
the direction of flow need be considered. Furthermore, it is 
well known that due to the low value of the laminar burning 
speed, the pressure variations are negligible. Therefore the 
pressure field is uniform and, denoting by a, 6 the velocity 
components in the z, 7 directions, the equation of motion is 


ay \Y ag 


Here yu is the viscosity of the gas so that if it is assumed that 
the molecular weights of combustible and combustion 
products are the same’, the viscosity is simply proportional to 
T, the local gas temperature. Now since the gas velocitie 
are very small in comparison with the velocity of sound, the 
density p depends only upon the local gas temperature, 
in fact p ~ 1/T. The two velocity components and the 
density are related further through the equation of continuity, 
expressing the conservation of mass as 


Now since the density is a function only of the temperature, 


5 The three simplifying assumptions underlying this choice of 
gas properties are: (1) The molecular weights of both gases are 
equal so that the molecular weight of the mixture is independent 
of concentrations; (2) the gases are composed of Maxwellian 
molecules so that », for example, behaves as 7; (3) the Eucken 
correction for thermal conduction holds for the mixture. 
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it is particularly convenient to employ the Howarth trans- 
formation which introduces a new vertical scale (5, 6) 


where p; is the density of the combustible free stream. Call- 

ing the new coordinate system z, y where z = z, and new 

velocity components u = a, v = + af, = (2 )as, it is easy 
pi 0 0% \p 

to show that the modified flow satisfies the boundary layer 

relations for a corresponding incompressible flow problem 


or oy 
ou 4 Ou (5) 
Or dy Oy? 


where v = u/f is the kinematic viscosity of the combustible 
stre:m. Formally then the velocity field u, v is independent 
of the temperature and composition fields and may be deter- 
mined separately. The calculation of the real field 7, +, how- 
ever, requires knowledge of the temperature field to invert 
the transformation given by Equation [3]. 

The temperature field is described by the energy equation 
(7,8) which equates the convection of thermodynamic enthalpy 
with the transmission by conduction and the heat generation 
by chemical reaction. Then assuming the value of ¢,, 
specific heat at constant pressure, to be the same for all gas 
compositions and independent of temperature, the differential 
relation is, neglecting the pressure variation 


de,pT 9d ( 
5 5 — 6 
+ pi ag +Q [6] 


where Q is the rate at which heat is being added per unit 
volume at a point z, 9, and \ is the thermal conductivity of 
the gases which will be assumed proportional to the local 
temperature. The heat addition rate Q is determined by the 
heat of reaction q,* the local mass concentration of the reacting 
molecules x, and the rate at which molecules of fuel are trans- 
formed into molecules of combustion products. This rate is 
usually accepted by chemists to be of the form (1/r)e~ 4/7, 
the Arrhenius factor, where 7 is the characteristic time asso- 
ciated with a transformation, A the activation energy which 
measures the kinetic energy a molecule must possess in order 
to react successfully, and # the universal gas constant. Fora 
unimolecular reaction, the rate of chemical transformation is 
simply the product of the local density of combustible gas 


with the Arrhenius factor, that is, xo(1/r)e~4/"7, Thus the 
rate of heat addition by chemical reaction is simply 
Q = qpx - e~A/RT (7] 


It is again convenient to apply the Howarth transformation 
to the energy equation so that it will be written in terms of the 
same quantities as the momentum equation. 
oT P 


oT + oT 
v— = — 
Ox Oy P, oy? cpr 


where P, = c,ui/\; is the Prandtl number of the gases, which 
isconstant. For nonvanishing temperature the heat evolved 
through chemical reaction vanishes only if the local concen- 
tration of combustible matter vanishes. This fact leads to 
the formal difficulty that a gas which flows from z = — @ is 
completely reacted unless a condition is added to disconnect 
the flow from x = — ©. The identical problem arises in the 


‘In view of the assumptions concerning the constancy of 
molecular weight and cp, g is to be interpreted as the average 
heat release per unit mass of the fuel gas. 
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study of the one-dimensional laminar flame and has recently 
been discussed by von Karman (3). In the present work it 
will be treated in a manner similar to that which von Karman 
has employed for one-dimensional flame propagation. 

Finally the distributions of combustible gas and products of 
combustion must be determined through the conservation 
relations for the chemical species. There are two species 
and hence two conservation laws; however, it is clear that the 
two must sum to the gross continuity relation, Equation [2], 
and hence only one is independent. It is convenient to con- 
sider the mass concentration «x of the combustible gas. The 
movement of this gas differs from that of the complete gas 
only by the diffusion velocity and the combustible mixture is 
destroyed locally according to the reaction rate. Then with 
good approximation 


OK Ox OK K 
pu + pi ¢ ) p- € [9] 


If the Howarth transformation is again carried out, and the 
Schmidt number S, = »/pD assumed to be constant, then the 
conservation of combustible material may be written 


[10] 
oy? T 
Equations [4], [5], [8], and [10] now describe the simplified 

combustion process which takes place in an equivalent in- 

compressible flow. By virtue of the Howarth transformation 
the momentum and continuity equations may be treated 
independently of the temperature and composition fields and 
hence constitute a conventional mixing problem. The energy 
relation and equation for conservation of combustible ma- 
terial must be solved simultaneously, employing knowledge 
of the velocity field. Since the equations are nonlinear it is 
clear that solutions may be expected in only certain cases or 
certain regions and only approximate solutions may be 
obtained with any generality. Because of the importance of 
the ignition region, a solution will be obtained for the early 
portion of the mixing zone. Then, by means of an approxi- 

mate treatment, the complete course of development of a 

laminar flame front will be examined. 


Initial Development of the Combustion Zone 


During the first portion of the combustion process, where 
the heat added by the chemical reaction is yet small, the 
velocity, temperature, and composition fields differ only 
slightly from those which would occur in the absence of chemi- 
cal action. The fields without combustion may be computed 
easily, and consequently the variation due to combustion may 
be treated as a perturbation. 

If the reaction were taking place at a temperature 7 and 
concentration x, the combustible material would be consumed 
at a rate x(1/r)e~4/"?. Due to the strong influence of the 
Arrhenius factor, the initial reaction takes place where T' is 
nearly 7’, and the concentration x is very small. The actual 
value of x where the reaction takes place most rapidly may be 
estimated by computing the maximum value of the local 
reaction rate. Thus | 


Ald? garig 


RT T dy 


dy r T 


But the temperature distribution (T, — T)/(T2 T;) is 


identical with the concentration distribution «(y) when the 
Then 


Prandtl and Schmidt numbers are equal. 


de /aT 
dy/ dy T, — 


and taking account of the fact that the reaction takes place 
where T ~ T>, the local concentration is approximately 
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Consequently the local rate at which ciateadindn matter is 
consumed is 


RT 


1 A/RT: “Ab. 


and the total heat release which takes place over a length x 
is proportional to 


e—A/RT: 


A 
T,/T>) 


RT? 
since the reaction takes place predominantly where the gas 
velocity is nearly u:. Then so long as z is such that the total 
heat addition is small, a perturbation procedure is in order. 
Therefore it is appropriate to define a characteristic length 
TU 
(1 — RT», 
which may be employed to make the formulation dimension- 
less; denote ¢ = 2/l. 

Because of its important role in the early stages of combus- 
tion, the temperature 7. of the hot stream is conveniently 
chosen as the characteristic temperature; denote T7/T. = ¥, 
T;/T: = and A/RT, = The energy and diffusion 
equations then may be written simply as 


associated with combustion, / = 


and 
2 
Ox oy Se Oy? T 


Now the momentum equation has a simple similarity solution 
regardless of the structure : the rest of the field. For, 
choosing a stream function (x, y) such that u = Oy/dy, 


it is easily shown that = where 
y V/u;/va and f(n) satisfies the differential equation. 


v= 
7 = 


with boundary conditions f’(~) = 1, f’(— ©) = — = A,and 


= ©) = 9. 


In the absence of chemical reaction 


Equations [12] and [13] also have simple solutions. The 

differential equations reduce directly to 

2 


e 


with boundary conditions ©) = = 3, ©) = 
1, = = Q; and similarly = 1, 
KO(— co) = 0, = = 0. 

The presence of chemical action destroys the elementary — 
similarity of this simple mixing problem so that the solutions - 
3 and « cannot be functions of 7 alone, but the distance z 
must also enter. Then, for example, if # depends upon both 7 
and x, Equation [12] becomes 


which reduces heel to Equation [15] when g = 0, for then. 
the dependence upon z is unnecessary. By writing Equatio 1 + P, erf vP 
[17] in the form & 
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do 


1 
the appropriate expression for 8 is 
= + + + +... [19] 


The situation is quite similar for the equations describing the 
conservation of combustible matter, for Equation [13] may 
be written 


1 
Se On? 


fo K 


so that likewise it is appropriate to express 


It is obvious that 8(n) and x(n) satisfy Equations [15] «nd 
[16], respectively, for the pure mixing problem, while the 
relations to be satisfied by 3(n) and «™(y) follow upon 
substitution into Equations [18] and [20], respectively, 


according to the classical perturbation technique. The first- 
order perturbation gives simply 
P, dn? 2 dn 2 
1 f dx -1) 
—— — = — d1) a0 —") (23 
S. dn? 2 dn fe 


where it is necessary that these perturbations and their de- 
rivatives vanish at 7 = +o. 

The perturbation solution may be worked out for various 
special cases, the simplest of which is that where the velocities 
of the two streams are equal, that is, where A = 1. Then the 
stream function is just f = », and it is easily demonstrated 
that the distributions of temperature and composition in the 
absence of combustion are just 


1— a 2 


KO(n)—1 = ; — erf (; V 8. [25] 


where erf(z) = (2/ dt. These are well-known 
solutions of the simple heat conduction and diffusion prob- 
lems. The first-order perturbation equation for temperature 
distribution, Equation [22], may now be written as 


P, dyn? 2 dn 
q TU 

KOI (1 — 300) . [26] 


where #(n) is given by Equation [24]. The homogeneous 
equation has two linearly independent solutions which may be 
written in the form 
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P, 


Then constructing the appropriate Green’s function as 


4 xP, 
Gn, 7; Pr) = [29] 
| =, 


the solution of Equation [26] may be written explicitly as 


C,T: 

f G(n, (Ae 


It is clear from Equation [23] that the first-order perturbation 
for the concentration may be obtained in a similar manner as 


1 


-1) dij. . . [30] 


v,) 
1 
f G(n, dy. . . [31] 


where, as indicated, the Schmidt number now replaces the 
Prandtl number in Green’s function. 

This form is very suitable for numerical computation since 
Green’s function and the rest of the integrand are easily 
calculable. The integration has been carried out for the 
particular case where P, = 0.91, S. 1.00, 3; = 0.286, 


J, = 23.96, va = 1.29, and 7, 300°K, where physical 

parameters are approximately those of the fuel azomethane. 
The resulting temperature profiles are given in Fig. 3 where 
itis clearly shown that the maximum perturbation occurs well 
low the center line of the mixing zone. This is due to the 
relatively preponderant effect of the stream temperature in 
increasing the reaction rate, so that although the concentra- 
tion is low, the temperature is sufficiently high to react almost 
immediately that combustible which has diffused into this 
region. The temperature profiles show that very early a 
bulge is formed where the temperature exceeds that of the 
hot stream. This region continues to grow, working its way 
toward the combustible stream. The tendency is therefore 
for the initial reaction to develop into a flame since the tem- 
peratures continually increase with z. 

The principal question is how long a time or, more realis- 
tically, how great a distance is required before the combustion 
takes place. The simplest significant length to measure is 
the value of z at which the temperature profile first exhibits a 
vertical tangent; that is, where the local temperature has 
increased sufficiently to form a ‘‘bulge” in the profile. Fortu- 
uately this point may be determined analytically, for con- 


dy 
dn 


dn 


Furthermore the smallest value of ¢ at which this condition 


Sc 
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may be satisfied is the dimension of interest and will be de- 
0 at &;, so that 


noted ¢,. It is necessary also that dt/dy 
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DURING THE EARLY STAGES OF REACTION 


calculating from Equation [32] 


i - = 0 b 
dn? dn? 
But now applying the differential relations for 3(n) and 


3(n), given by Equation [15] and Equation [22], respec- 
tively, it follows that the value y; for which a vertical slope in 
the temperature is first observed is given by the simple alge- 
braic relation 


1 
= —— — de? [33] 
Then, according to Equation [32], the value of &; is just 
dyn dn /ni 


Now the expressions involved in Equations [33] and [34] are 
not simple to handle, but fortunately the values of n; which are 
involved may be seen, from Fig. 3, to be < — 0, so that good 
approximation may be made using asymptotic expansions for 
3 (ny) and for the right-hand member of Equation [33]. For 
example, when the Prandtl number is less than unity, the 
asymptotic expressions for «(n) and 3(y) at large negative 
values of » are 
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Then from Equation [33] it follows that the value of y; is determined by the relation 
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where the asymptotic expansions have been used for dd“ /dy. 
Now the principal variation of z; is due to the exponential 
term so that it is clear that as the temperature of the hot 
stream decreases, the detachment distance of the combustion 
zone increases exponentially. Using the same values of the 
parameters as those employed in constructing the temperature 
profiles, the detachment distance x; has been computed for a 
range of hot stream temperatures. Calculations incidentally 
showed that 7; was about —6, thereby justifying the use of 
asymptotic expansions in its evaluation. It appears from 
Fig. 4 that the value of x; increases enormously as the hot 
stream becomes cooler. Therefore it is seen that although 
this process of combustion in a laminar mixing zone shows no 
distinct blowoff, the flame detachment becomes so great for 
low stream temperature that it exceeds the physical dimen- 
sions of any apparatus. 

Because the initial reaction occurs at such large negative 
values of 7, it is a simple matter to extend the analysis to the 
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After numerical solution for n;, the value of &; is obtained through direct substitution into hii [34], that is 


(8. — -1 
12] 1 1 


case where A ¥ 1 by employing the asymptotic expansions of 

the solution for the velocity field. However, this additional 

consideration modifies Equation [35] only to the order 

(e 4)?, and therefore Equation [37] is valid, to the present 

approximation, for any value of A that is not too small. 


Development of the Laminar Flame 


The perturbation scheme employed in discussing the initial 
combustion zone is, of course, inappropriate for following the 
transition to a laminar flame front. However, since the 
initial region is known from the foregoing analysis, an approxi- 
mate procedure is employed in extending the solution 
throughout the remainder of the flame zone. The technique 
that allows best use of the information which has so far been 
accumulated is the integral method introduced by von Kar- 
man into the study of laminar and turbulent boundary layers. 
In this method the differential equations are integrated across 
the stream utilizing a descriptive knowledge of the distribu- 
tion of velocity, temperature, etc. The resulting ordinary 
differential equations then describe the manner in which the 
geometric widths of assumed profiles vary along the direction 
of flow. Extension of von Kdérmdn’s integral method to 
problems of combustion simply necessitates addition of the 
integrated equations for species conservation to those for 
momentum and energy employed in the analysis of com- 
pressible boundary layer. 

Since the success of the method is directly related to the 
reasonable choice of approximate profiles, it is necessary to be 
assured that profiles selected exhibit the essential physical 
characteristics. The velocity distribution is probably the 
simplest, since after the Howarth transformation it is inde- 
pendent of temperature or composition variations. The two 
fluid streams are initially separate and it is convenient to 
describe them by separate integral relations. The streamline 
which divides the two fluids is not straight but deflects during 
the mixing process. This deflection is negligibly small, how- 
ever, and the division will be taken at the line y = 0; the 
velocity uw is constant along this streamline (9). Now let 6 
represent the extent of the mixing zone into the combustible 
mixture and —é, the extent into the hot combustion products. 
The momentum integral for the upper stream is then obtained 
by integrating the momentum equation, Equation [5], over 
the range 0 Z y Z 6, in the conventional manner (10), to 


obtain 
dx Oy \uiJy=0 


Similarly in the lower fluid region 


ll 


d 0 
dz J —& (z) w uy oy \u/y=0 


The results of these two relations must match at the dividing 
line so that the viscous shear be continuous. Since all fluid 
properties are continuous at this point, 
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The energy integrals are established through Equation 
[12]; suppose the extent of the thermal layer above the z-axis 
is 8:(a), while that below is —6.(x). It is convenient to make 
a division between the ‘‘upper” and “lower” streams along a 
line 8o(r) above which the reaction rate may be neglected. 
If the temperature of the hot stream is not too great, this may 
be taken as the isotherm having a temperature 3 = & 
corresponding to that of the dividing streamline (y = 0) in 
the absence of combustion. Then integrating over the region 
y > Bo(x) gives 


d (A: (x) 
dz JBo(z) 
d Bo (x) u v Ov 
ff -dy = — — 
dz P,U; Oy'y =8, (2) 


and for the region y < By(x) 


d Bo(xz) u 
dy— 


dz Be (x) 
d Bo(x) u v 
1 rll) = Bo (zx) 


_ 4 Bo (x) 
For these equations the matching between the two regions 
expresses the physical fact that the heat conduction is con- 
tinuous across f(z), that is 


OY B+ 


[42] 


Integration of the concentration relation [13], presents the 
additional complication that the concentration «(8) is not 


constant. Integration above the line (7) leads to the 
relation Je 


and below the line By(x) 


dx J —As (x) 


v_ 1 (x) 
where A,(z) and —A,(x) are the thicknesses of the diffusion 
zone in the upper and lower streams respectively. The con- 
tinuity of diffusive transport of combustible matter at o(z) 
gives the matching relation 
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The perturbation analysis of the previous section disclosed 
the fact that the temperature profiles have a different charac- 
ter accordingly as the value of z is less than or greater than 2; 
where a temperature maximum first appears. To illustrate 
the manner in which this fact exhibits itself in the present 
approximate treatment, consider the situation where the 
Prandtl number is less than the Schmidt number so that the 
diffusion of combustible matter into the lower stream takes 
place more slowly than the increase of the thermal boundary 
layer thickness, that is, 82 > As for the initial reaction zone. 
Clearly then, the reaction takes place most rapidly at the 
highest temperature and consequently near the boundary A» 
of the diffusion zone. In spite of the low concentration of 


Marcu-Aprin 1954 


[45] 


xe dy. 


- quently it is clear that A, = 


ouminedibi material, the temperature increases in this region 


causing a reduction in thickness f» of the lower thermal layer. 
In order that a combustion wave develop, it is necessary that 
the reaction eventually take place at a high temperature so 
that the thermal layer must decrease in thickness at a rate 
exceeding that at which the diffusion layer decreases due to 
consumption of combustible. Consequently the diffusion 
layer must exceed the thickness of the thermal layer after 
sufficient distance has been covered. As soon as this occurs, 
there exists a zone of reaction in the region where the tem- 
perature of the gas already is equal to 72 and hence the local 
temperature will exceed 7. This corresponds to the first 
appearance of a maximum in the temperature profile which 
was observed in the perturbation calculation, and conse- 
Be. at a distance x; downstream. 
Furthermore, a bifurcation occurs in the thermal layer thick- 
ness at x,;, and for x > x;; one branch lies above the maximum 
temperature, the other at the edge of the lower thermal zone, 
below the maximum. It is only the lower branch which 
possesses particular physical significance. Therefore, 
although an approximation to the velocity and concentration 


profiles may be obtained which is valid over the whole mixing 


the temperature profile must be approximated 


or downstream of the distance where a maximum is first 
obtained in the temperature profile. 

It is both convenient and reasonably accurate to represent 
the various profiles by appropriate portions of trigonometric 
functions. For example, Lock (9) has employed this approxi- 
mation in one method of treating the viscous layer between 
two parallel streams of gas. For the velocity profile the 
distribution will be 


Uo Uo 
(1—") sin y>0 
— —— }sin - 
uy uy — bo 


The concentration distribution is also relatively simple except 
that it must be written with respect to the concentration «(8») 


els 


and with respect to the limit o(x) of the combustion zone. 


Therefore 
= w(Bo) + (1 — «(8o)) sin y>Bo 


Now for the region x < x, the temperature profile may also be 


approximated in this elementary manner 


y> 


v do + = do) sin — 


do + (1 do) sin 


but for z > x; it must be possible for the profile to show a 
maximum temperature. To accomplish this in a relatively 
simple manner it is sufficiently realistic to assume that for any 
value of z, the maximum temperature occurs at the lower 
boundary of the diffusion zone, that is, where y = —A,(z). 
The maximum temperature is then just #(—A:) and the 
= + — %) sin y > Bo 


temperature profiles are 


); y < Bo... [49] 


= Jo + (#(—Az) — sin = 


2 As — Bo 
2 2 — 
2 Sy < [50] 
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Substitution of these profiles into the preceding six integral 
equations yields six simultaneous ordinary differential equa- 
tions for the functions 6,(x), A(x), Ac(x), B:(x), and 
6.(x). As an illustration it is instructive to use the same 
of example employed in the perturbation analysis, that is, where 

= Ue =u. Then only four functions are involved, 6,(x) 

7 and 62(xz) being absent since no shearing stresses are present. 
ye Furthermore the temperature Jp is just (8; + 1)/2, since it is 
that value which would occur at the interface in the absence 
of combustion. From the trigonometric relations themselves 
it is easily shown that f(x) = '/2(8: — Be) while the concen- 
tration along Bo(x) is = (Bo + Ae)/(Ar + Ae). Upon 
evaluation of the integrals in Equations [41], [42], [44], and 
[45], the following differential equations are obtained for the 
region x < 2;, after some simplification 


— (A; + = ———_ - - 
ae — 2) Seu + =) 
K(Bo)( Ar A2)I 
—2(1 —d,)ure—™ 
1 2x(Bo) 1 
——- 
(80) = — 2) Seu «(Bl — u(Bo)) (Ar + Ae)? 
ad 
(Bi Bo) = 
i 4 


de — (1 — 


where J is a definite pre over the reaction zone given by 


Dp 


Now since a numerical integration is obviously called for, it is 
necessary to investigate the behavior of the solutions to 
Equations [51] near the origin where integration must begin. 
The expansions are naturally in powers of — and may be com- 
puted as 


Actually these expressions give a fairly accurate idea of the 
behavior of the mixing and reaction zones for « < x;. For 
example, the value of x; itself, determined by the condition 
Ao(x;) = Bo(x;), may be calculated from these results by noting 
that since A, + = «(Bo)(4i + Ae) and + Bo = — Br, 


e+ 
2° — 31) 3(x — 2) 
Bo Qn q 213, Der 
u 
u 
if q 


— = 
Then from Equations [53] it follows 


the value of x; satisfies the relation 8,(2;) 
[Ai(a;) + Ac(z;)]. 
directly that 


1= ; ( ) [54] 
Cp (1 


The complete calculation is shown in Fig. 5 for this region for 
numerical values corresponding to azomethane which have 
been employed earlier. The results are essentially the same 
as those found by the perturbation analysis. 

The differential relations which describe the development of 
combustion for z > x; are more involved than those describing 
the initial region due both to the more complex temperature 
profiles and to the initial conditions which must be used at 
x = x; where integration is started. Furthermore, there are 
now five of these differential equations and they may he 
written as 


d A A ) 1 1 
™ — 2) Seu «(Bo)(1 — «(Bo)) + Az 
4 = I( As, Bo) 
bs) UT ur(1 — 
A 


+A 

1 2 

+ 
Ath) Bi — Bo 


— 2) Seu «(Bo(1 — + Ae) 


Bo) 
| A(x q — 
dx( Bo) 1 d 
dx 2(A; + As) + As) (Bi — Bo) + 
2 
Seu (A; + As)(1 — «(Bo)) 
E (A; + A2)(1 — x(Bo)) 
P (8: — Bo) 


1 


P,u (61 — Bo) 


d k(Bo)( Ar + 
~ 
A Ao) 
1 BPO 
= (55) 
where now 
— % 
Bo) = \ 42) — G -1) dg 


The development of the combustion zone into a laminar flame 
front was calculated by numerical integration of the differ- 
ential Equations [51] and [55]. The resulting value of the gas 
temperature at y = —A, is shown in Fig. 6 from the value 
x = 2; on to the point where a normal flame front 
appears. The exponential increase of the maximum tem- 
perature is shown as the gas moves away from the point of 
initial mixing. 
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OF FLOW DURING THE LATER STAGES OF FLAME DEVELOPMENT 


The detailed progress of the various regions and zones are 
traced in Fig. 7, in which the vertical scale is exaggerated, and 
demonstrates well the relatively sudden appearance of the 
laminar flame propagating into combustible mixture. This 
is indicated by the rapid convergence of the lines representing 
different temperatures and hence the appearance of strong 
temperature gradients. Also of interest is the observation 
both from Fig. 6 and from Fig. 7) that the small perturbation 
analysis is probably adequate to describe the major portion 
of the reaction zone preceding the appearance of a lamin 
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flame front. The temperature distribution profiles which 
occur near the end of the ignition zone, plotted in Fig. 8, 
show the rapidity with which the temperature gradient 
steepens to the laminar flame front. It is particularly striking 
that the major steepening of the gradient takes place in the 
last 10 per cent of the reaction zone. 
Concluding Remarks 

Through application of a boundary layer type of analysis to 
the problem of ignition and combustion in the mixing zone 
between parallel streams of combustible gas and combustion 
products, some of the essential features have been deduced 
with relative ease. In particular, although the stream of 
combustible gas eventually ignites, the distance (or time) 
required is excessive where the temperature of the hot stream 
is too low. Then the flame develops so far downstream that 
it is essentially “blown off” any finite apparatus. 

In a broader sense the analysis shows that a new class of 
combustion problems is open to investigation through ex- 
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to the present one, this class includes such important prob- 
lems as the theory of thermal quenching near a cool wall, the 
ignition and flame stabilization on a heated surface from 
which the boundary layer is unseparated, the erosive burning 
of solid propellent grains, and a great many others. There is 
no essential difficulty in proceeding to problems with axial 
symmetry and it seems quite possible that the process may be 
extended to include cases of turbulent mixing. Thus the 
theory of the plane laminar flame together with the descrip- 
tion of such boundary layer regions as may be treated by the 
methods just described, seem to provide an adequate descrip- 


tion of a relatively large group of problems which involve 
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Ramjet Diffusers at Supersonic Speeds 


(Continued from page 84) 


The preceding examples have assumed fortuitous and some- 
times improbable combinations of conditions to promote os- 
cillations, but, as missile experts all know, Mother Nature has 
a perverse habit of assembling just such circumstances to 
confound the luckless experimenter. In the foregoing we have 
considered only simple configurations. The possibilities be- 
come almost unlimited when we introduce additional mecha- 
nisms such as fuel and/or area controls, each with its own 
lag and inertia characteristics. Then truly, Missiles Engi- 
neer, “there are more things in heaven and earth than are 
dreamt of in your philosophy.” 
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Low-Speed Combustion Aerodynamics 
: ROBERT A. GROSS! and ROBIN ESCH? 7 
Combustion Aerodynamics Laboratory, Harvard University, Cambridge, Mass. =" 
The problem studied is the two-dimensional, low-speed, Introduction 
inviscid, steady flow field containing a finite length dis- ‘ 
continuity which represents a flame. The flow takes place HE problem of determining the flow field associated 
in unbounded space and has uniform parallel flow at in- with the burning of combustible fluids is an extremely 
finity. It is shown that if the density ratio \, the pressure complicated task. There appears to be little hope at pres- 
ratio p/p: across the flamelike discontinuity, and Q, the ent of solving the entire real problem which encompasses 
heat of reaction are constant, the flow field is everywhere the realm of compressible viscous fluid dynamics, chemical 
irrotational. For such cases the flow field can be described kinetics, heat and mass transport, etc. To obtain a better 
by classical potential theory where the flame discontinuity understanding we desire solutions of simplified problems 
is represented by a source sheet. The fundamental equa- based upon the real phenomena. To be useful, these prob- 
tions for the flow field and the discontinuity are formu- lems must be simple enough to permit analytical treatment 
lated for low Mach number flows. The source sheet equa- yet must at least contain the essential features of the phe- 
tions are given in detail. Two approximate solutions to nomena, 
flow fields are presented; one is a finite length straight- For many investigations in the field of aerotkermodynamics 
line flame inclined at an angle to the free stream; the it is reasonable to assume that the chemical reactions take 
second is the V-flame shape. Dimensional considera- place in such a short time and space that this region may be 
tions are discussed and some experimental results are re- treated as a discontinuity in the flow field. The thickness 
ported. of such reaction zones is of the order of 1 mm at atmospheric 
r pressure. On either side of this discontinuity the flow may 
a be considered as laminar, inviscid, and steady. Assuming 
Nomenclature the existence of a normal flame velocity permits us to com- 
pute and describe by known methods the flow phenomena 
a = local speed of sound occurring in many combustion problems. 
f(r),g(7) = parametric quantities for flame shape In particular we shall consider here the case of two-dimen- 
ho = stagnation enthalpy = h + (q?/2) sional, low-speed flows. By low speed we mean that the 
i,) = unit vectors in the z and y directions Mach number is so small that compressibility effects may be 
l = coordinate along source sheet ignored. Thus, the reactants (upstream of the discontinuity) 
L = length of source sheet Te and the products (downstream of the discontinuity) are con- 
M = flame-speed Mach number = S/V yRT = S/a sidered to be of constant, though different, densities. The 
" = unit vector normal to flame discontinuity fluid will be considered as an ideal gas. 
p = pressure 
Veo = velocity at infinity wl 
Q = heat of reaction from combustion ; aie The flow field may be considered as two regions separated 
r = distance from a point in flow field to element of the by the flame discontinuity. We seek possible fluid motions 
R in each region such that these regions are connected by con- 
servation of mass and momentum across the discontinuity. 
Two equations express the principles of conservation of 
. = zcomponent of velocity 7 mass and momentum. Under the simplifying assumptions 
a" = ycomponent of velocity previously mentioned, these equations are 
= coordinates 
= complex variable 1] 
Y = ratio of specific heats ~ Or oy 
r = source sheet strength per unit length of discontinuity — ' 
v = gradient operator 
6 = angle between flame discontinuity and 7, =¥ 
d = density ratio = p;/p2 where q, p, and pare the velocity, density, and pressure of the 
M = dipole strength =? gas at the point x,y, and u,v are the velocity components of @ 
4 a density ro in the x and y directions, respectively, and V represents the 
Fig.7 gradient operator. Equations [1] and [2] must be satisfied 
& = real part of complex potential everywhere in the reactant and product flow fields. 
¥ = stream fenction The condition to be satisfied across the discontinuity are 
that mass and momentum be conserved. Mass conservation 
yields 
_ Presented at the 8th ARS National Convention, New York, - 
N. Y., December 3, 1953. 
n = normalcomponent = = =~ a 1 Research Associate. Now associated with Fairchild Engine 
t = tangential component and Aircraft Corporation, Farmingdale, L. I., N. Y. 
8 = source sheet 2 Research Assistant. 
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where gz, is the velocity component normal to the disconti- 
nuity on the product side, S is the flame speed or normal ve- 
locity component on the reactant side, and X is the density 


ratio of reactants to products eure 


the discontinuity, respectively. 
Momentum conservation equations tangential and normal 
to the discontinuity are 


The kinematics of the discontinuity, expressed by Equations 
[3] and [5], are shown in Fig. 1. Defining the flame speed 
Mach number as 
De 
—=1- 1) M,2 


Nee 

Equation [6] may be written as 


where 7 is the ratio of specific heats. Tsien (1)? has set down 
the analogous equations for the compressible case, where one 
needs the addition of the energy conservation equation. The 
fact that for low-speed flows one may consider the flow as 
incompressible and the density ratio \ as a constant leads to 
considerable simplifications. 

Equations [1] and [2] plus the boundary conditions specify 
the reactant and product flow fields, while Equations [3], 
[5], and [7] connect these two fields across any flamelike dis- 
continuity. 

The reactant flow fields considered here will have uniform 
and paralle] streamlines far upstream of the discontinuity. 
Therefore, the reactant flow field is irrotational everywhere. 
This follows from the theorem of Lagrange and the fact that 
there is a unique pressure density relationship. Thus 


Equations [1] and [2] are satisfied for the irrotational reactant 
flow field if one finds a solution to Laplace’s equation 


Since the flame discontinuity can generate vorticity, it is 


necessary to examine what happens to w as the discontinuity 
is crossed. Tsien (1) has shown that in general, the combus- 
tion products are rotational.4 To see why this is so, we make 
use of a theorem by Crocco (2). For steady inviscid flows, 
the equations of motion may be written as 


or for two-dimensional flow (3) 
dho\ 


3’ Numbers in parentheses refer to the References on page 101. 
‘It should be noted that the » defined by Tsien is different 
than that used in this paper. 


FLAME - LIKE 


FIG. 1 FLAME FRONT KINEMATICS 


where ho is the stagnation enthalpy, 7’ the temperature, and 
o the entropy. On the reactant flow side, w, = 0 and there 
are no entropy or total enthalpy gradients. The energy 
equation which relates the product and reactant flow field is 


where Q is the thermal energy introduced into the flow at the 

discontinuity by means of the chemical reaction. This heat, 

released by combustion Q, may be considered as a constant. 
The entropy change across the discontinuity can readily be 

computed 

C,In 


— = C,In 
p 


Pr 


The entropy change across the discontinuity is a function of 
the density and pressure ratios at that point. Thus, if ) is 
not a constant (which is the case in compressible flow), or if 
the flame speed S, or the heat of reaction Q is a variable along 
the flame front, the product flow is indeed rotational. This 
is the case whenever there are entropy or energy gradients 
in the flow field. The product flow field is therefore deter- 
mined by 


V? 


where the flame boundary values of @ are given by Equations 
[10], [11], and [12]. 

It must be noted that the stream function y in Equation 
[9] is different from that in Equation [13], since the density 
of the reactants is different from the products. The solution 
of Equations [9] and [13] coupled by the flame front Equa- 
tions [3], [5], and [7] plus boundary conditions, is quite 
difficult. However, several special cases of interest have al- 
ready been treated. Scurlock (4) has numerically solved the 
case of the flame shape and flow field in a two-dimensional 
combustion chamber of constant width. He neglects the 
pressure drop across the flame front, treats \ as a constant, 
and reduces the problem to a quasi-one-dimensional calcu- 
lation. Tsien (1) examined the same problem and showed 
that, with these same assumptions, the calculations could 
be made very much simpler. Tsien also examines the com- 
pressible case and treats the vorticity production by the 
flame front in detail. Ball (5) solved the same problem, in- 
cluding the pressure change at the flame front, by the relaxa- 
tion technique. 
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The Problem ~_| 


Much of the present-day low-speed combustion research is 
Jone on burners in which the gas is flowing into essentially 
ynbounded space. This combustible jet is often surrounded 
by moving secondary air to stabilize the flow. For example, 
at the Harvard Combustion Aerodynamics Laboratory, 
fames are produced in the combustible jet at a distance well 
ibove the nozzle in the combustion tunnel (see Fig. 2). The 
fame is thus burning in what is nearly equivalent to un- 
bounded space. The primary (combustible) jet is sur- 
rounded by moving secondary air. The flow is nearly two- 
dimensional and steady. What we desire to find are analyti- 
cal solutions to flow fields of the type which we observe under 
these experimental conditions. 


vIG. 2. V-FLAME IN HARVARD COMBUSTION T 


The problem thus consists of a two-dimensional, steady, in- 
viscid flow in unbounded space. The flow field contains a 
combustible stream surrounded by secondary air. In actu- 
‘lity, the flame is stabilized on a thin wire or flat plate above 
the primary nozzle. The flow is low speed, with velocities 
well below 100 ft/sec. The reactants leave the nozzle with 
a uniform flat velocity distribution. Thus, Equation [9] 
determines the reactant flow field. Since the reactant flow 
field is irrotational, the vorticity of the product flow field is 
determined by Crocco’s theorem expressed by Equation 
10a]. Since there are no entropy gradients in the reactant 
fow, Equations [12] and [7] permit us to find the entropy 
gradients on the product side of the flame discontinuity. 
Equation [11] gives the stagnation enthalpy on the product 
ide of the discontinuity. Therefore, Equation [10a] permits 
us to compute the vorticity on the product side of the flame 
front. However, we shall consider \ and Q as constant. If 
the flame speed Mach number is also a constant (implying 
constant flame speed and constant 7) then the pressure ratio, 
h/p; is a constant. These three assumptions, that \, Q and 
Mf, are constant, are good ones for low-speed flows and they 
lecessarily imply from Equation [10a] that the product flow 
feld is irrotational. Therefore, the entire flow field is irro- 
tational and the powerful tools of classical potential theory 
tay be used. 

It is interesting to note that one cannot neglect the pressure 
thange across the discontinuity and still satisfy the mass 
conservation Equation [3]. The assumption of p2/p; = 1 is 
quite good for high-speed flows where @; ~ 92 but cannot be 
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used in low-speed flows. This small pressure change across 
the flame front arouses some interesting questions concerning 
the stability of such flamelike discontinuities, but we shall not 
treat such aspects here (6). 

With \, Q, and p2/p; constant, the entire flow field is irro- 
tational. The different density of the reactants and products 
still requires the use of separate ¥ functions. However, if 
one used a single ¥ function for both reactant and product 
fields, some interesting progress toward an approximate solu- 
tion may be made. Thus, the entire flow field will be de- 
scribed by a single ¥ function and the influence of the flame 
discontinuity will be exactly accounted for by an appropriate 
source sheet. By considering the reactants and products to 
be of different densities, the use of a single ¥ function will 
introduce a difficulty along the two free streamlines which 
separate the products from the secondary air. However, this 
simplification permits us to arrive at an exact solution to the 
flow field everywhere except along these two free streamlines. 
Although this difficulty may be serious, this simplification 
does permit us to arrive at a solution which describes the type 
of aerodynamic influence that the flame exerts. 

Thus the flow field is a solution to Equation [9] with uni- 
form parallel flow at infinity. The flamelike discontinuity 
will be represented by an appropriate source sheet. An ex- 
amination of Fig. 1 or Equations [3] and [5] will show that by 
choosing the source strength per unit length of the discon- 
tinuity to be 


(A — 1)S 
2r 


the kinematic conditions across a flame front are satisfied. 
Each element of the source sheet will influence the remainder 
of the flow field. The solution of Laplace’s Equation [9] in 
two dimensions, which gives the velocity induced by an ele- 
ment of the source sheet at a distance r, is 

dl (x—1)Sdl 

T 


r 2r r 


Consider Fig. 3.. An arbitrary shaped discontinuity is shown 
and it is specified by the parametric equations x = f(r) and 
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y = g(r). The flow field induced by the element of the 
source sheet di at the point (¢,7) is given by 


r2 


ic 
a 
( 
if 
f(T) 
y= g(t) 
ry 
if 
r 
‘ 


where 
dl = + [18] 
The primes indicate differentiation with respect to r. The 


velocity induced by the entire source sheet at the point 
(t,n) is therefore given by 


X 
2r 1 r2 


— 5) + j(g — 


where 7 and 7 are the unit vectors in the z and y directions, 
respectively. Equation [19] gives the contribution to the 
velocity from the source sheet only. This induced velocity 
approaches zero as r becomes very large. It is desired to 
have uniform parallel flow far from the discontinuity. Since 
Laplace’s equation is linear, solutions may be added. The 
only boundary conditions to be satisfied are those at the dis- 
continuity and at infinity; thus, addition of a uniform ve- 
locity is permissible. The complete velocity field is given by 


(5,9) 
[19] 


where 7. is the velocity at infinity. 

The fact that the velocity component normal to the dis- 
continuity on the reactant side must be equal to a pre- 
assigned quantity (the flame speed) has not yet been taken 
into account. Therefore, a constraint equation which limits 
the shapes of permissible flamelike discontinuities is given by 


=S........ [21] 


where (20,0) represents a point on the discontinuity surface, 
and n is the unit vector normal to this surface. That is, 


—tg"( 10) + 
+ 


Designating the velocity at infinity as 


~ 
the desired permissible flame shape constraint equation can 


be written out in detail. It is 
s — 90'Gz —(A — 1) 
(fo’? + go’?) 10) 2r 
n—>g(70) 


+ g(r)? 


(g(r) — ant . .[23] 


This is a nonlinear integrodifferential equation where the in- 
tegral must take on preassigned values along the unknown 
function—namely, the shape of the discontinuity. The 
uniqueness of the problem is uncertain and appears to require 
additional physical conditions to be specified, such as stability, 
etc. 


Some Solutions 


There are no known general techniques for finding a solution 
to equations such as [23]. However, one exact solution has 


been found. Itis 
r=ar y=br 


where a and b are constants. This solution represents a 
finite straight line inclined at a fixed angle relative to the ve- 
locity vector at infinity. By choosing the coordinate sys- 


tem so that q,, lies along the y-axis, the solution may be veri- 


98 


fied by substitution into Equation [23]. Thus — 
1-1 (a? + 
2 

lim an — 

tare (ar — $)? + (br — 

n—>bro 


Upon evaluating the integral and then taking the limit, the 
right side of the above equation approaches a value of —r. 
Therefore 


(a? + 2 


where @ is the angle between the discontinuity and the ve- 
locity vector at infinity. This same result may be obtained 
simply from the geometry at the flame front; see Fig. 8. 
This gives the straight-line flame-front discontinuity position 
in terms of the dimensionless quantities \ and S/q.. 

Flame fronts that are extremely close to a straight line have 
been observed and, in fact, can be produced at will in the 
Harvard Combustion Tunnel. The details concerned with 
finding the equation of the streamlines for such a flame front 
are treated in Appendix 1. The result of this analysis gives 


the source sheet contribution as 


where ¢ and 7 are coordinate directions, and ¢ and 2 are 
angles as specified in Appendix 1. The — for the 


streamlines in the flow field is 
= (x cos @ + y sin 6) cosa) + 
(x sin @ — y cos @) - — fe 


(x, cos 6 + y: sin 6) log + (zz sin — ye cos — 
(x; sin — y; COs . . [27] 


v is the stream function, and the subscripts 1 and 2 indicate 
the end points of the source sheet. Since lines of constant 
are streamlines, Equation [27] permits mapping the entire 
flow field by specifying the values of S/q,. and A. Such a 
flow field has been mapped out and is shown in Fig. 4. _ This 
represents a flow for S/q.. = 0.1, and \ = 7; values which are 
common for deflagration. 

At the discontinuity, the tangential component of velocity 
gr is small compared with usual values of g, except for points 
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close to the ends of the source sheet. Here, the solution has 
singular points and the values of g become very large. This 
of course violates the low Mach number restriction and the 
solution is not valid in the near neighborhood of the end points 


of the discontinuity. The pressure distribution is given by 


since the stagnation enthalpy is constant in either the reactant 
or product flow field. However, the value of ho changes across 
the flame discontinuity as indicated by Equation [11]. By 
the nature of the solution, 7 is continuous across the free 
streamline which bounds the products and the secondary air. 
However, /, > ho,, and this results in the difficulty previously 
mentioned, namely, that these two bounding streamlines 
have a resulting pressure drop across them and hence could ff: 
not be steady. This results in an error in the position of the 
streamlines. The magnitude of this error will be compared 


with experiment in a later section. 
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FIG. 4 STRAIGHT-LINE FLAME FLOW FIELD 


Examining Fig. 4 shows some interesting features. The 
streamlines which pass through the ends of the source sheet 
are considered to be boundaries of the combustible core. The 
fow outside this core is assumed to be noncombustible second- 
ay air. The effect of the flame front on the remainder of 
the flow field is evident. The entire combustible core is de- 
fected from its initial path because of the presence of the 
| tame front. Such stream deflections are evident in schlieren 
photographs of real flames. One should recall that in this 
malysis these effects are not because of such things as heat 
transfer, but are the aerodynamic interaction of the flame 
ont and the flow field. 

A commonly observed flame is the so-called V-flame. Such 
iflame is shown in Fig. 2. Since this appeared so similar to 
the straight-line solution already given, it was hoped that 
ve could generate simply a solution to this V-flame by reflect- 
ing the first solution about its y-axis. Although the resulting 
fow field still satisfies Equations [3], [5], and [9], the flame 
weed constraint Equation [23] is violated because of the 
nutual interaction of the two sheets. 


: An examination of the induced normal velocity to these 
heets indicated that q:, varies about as /? where / is the length 
tlong the discontinuity, measured from the vertex. By plac- 

' ig a proper strength dipole at the origin, the flame speed 


may be made nearly constant over the major length of the V. 
This problem of the V has been worked out and is shown in 
it Hig. 5. The area within the circle at the vertex contains the 
lomain where the singularity (resulting from the dipole) 
iolates the basic assumptions of the problem. The details 
{the V solution and dipole addition are given in Appendix 2. 
A graph of qi, as a function of / is shown in Fig. 6. 

In comparing the V solution to the single sheet, it is evident 
at for the same values of \ and 6, the V requires a smaller 
alue of S/q... This is physically evident from the mutual 
iteraction of the two flame surfaces. The V solution has 
“arly a constant flame speed except close to the apex where 
he solution breaks down ordinarily. 
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FIG. 5 V-FLAME FLOW FIELD 


Dimensional Considerations; Some Experimental 
Results 


The previous analysis permits an interesting insight into 
the important physical parameters present in such aerother- — 
modynamic problems. For those cases where the flame may | 
be treated as a discontinuity and no solid boundaries are | 
present, the following dimensionless groups determine the | 
problem: 

8/d. = flame speed/velocity at infinity a 
pi/p2 = X = density ratio 
M, = s/a = flame speed Mach number 
Q/cpT; = heat of reaction/enthalpy of reactants 


y = ratio of specific heats = cp/cy 
In the presence of boundaries it would appear that the vis- 
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cosity and thermal conductivity of the gases must be in- 
cluded. 

In the development of the present problems, some inter- 
esting questions arose. In reality, the flame speed Mach 
number M, varies slightly and the product flow is rotational. 
If the flow is slightly rotational, when can this rotation be 
justifiably neglected? The characteristic length present in 
this problem is the finite length of the flame front, L. It 
would appear from dimensional considerations that the vor- 
ticity does not appreciably influence the shape of the flame 


front 
In the combustion tunnel at Harvard University it is pos- 
sible to experimentally check some of the conclusions drawn 
from this analysis. In the presence of a thin wire flame holder 
held above the primary nozzle, flat sheets and V-type flames 
can easily be produced. The easiest quantitative prediction 
to check is that of Equation [24], the flame angle. This was 
done for propane-air flames. Some typical results are given 


below. 
0.0490 10° + 2° 
7.23 0.0494 10° + 2° 
7.00 0.100 — 23° 16° + 2° 
7.30 0.0962 23° 15° + 2° 


The above numbers represent values of s and q., over the range 
1.15 < S < 1.47 ft/sec, and 13.8 < q. < 30 fps. The 
agreement appears good for small values of s/q.. but is poor at 
higher values. 

The neglect of the product vorticity does not appear to 
adequately explain the disagreement. However, the use of 
a single ¥ function for the products and the secondary air 
belies the fact that physically they have different densities. 
Thus, in our solution if we consider the density of the prod- 
ucts and secondary air to be different, we have violated mo- 
mentum conservation along the bounding streamlines. This 
may partially explain the difference in angles measured ex- 
perimentally and predicted from our equations. 

If one neglects the aerodynamic effect of the flame, the 
value of @ would be given by sin 6 = s/q,. Using the above 
data, this would predict values of 6 less than 6°, which are 


far too small. 
@ 
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he APPENDIX 1 
‘ Field Equations for Straight Line Source Sheet 


We consider an element of length di of the source sheet, as 
shown for the general case in Fig. 3. The sheet extends in- 
finitely in the direction normal to the page. The potential 
arising from such an element is 


generating velocity 
dv = V(d®) = 


The entire potential and resulting velocity are obtained by 
integrating these expressions over the entire source sheet, 
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and superimposing a constant flow field. For any flame 
shape, the resulting field satisfies the laws of mechanics both 
throughout the reactant and product regions and across the 
discontinuity. The one remaining condition, that the nor- 
mal velocity of the reactants into the source sheet be equa! 
everywhere to the flame speed, has been shown previously to 
be satisfied in particular for a finite straight line source sheet 

The integration of the above equations for this case is 
straightforward; it may be checked that Equations [25| 
and [26] result. A more convenient technique is the use of 
the complex potential 


W(z) = 


since it gives also the stream function y, by definition re. 
lated to the potential @ by the Cauchy-Riemann conditions 


her 
oy ox 


which, of course, imply that y is a solution of V?y = 0. Since 
Ve-Vy =0, the lines ¥ = constant are everywhere norma! 
to the equipotential lines, hence being streamlines. 
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FIG. 7 SOURCE SHEET 


The two-dimensional complex potential for a source element 
of real length dl is 
dW = T log (z — 2’)dl 


(this is the analytic function having d@ as its real part, where 
= The entire flow field is then described by 


W(z) = hes (x 
i 
where the configuration is as shown in Fig. 7; the integration 
is taken in the real direction. Hence we have ; 
2: 
W(z) = — 2’) log (z — 2’) +(z- “| 
21 
(z — a) 
sis + at #2 + log (2 24) — 21 log (2 2) 
ting 
As shown in Fig. 7, we define z = ¢ + in 
z-a4=ne™ 
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giving the real and imaginary parts of W (z) 


r 
Flog” — — +h — + log — 
r2 


log ri + — 2) 


r re 
+ — +m log + S262 — 
r re 


where we use nz = m1. 
By differentiation of either of these (expressing 41, $2, "1, 
and rs in terms of ¢ and »), we find 


=Tl 25 
= — [26] 


and finally, by a counterclockwise rotation of the coordinate 
system by (2/2 — 6) (see Fig. 8), and by adding —2xq.,, the 
stream function of a constant velocity in the y direction, we 
obtain the final stream function 


#8 

q 

= (x cos + y sin (108 + 
re 


r 2 
qo. 
(x sin 6 — y cos @) | — — sin | — 


r 
(x; cos 6 + y sin 6) log + (2 sin 8 — ye cos 8) d: — 
Te 
(x; sin — y; Cos . . [27] 


We may check that this solution has the properties claimed ; 
thus, examining Equations [25] and [26], we see that the in- 
duced velocity goes to zero as r grows large in any direction. 
Further, the normal velocity q, is constant along the source 
sheet, and discontinuous by the constant amount (A — 1) S = 
27T across the source sheet, as desired. 


APPENDIX 2 
V-Sheet Equations 

Making use of the line of symmetry, we put the origin at 
the vertex of the V. Let ¥(2,y) be the stream function for the 
single sheet (forming one leg of our V), as derived in Appendix 
1. Then —y¥(—2z,y) will give the stream function contribu- 
tion at the point (z,y) from the other half of the V. The 
minus sign arises as, in reflecting about the y-axis, the co- 
ordinate system becomes left-handed. Then, in order that 


qz = Oy/dy 
qy = — Oy/dxr 


still be true in the original coordinate system, we must replace 
¥(—a,y) by —¥(—2,y). 
Hence, proceeding as before, our complete stream function 
is 
= Wry) — — 


This of course satisfies the laws of mechanics throughout the 
plane and across the discontinuity; however, it gives a: qin 
(normal velocity of the reactants at the source sheet) which 
differs badly from a constant (see Fig. 6). This discrepancy 
is not remedied by changing @ or our fundamental parameters 
\ and s/qg,. Thus we are led to the conclusion that, for our 
simplified model, the V is not a possible flame shape. 

The fact remains, however, that flames which are fairly 
two-dimensional and fairly V-shaped are observed. (See 
Fig. 2.) Though conditions in the Harvard Combustion 
Tunnel are not the ideal conditions assumed here (the second- 
ary air region does not extend to infinity, the problem is not 
exactly two-dimensional, a flame holder and nozzle boundary 
layer have been introduced, etc.), it is felt that experimental 
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é 
FIG. 8 VELOCITY CONTRIBUTIONS AT FLAME FRONT 


perturbations do not explain the disagreement. Perhaps 


- some assumption (such as that of constant flame speed) is 


unjustified, or perhaps the influence of a flame holder is impor- 
tant throughout the real flow field. However, as remarked 
above, the use of a single ¥ function is unrealistic, and proba- 
bly this is responsible. 

Putting aside for now the problem of an exact solution of 
the simplified problem, we seek an approximate solution to 
the V-flame which does not involve an exorbitant amount 
of computation. Let us add a dipole at the origin, oriented 
with source in the negative y direction, which induces a ve- 
locity field with stream function 


where u is a constant. (It may be readily verified that this 
‘isa solution to V?y = 0.) The total stream function then is 


= (2,y) — W(—2,y) — (do + + + y?) 


The addition of the dipole and Aq. has introduced a change 
iN in given by 
sin @ + sin 


Choosing arbitrarily a tolerance of 25 per cent, it was found 
that by adjustment of the parameters u and Aq.., gin could be 
made approximately constant over the upper three fourths of 
the source sheet. (See Fig. 6). The resulting approximate 
solution is presented in Fig. 5. Within the circle about the 
origin the dipole dominates the solution. The effects of the 
dipole are noticeable just outside of the circle; fortunately 
they fall off as the square of the distance from the origin. 
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“High. Frequency Co Combustion 


= 

A theory of unstable high-frequency oscillations in 
solid propellant rockets is advanced with the mechanism 
of self-excitation based on the following simplified model. 
Both the rate of primary decomposition of solid propellant 
element and the rate of activation of the intermediate prod- 
ucts during the time lag are assumed to depend upon 
the gas oscillations. Both interactions are expressed in 
terms of the instantaneous pressure according to power law 
with exponent n and m, respectively. Thedependenceof the 
rate of decomposition on drifting velocity is shown to be 
not of fundamental importance in determining the sta- 
bility of oscillations in the combustion chamber. Analysis 
shows that a simplified over-all pressure index of inter- 
action S = m — (n/2) must be bigger than certain mini- 
mum value if unstable oscillation is to be possible. Un- 
stable oscillation further requires the time lag to be in 
proper range. The most unstable mode is shown to be 
the fundamental spiral mode, and its tangential com- 
ponent with exactly the acoustical frequency is the one 
which is most easily observed at finite magnitudes because 
of its rapid amplification. Several configurations of the 
propellant grain with radial burning surfaces have been 
investigated. Rod grain becomes more stabie while other 
grain shapes become less stable in the course of firing. 
Rod grain is the most stable while tubular grain is the 
most unstable configuration. The stabilizing effect of 
inserting a nonburning solid rod into tubular grain is 
verified by the theory. Several other aspects have also 
been treated. Solid rocket with end burning grain also 
becomes less stable in the course of operation. Both the 
spiral and the longitudinal modes in rocket with end- 
burning grain are briefly investigated. 

reve 


Rockets with Tubular Grain 


allele first a solid propellant rocket using perforated 
cylindrical grain with interior burning surface, that is, the 
tubular grain. In this case there is no solid nonburning 
surface. The radial boundary condition [14] is replaced by 
the requirement that all disturbances must remain finite at 
the axis of the cylinder. Thus, the constant c in Equation 
[8] must vanish. When the radius of the burning surface is 
taken as the reference length, it can be found either directly 
or by taking the limit of r; ~ 0 with r2 = 1 in Equation [20], 
that 


if! 


= 


For the fundamental spiral mode, 8 is equal to unity. 
Then mo = 1.84 which is the lowest zero of J1’(n). Equation 
[27] gives x = 0.70 which is a constant throughout the entire 
period of firing of the rocket. Equation [25] for the mini- 


5 Part I of this paper appeared in the January-February issue 
of JET PROPULSION, pages 27-32 


Rockets. 
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“Instability in in 
Part 2 = 


mum index of interaction compatible with unstable oscilla- 
tions then becomes 
9( ?) 
= [28] 


1 
Suin = 27 E + 0. 70 


During the course of firing of the solid rocket after steady 
state has been reached, the burning rate at stable operation 
is constant and, therefore, the normal velocity 0, of the 
burned gas at the burning surface is also constant. For the 
interior burning surface, 9, is negative. When the propellant 
burns, the radius of the burning surface is increasing. As 
the length of the grain is constant, the value of \ is decreasing 
with time. Whether the value of Sin is decreasing or in- 
creasing with decreasing 2 is of particular interest because 
this indicates whether the system becomes more stable or 
more unstable in the course of firing. 

By expanding the left-hand side of Equation [16] into real 
and imaginary parts, one obtains 


tan Az] = 
1 dé, Sin 2 AE, — AE; sinh 2 AE; 


= 
2 cosh? cos? AE, + sinh? sin? AE, 
9[AE tan AE] = 
1 AE; sin 2 AE, + AE, sinh 2 AE; = 
2 cosh? dé; cos? At, + sinh? AE; sin? AE, 


tan Ag] 
9[AE tan AE] 


The values of J[Aétan AE] and R[Aé tan AE]/ J[AE tan AE] 
for several values of \é, and dé; are plotted in Fig. 3.6 Equa- 
tion [30] shows that the value of ®[Aé tan AE]/ J[AE tan AE} 
is always negative with its magnitude equal to A;,/A, as 
determined from Fig. 2 by the value of the reduced frequency 
parameter Q = Xé,.e/[(2/y + — for the axial os- 
cillation. For the fundamental mode of axial oscillation, 
\£, is of the order of unity and the reduced frequency parame- 
ter is of the order of e which is in practical cases a small 
fractional value. Consequently A;/A, is of the order of 3 or 
even larger. Under such circumstances, Fig. 3 shows that 
the values of ®[AE tan AE] and J[Aé tan AE] can be repre- 
sented by their asymptotic values — dé; and \é,, respectively, 
with reasonable accuracy. Hence, using Equations [29], one 


9(£?) ie? 

Ue 


Since viscous effect is neglected, 7%, represents the mean axial 
velocity at the exit x = \, and the continuity of mass flow 


gives 
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6 Figs. 1 through 5 appeared in Part 1. sinc 


Equation [28] for the minimum value of S of rockets with 
tubular grain, where r; = 0 and rz = 1, becomes 


1 
Smin = — [1 + [33] 
2y 


where 9 = 1.84 for 8 = 1 has been substituted. 

For a given solid rocket e, |é,|, and \ are known. The ratio 
of the reduced frequency parameter 2 and A, can be directly 
calculated from 


= — = 4, 
Ay eL\y + :) ] 
The values of A, and A; can thus be read from Fig. 2, and 
Swix for the particular configuration can be calculated from 


Equation [33]. For example, with 1/e = 7.5, |i:| = 0.01, 
and \ = 15 we have for y = 1.20 


A, = 0.14 A; = 0.30 


Smin = 0.62 


i, = 2d|6,| = 0.30  /A, = 1.78 
2 ySmin — 1 = 0.49 


This calculation is repeated for different values of \ of two 
rockets with 1/e = 7.5 and 10, respectively, using a solid 
propellant with |d,| = 0.01. The results are plotted in Fig. 4. 
In the determination of A, from the ratio of 2/A,, it is found 
that in certain range of values of the ratio Q/A,, A, is multi- 
valued. The lowest value is selected because we are only 
interested in the fundamental axial mode which gives the 
lowest value of Sin for the time being. If \ is kept increasing, 
both the value of # = 2 |d,| and Q/A, increase monotoni- 
cally. When the value of Q/A, is sufficiently large, no solu- 
tions of A, can be found from Fig. 2. Solution, however, 
does exist in the region of 2 beyond that covered by Fig. 2, 
but this solution will give an Smin much bigger than unity. 
Furthermore, the approximate method is no longer valid in 
this region. Therefore, the solutions in regions of large Q 
are discarded. The results of calculation as given in Fig. 4, 
indicate a very steep rise of Smin with increasing \ before 
\ is big enough to invalidate the present method of approxi- 
mate calculation. The solution for Smin with larger values 
of \ is hence of little practical interest. It is, therefore, in- 
ferred that during the course of firing of a solid rocket using 
tubular grain, the solid rocket destabilizes at a very fast rate 
in the early stage of operation until the value of » of the 
rocket becomes smaller than certain critical ranges of values 
which depends on the axial geometrical parameter 1/e 
(A = 15 to 18 as shown in Fig. 4 for the cases calculated). 
After passing through this critical range, the rate of de- 
stabilizing becomes much slower. The small high-frequency 
disturbances which are stable at the earlier stage of operation 
may thus become unstable at a later stage when the value 
of \ of the system becomes sufficiently small as to give a 
Smin lower than the value of the index S of interaction of the 
particular propellant used. This destabilizing tendency of a 
solid rocket with tubular grain agrees with the fact that 
abnormal pressure peak occurs in the later period of opera- 
tion of a solid rocket preceded by a period of smooth operation. 

From Fig. 4, the importance of the axial geometry of the 
solid rocket is easily recognized. The destabilizing effect of 
using shorter “effective nozzle,’’ i.e., larger value of 1/e; is 
very significant. 

Suppose an end seal of radius 7 is installed at the end of 
the propellant grain; then Equations [32] and [33] become 


le = 2 34 
Ue = ro? = Ala, ] 
and 
1 
Smin = E + 5.15 [35] 
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since 79 <1, # is greater than 2Ald,|. For the same value 


of \ the product A;A, would change slightly from the corre- 
sponding value without end seal. However, this change of 
A,A, is rather small as compared to the increase of S min, due 
to the factor 1/7.4 The use of end seal with rp < 1 is thus 
stabilizing, and the stabilizing effect is larger when 7 is 
smaller. Some calculated results are shown in Fig. 5. This 
stabilizing effect of using end seals of smaller inside diameter 
as obtained from the present investigation agrees with the 
experimental results in (1). They find rough burning in 
rockets using end seal of larger inside diameter, but find 


- smooth burning in rockets using end seal of smaller inside 


diameter. 

Now if the end seal is nonerosible, that is, the inside 
radius of the end seal is constant throughout the operation, 
the value of rp decreases continuously with time during opera- 
tion because the burning surface radius or the reference length 
increases with time. Hence, the stabilizing effect of the non- 
erosible end seal increases with time. Whether this sta- 
bilizing effect can offset the destabilizing effect of decreasing 
\ toward the later stage of operation can be determined by 
further calculations for the particular geometrical configura- 
tion. If the end seal is erosible the value of ro decreases less 
fast than that for nonerosible end seal and, in extreme cases, 
ry May even increase depending upon the rate of erosion 
of the end seal. The rocket with nonerosible end seal is 
thus more stable than a similar rocket with erosible end 
seal of the same initial inside diameter. 

Now we investigate the stability of the fundamental non- 
spiral mode in rockets with tubular grain as compared to the 
stability of the fundamental spiral mode. For nonspiral 
mode 8 = 0, and the first or the lowest nonzero value of 
that makes Jo’(m) = 0 is 3.83 and x = 1. The expression 
for Smin for the fundamental nonspiral mode becomes 


1 
Smin = [1 + [36 | 


which is considerably larger than the values of Smin for 
fundamental spiral mode as given by Equation [33] under 
any circumstances. This shows that the fundamental high- 
frequency radial mode of oscillation with no circumferential 
component and small axial component is much more stable 
than the fundamental spiral mode where circumferential 
component is predominating with negligible axial component. 
(The case of 8 = 0 and m = 0 is essentially an axial mode 
which cannot be analyzed with the approximate scheme used 
in the present investigation.) The relative stability of high- 
frequency oscillations of different modes can be compared 
with Equations [28] and [31]. If the axial component of an 
oscillation is not the fundamental mode, but a higher mode, 
then the value of \, corresponding to this higher axial mode 
is larger than the value of Ag, corresponding to the funda- 
mental axial mode. The reduced frequency errs for 
the flow in the de Laval nozzle, 2 = Xé,-e / ) 

will also increase proportionately. From Fig. 2, it is seen 
that while Q increases, both A, and A; increase so long as 
Q is not too big. Equation [31] then shows that J(&)/—2noi, 
is bigger in the case of higher axial mode. The value of Smin 
for oscillations with nonfundamental axial component is 
hence larger than the value of S min for oscillations with funda- 
mental axial component. In other words, oscillations with 
nonfundamental axial component are more stable than the 
oscillation with fundamental axial component. 

For oscillations with fundamental axial component and 
successive higher modes of circumferential component, i.e., 
when £ takes successive higher integral values bigger than 
unity, the lowest zero 7 of Jg’(n) increases with 8. Thus, 
the ratio Q/A,, the reduced frequency Q, the product A,A,, 
and ultimately Sin, increase with increasing 8. Oscillations 
with nonfundamental circumferential component are, there- 
fore, more stable than the oscillation with fundamental cir- 
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cumferential component. Likewise, oscillations with non- 
fundamental radial components, i.e., when the zeros of Jg’(n) 
larger than 7 are taken, can be shown to be more stable than 
the oscillation with fundamental radial component. It is, 
therefore, concluded that the high-frequency oscillation com- 
posed of the three fundamental components in the circum- 
ferential, the radial, and the axial directions is the most un- 
stable mode. Let us call this mode the fundamental spiral 
mode. The investigation of the relative stability of this 
fundamental spiral mode alone in different rockets will indi- 
cate the relative liability of the occurrence of large amplitude 
oscillations in these rockets. 


Rockets with Different Radial Configurations 


If a nonburning rod is installed axially inside a solid rocket 
with tubular grain, the expression for Smin, as given in Equa- 
tion [25], is more conveniently written with r; = a and 
ro = 

1 
Smin = — [1 + ........... [37] 
where o = —(#,/2 é,)x which is a parameter of radial dimen- 
sions only. « is given by Equation [20], and 2, is negative 
on an interior burning surface and is positive on an exterior 
burning surface. 

¢ = mA;A,/(1 — a’) is the parameter critically dependent 
on the axial cliniers condition at x = A. 

no &w is given by the solution of the equation 


For the case of fundamental spiral mode where 8 = 1, the 
smallest positive values of for different values of a are 
determined graphically, based on the tabulated values of 
Bessel’s functions in (7). They are presented in Fig. 6 with 
no and Y,’()/Y1'(moa) plotted against a. It should be noted 
that when a approaches zero the function Y1’(mo)/¥1'(ma) 
behaves like a?. The parameter o is also calculated and 
plotted as shown in Fig. 7. When the radius of the non- 
burning rod increases, o increases monotonically and _ is, 
therefore, stabilizing. The actual value of Sin further de- 
pends on the parameter ¢, which depends on a number of 
factors involved in the boundary condition at the axial exit; 
for example, the nozzle geometry, the extent of the stagnant 
gas region, the shape of the end seal, and so on. Suppose 
the approximate boundary condition is taken and A; and A, 
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are determined from Fig. 2 in the manner explained in pre- 
vious section, it is found that when a increases from zero, 
the value of ¢(a)/¢(0) shows a slight decrease and then in- 
creases rapidly, as shown in Fig. 8 for the particular case of 
\d,| = 0.01 and 1/e = 10 with three different values of ). 
The value of o(a)/o(0) is also plotted in Fig. 7 which shows 
a rapid increase when a begins to increase from zero, but 
becomes nearly constant when a is bigger than 1/2. The 
product of the two parameters {-o for given value of A and, 
therefore, the value of Sin will increase steadily with in- 
creasing radius of the nonburning rod. In Fig. 9 the value 
of Smin is plotted against \ for a = 0, a = 1/4, and a = 1/2. 
Therefore, by inserting a nonburning rod into the tubular 
grain of a solid rocket, the system is stabilized. The stabil- 
izing effect expressed as the increase in Smin is larger if the 
radius of the nonburning rod as a fraction of the burning 
surface radius is larger. From Fig. 9 it is seen that the in- 
crease in Smin When X is less than 15 is rather small, 
and the destabilizing effect does not seem to be signifi- 
cant. However, if we look for the critical value of \ 
corresponding to a given value of S of the propellant, 
the stabilizing eftect is considerable in reducing the m:xi- 
mum value of \ compatible with unstable oscillations. In 
practical cases, the critical value of \ is only reached near the 
final configuration of the tubular grain. The reduction of \ 
by 1 or 2 can easily suppress the instability completely. 
The stabilizing eftect of the insertion of a nonburning rod 
into a tubular grain has long been confirmed by experiments 
and experience. 

Consider now, solid rockets using rod grain with exterior 
burning surface of radius a which decreases continuously in 
the course of operation. The radius of the combustion 
chamber is taken as reference length so that rz = 1 and 
\ is constant throughout the operation and 27, is 
Thus, one obtains from Equation [20] 


r; = a. 
positive. 


2 2 2 
1 —a? no2a? no? JL Ye'(n0) 


where is determined by = 

Ys'(mo). The values of y for fundamental spiral mode with 
8 = 1 are given in Fig. 6. The values of o calculated from 
Equation [38], are plotted in Fig. 7. It is seen that o is 
always bigger than unity and increases monotonically without 
limit as a approaches zero. During the period of operation 
of a rocket with rod grain, \ is constant but a decreases from 
the initia] fractional value toward zero. The decrease of the 
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parameter ¢ as given in Fig. 8 may decrease the value of 
Smin When a is reasonably large during the initial stage of 
operation, but is very soon overbalanced by the rapid increase 
of the parameter o toward the later stage of operation. 
Some calculated results are shown in Fig. 10 with S.nin plotted 
against a for several values of \. Rockets with rod grain will 
hence become much more stable toward the later stages of 
operation, and if a rocket with rod grain is stable during the 
initial and the early stage of operation, it will remain stable 
throughout the operation. In Fig. 10, it is seen that the value 
of Smin during the early stage of operation, for the value of 
\ as low as 12.5, remains bigger than unity, which is twice 
the value of Snin for rockets with tubular grain of equal \ 
and stabilized by nonburning axial rod. Therefore, it seems 
that a rocket with rod grain is unlikely, even if not impossible, 
to become unstable. So far as it is known, abnormal pressure 
peak has not been demonstrated in rockets with rod grain. 

Now consider solid rockets using rod-in-tube grain with 
both exterior and interior burning surfaces at r; = 1 and 
ry <1. The boundary condition given by Equation [14] is 
dropped and the boundary condition as given by Equation 
[15] is applied at both 7, = 1 andr, =a<1. Thus 


Jp'(n) + c¥ p(n) _, + c¥p'(na) 
Jan) + c¥ p(n) Jena) + c¥ gna) 


where %, is the normal velocity of the burned gas on the 
interior burning surface r, = 1 and is, therefore, negative. 
Let ¢o = = and ¢ = + 
(de/dn)y=m dn. Using the scheme of finding approximate 
solutions in previous sections, we can solve Equation [39] 
simultaneously with Equations [13] and [16]. The expres- 
sion of Sin is 


Smm = fl + 4» 
where ¢ = 9A;A,/(1 — a) 


nora?) Yg'(m)? + a¥ (ma)? 


. [40] 
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FIG. 10 MINIMUM PRESSURE INDEX OF INTERACTION Spin FOR 
ROD GRAIN WITH DIFFERENT LENGTH TO RADIUS RATIO \ OF THE 
NONBURNING SURFACE 


Both \ and a decrease with time during the operating period 
of a given rocket with rod-in-tube grain. o is a function of a 
only and is plotted in Fig. 7 when 8 = 1. We see that the 
value of o of solid rocket using rod-in-tube grain is in general 
slightly larger than the value of ¢ of solid rocket using tubular 
grain with nonburning rod having the same value of a. This 
slight difterence in ¢ does not seem to justify a conclusion 
about the relative stability of rocket using rod-in-tube grain 
as compared to rocket using tubular grain with nonburning 
rod, because the variation of the parameter ¢, being a func- 
tion of the rather uncertain axial boundary condition at 
x =X, may be a decisive factor in the comparison of these 
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two critical cases. But a rocket using rod-in-tube grain 


seems to be much more stable than the corresponding rocket 
using tubular grain having the same value of }. Therefore, 
if the initial radius of the central rod is sufficiently larger than 
the thickness of the cylindrical shell, so that the shell will be 
burned out before the rod, the system is less likely to become 
unstable. On the other hand, if the radius of the central rod 
is smaller so that the rod is exhausted before the shell, the 
system will become unstable just as the corresponding rocket 
using tubular grain would near the end of operation. 


Rocket with Restricted End Burning Grain 


It is also interesting to investigate the case of solid rockets 
with restricted end burning grain. This system looks very 
likely to be a one-dimensional system and can be analyzed 
by the method used in (4). We shall see whether the presence 
of the transversal components will be of importance in the 
stability analysis of such systems. 

The acoustical solution satisfying the periodicity in the 
circumferential direction and the radial boundary conditions 
is given as 


50) = nor)(e%* + [41] 


with a? = ¢ + ?, where ¢ and mp are the characteristic con- 
stants in x and r coordinates, respectively, 8 = integral char- 
acteristic constant in @ direction, and m = the lowest zero 
of Jg’(n). 


The boundary condition at axial exit x = d is 


.. [42] 
The boundary condition at burning surface x = 0 is - s 
ke 
(0) (0) 


By separating the real and the imaginary parts in Equa- 
tion [44], one obtains for neutral oscillations 


y(m — n)sinwr = —Y 
y(m — n) cos w7 = ym — (1 — 


where 


= 
2 


where the value of sin~! is taken in the proper quadrant con- 
sistent with the sign of cosw7. Equations [44] to [47] are 
closely analogous to Equations [4.1.1] to [4.1.4] in (4) with 
necessary changes of notations. 

If mo is taken to be zero, that is, if there is no radial com- 
ponent of oscillation,’ one has a = w = — or a/g = 1, and 
only longitudinal oscillations can exist. Equation [44] with 
a/t = 1 can also be derived from one-dimensional considera- 
tion. Since m — n = S(n/2) is of the order of unity or 
smaller, one sees from Equation [45], for this one-dimensional 
case, that Y must be of the order of unity, and that sin \¢ 
must be of the order of #, (cos 4 must not be small because 
the denominator is of the order of cos? Ag). Hence, dé will 
be close to integral multiples of + with a deviation of the 
order of %,. While \¢ varies in the neighborhood of 7 or 
integral multiples of x, S will have its minimum value, Sm, ., 
when Y is approximately zero. Thus, Smin is approximately 

iven by 


Sain = & 5 (1 + A,) {48] 


27 
If the deLaval nozzle of the rocket is “very short,’ and the 

“short nozzle boundary condition’’ applies as explained in (4), 


then A, = (y — 1)/2 and 


If the nozzle is not “very short,’ Equation [48] gives the 
value of Smin for pure axial neutral oscillation of frequency 
\é with A, determined from Fig. 2. Since A, increases with 
Ag, we see that Smin for higher axial modes are larger, that is, 
the higher axial modes are more stable than the lower modes. 
Qualitative stability behavior of pure axial mode in rockets 
with end burning grain is, therefore, closely analogous to 
that in liquid rocket of similar geometry with concentrated 
combustion at the injection end. However, there are two 
fundamental differences between the present results and 
those in (4). Since the length of the combustion chamber of 
solid rocket increases with time during operation, the effec- 
tive length of the nozzle e is decreasing. Therefore, the 
values of A, and Sin for a given axial mode, decreases during 
operation. In other words, a solid rocket with end burning 
grain destabilizes itself. The unstable ranges of the time lag 
for given value of S of the propellant also varies, due to the 
change of w = ~ as can be seen from Equation [47]. There- 
fore, if the value of Smin of the rocket with end burning grain 
becomes smaller than the value of S of the propellant, the 


(1 + cos 2 + (1 — cos 2 dé) (“*) (A,2 + A?) — 2A, sina | 


7 2 
2A, cos 2 + E (“*) (A? + 4] 
Wie 


Y=- 


The value of the average pressure index S of interaction corre- 
sponding to neutral oscillation with axial component ¢ and 
radial component ) = (w? — ¢)'/? can be found from Equa- 
tion [45] as 


‘eon n 1 Y2 
S= = 1— X 


and the corresponding critical time lag 7 as 


(1 + cos 2 + (1 eos 26) + A,?) — 2A; (“*) sin 2 XE 


s 


axial mode will very likely become unstable. Another funda- 
mental difference is that the analysis in (4) deals with in- 
trinsic instability only where the rate of supply of the com- 
bustible mixtures is constant, while the present analysis deals 


‘It can be shown from the solution of the wave equation, that 
when » = 0, 8 must be zero in order that both of the two radial 
bound: conditions at r=1 and r =0 may be satisfied. 
Physically, this means the absence of any transversal com- 
ponents of oscillation or the flow is strictly one-dimensional. 
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with systems with varying rates of supply of the combustible 
gases. 

Suppose now a radial mode np associated with circumfer- 
ential mode of order 6 = 0, 1, 2, etc., is superposed = 

axial mode of frequency \é, then w/t = I + (*)'] is 
bigger than unity. If 9/¢ is of the order of unity, for example, 
in the case of fundamental spiral mode with 6 = 1, 7 = 1.84, 
ht < mw, and A = O(1), then w/¢ is still of the order of unity. 
Following the reasoning for the case of pure axial mode, one 
sees that A€ is still close to integral multiples of x with devia- 
tion of the order of %. The value of Smin for the neutral 
spiral mode as specified is still given by Equation [48]. 
However, A, is not known even approximately. Because 
when o/é is of the order of unity, the time-wise variations 
of the flow properties due to the transversal and the longi- 
tudinal components of oscillations are also comparable. 
Thus, the approximation of taking over the one-dimensional 
results in (5) is invalid. Until an analysis of nozzle flow 
with transversal oscillations is made, very little can be said 
about A, as a function of dé, 7, and probably 8 as well. The 
investigation of the qualitative stability behavior of this type 
of oscillation should, therefore, be studied more carefully. 
If, on the other hand, m/¢is large, either in the case of funda- 
mental spiral mode with large \ or in the case of higher 
members of the radial modes superposed on fundamental 
axial mode, then w/é is much larger than unity of the order 
of 1/a%,. With m — n of the order of unity, it is no longer 
necessary that Az is in the neighborhood of some integral 
multiple of x. The frequency of the axial component of the 
neutral oscillations and, therefore, of the unstable oscillation 
may be quite different from and cannot possibly be identified 
with any of the natural organ-pipe frequencies of the com- 
bustion system. 

The simplifying scheme of calculating Simin by Equation 
[48] as 1/2y(1 — X) is no longer valid. Numerical deter- 
mination of s and 7 for different values of mo, £, #, and n can 
however, be carried out with Equations [46] and [47] in a 
manner analogous to that in (4). Fig. 2 can still be used in 
determining A, and A; as a rough approximation. The 
qualitative behavior of the variation of Sinai. will be expected 
to be not much different from that of pure axial oscillation 
in solid rocket with restricted end burning grain. Al 
Summary of Conclusions 


1 High-frequency unstable oscillations in solid propellant 
rockets are excited by the sensitivity of the rates of the com- 
bustion processes to oscillations of the burned gas in the 
combustion chamber. The following simplified model is 
assumed. The solid propellant first decomposes into inter- 
mediate gaseous products without appreciable energy release. 
These intermediate gaseous products are then activated 
during a time lag 7, near the end of which these intermediate 
gases undergo a series of combustion reactions at a very large 
rate which releases practically all the chemical energy made 
available through combustion. The interaction between 
the rates of the combustion reactions and the gas oscillations 
is expressed in terms of the pressure dependence of time lag r 
through the following integral relation 


t 
p™t’)dt’ = const. 
t-—r 


where m represents the average extent of interaction including 
temperature effects. The rate of primary decomposition of 
the solid propellant on the burning surface is assumed to be 
proportional to the nth power of the instantaneous gas pres- 
sure exerted on the solid surface. Both of the two constants 
m and n can only be determined by experiments at the present 
time. The effective over-all interaction parameter S is found 
to be the difference between the two indices S = m — (n/2). 
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Small high-frequency oscillations of a given mode of either 
spiral or nonspiral type can become unstable when both the 
value of S of the propellant is larger than a certain minimum 
value, and the value of 7 lies in one of the unstable ranges of 
7 for that mode. The minimum value of s of a given mode 
is in general of the order of or smaller than unity, and depends 
on the rocket geometry and the burning rate of the propellant. 
The unstable ranges of 7 of a given mode depend on the 
rocket geometry and the value of S of the propellant. 

2 For solid rockets with cylindrical burning surface or 
surfaces, the frequency 6 of the circumferential component 
of the unstable oscillations of the spiral mode is always equal 
to the frequency of the acoustical circumferential mode of the 
combustion system, i.e., 8 = integral values. The frequency 
n of the radial component is close to the frequency 7 of the 
acoustical radial mode with a fractional deviation of the 
order of |d,|, the dimensionless normal velocity component of 
the burned gas on thé burning surface. The frequency ¢ of the 
axial component is, however, significantly different from 
the natural organ-pipe frequency kr/d. The frequency a of the 
variation of flow properties at a given point is close to 
the frequency mo of the acoustical radial mode. When the 
unstable spiral mode of oscillation amplifies, the circumfer- 
ential component amplifies at a much faster rate, and meas- 
urement of the unstable oscillations at finite magnitude may 
not be able to reveal the associated radial and axial com- 
ponents of oscillation. 

3 Despite the fact that the burning surface is in the 
radial direction in solid rockets with cylindrical burning 
surface, which seems to be most susceptible to the oscillations 
in the radial direction, the nonspiral oscillation with pre- 
dominating radial component but no circumferential com- 
ponent is more stable than the fundamental spiral node. 
Among all the possible unstable high-frequency oscillations 
either of the spiral or of the nonspiral mode, the fundamental 
spiral mode is the most unstable mode. 

4 Solid rockets using tubular grain or perforated pro- 
pellant grain with interior burning surface are found to be the 
most unstable grain shape as compared to other shapes 
studied in the present paper. It is shown that: (i) A solid 
rocket using tubular grain become more unstable, that is, it 
destabilized when the ratio \ decreases with time of opera- 
tion. The minimum value of S of such systems decreases 
very rapidly when ) decreases until certain critical range of 
d is reached, after which S,,;, decreases at a much slower rate. 
The rocket that is stable for most part of its operation will 
suddenly become unstable when the value of Smia of the 
system becomes sufficiently below the value of S of the pro- 
pellant. (ii) The rocket with tubular grain is stabilized by 
introducing a nonburning rod into the axial cavity. Both 
the value of Simin is increased, and the critical range of 2 is 
reduced. If the value of \ of the tubular grain near the end 
of its operation is close to the critical value of \ of the rocket 
system, the introduction of a nonburning rod can suppress 
completely any instability. The stabilizing effect increases 
when the radius of the rod increases. (iii) Rockets with non- 
erosible end seal are likely to be more stable than those 
with erosible end seal of the same initial inside diameter. 
The effect of increasing the inside diameter of the end seal 
is destabilizing. 

5 Rockets using rod grain with exterior burning surface 
is the most stable configuration among those studied here. 
This system becomes more stable when the radius of the 
burning rod decreases in the course of its operation, especially 
near the end of the operation. 

6 Rockets using rod-in-tube grain with both interior and 
exterior burning surfaces are more stable than the corre- 
sponding rockets with tubular grain. If the initial radius of 
the central rod is sufficiently larger than the thickness of the 
shell, this system will remain more stable than those with 
tubular grain. 
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= 7 By decreasing the parameter e which represents the 
>» effective length of the nozzle, the rocket system is destabilized. 
8 The stability behavior of high-frequency small dis- 


: lyzed strictly by one-dimensional model only when there is 

eS | no transversal component of oscillation present. The re- 
oes stricted end burning rocket destabilizes itself toward the 
pure longitudinal oscillation during the course of its opera- 
tion. If the ratio of the length of the combustion chamber 
space to its radius is large, the frequency of the axial com- 
ponent of the unstable fundamental spiral mode will not | 
necessarily be close to the fundamental natural organ-pipe — 
frequency of the combustion system. 


Postulate series solutions in terms of the mean exit velocity 

ii, as follows 


p= pO + + +... [Al] 


Substitute these into Equations [7] and equate terms of the 
same power of #,. The ratios of the steady-state velocity 
components to #, are less than unity and are considered as 
terms of the order of unity or smaller. The equation for the 
Oth order solutions are 


1 
- xp =0 


i 


= —§, wile 


Equation [A2] can be reduced to the single partial-difterential 
equation 
(a? + = +4 5,, ++ +5 = 0..[A3] 
where V? is the Laplacian operator in cylindrical coordinates. 
The solution of Equation [A3]is easily obtained by the method 
of separation of variables and is well known as the cylindrical 
acoustical wave solution. This solution satisfying particular 
sets of boundary conditions are given as Equation [8] and 
Equation [41], respectively. 

The differential equation for the first correction 6 is ob- 
tained as 


} 
(a? + = ia + 
Ue z 


The inhomogeneous part of Equation [A4] is a known func- 
tion of x, r, and @ with a magnitude of the order of |é| or = 


which is considered as small compared to unity. Call this 
known function F(z, re’ and expand it into Fourier series 


whose jth term is F;(r) exp [=| The solution 
can be obtained from Equation [A4] in the form of infinite 
series 


j=l 


turbances in restricted end burning solid rocket can be ana- na 
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Jp(r0’) Y )F; (r) i| 
1/2 


i g and Yg are #th-order Bessel functions of the first and 


the second kind, respectively. 

Consequently, the first-order corrections 7,6“ can be found 
at least numerically after the determination of the steady- 
state solution and the acoustical solution of the particular 


system. 
APPENDIX 2 


It is experimentally found that the rate of regression of 
the solid propellant surface during the rough combustion 
period is abnormally high even if the increase of burning rate 
due to the increase of mean chamber pressure is taken into 
account. This phenomenon is hardly compatible with any 
linearized small perturbation theory. For example, if the 
burning rate is sensitive to the pressure variation, then a 
small sinusoidal pressure disturbance will decrease the burn- 
ing rate during pressure defect period to such an extent as 
will cancel to a large part the increase of the burning rate 
during pressure excess period. The mean burning rate and 
the mean combustion chamber pressure should not increase. 
If the oscillation is of sufficiently large amplitude, the non- 
linear eftect may produce significantly larger increase of 
burning rate in pressure excess period than the decrease in 
pressure defect period, and result in an increase in mean 
burning rate and mean chamber pressure. An alternative 
explanation of the increase of mean burning rate is the de- 
pendence of the burning rate on drifting velocity along the 
burning surface. Under the same pressure, the burning rate 
is smallest when there is no drifting velocity along the surface. 
If an oscillation of the drifting velocity is introduced into the 
gas system, the burning rate will be increased both in the 
velocity excess and in the velocity defect period. It seems, 
therefore, that a small perturbation analysis based upon this 
velocity mechanism can result in increasing mean burning 
rate and mean chamber pressure. However, the drifting 
velocity in the axial direction along the burning surface is 
not zero in steady-state operation. It is only the circum- 
ferential velocity component that vanishes in steady state- 
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and only a small velocity oscillation in circumferential direc, 
tion that can possibly increase the mean burning rate. Since 
the burning rate depends on the magnitude of the drifting 
velocity, i.e., (u? + w?)'”, and since # is not zero while w = 0, 
the instantaneous drifting velocity is [#2 + 2au’ + + w’ 2)'/: : 
where prime indicates small perturbations. Both of the two 
nonlinear terms u’? and w’?, which are independent of the 
sign of the perturbations, contribute to the increase of the 
mean burning rate and, consequently, the mean chamber 
pressure, just like the nonlinear terms of the pressure per- 
turbation do. These terms of nonlinear character should not 
be considered in a linearized stability analysis of small dis- 
turbances because they are small quantities of higher order. 
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It is, therefore, inferred that the increase of the mean burning 
rate and the mean chamber pressure is a “nonlinear’’ property 
and should not be considered in the linearized analysis for the 
stability of small perturbations superposed on the steady flow. 

A stability analysis based on the exciting mechanism as 
represented by Equation [4] can, of course, be made without 
much complication in analytical procedure. But due to the 
introduction of a second empirical interaction parameter b, 
the ultimate stability limit would involve two parameters, 
§ and b, which have to be presented in a plane, and the deter- 
mination of such stability limit would involve much more 
numerical work than is intended in the present analysis. 
The following analysis indicates that the velocity mechanism 
alone is not very likely to induce instability and, therefore, 
the stability analysis is made on the pressure mechanism only 
as given in Equation [5]. 

When the velocity mechanism only is considered, the in- 
stantaneous burning rate m,(t) is 


mi(t) = = po(l +b7)... [A6] 


where tho = pi is the burning rate at zero drifting velocity, 
and = + 6V) under the approximation # < 1. 
The boundary condition at the burning surface r. = 1 in the 


ease of a rocket with tubular grain is 


v 
= 1 1... 
5 5 | 


Equation [14] is replaced by c = 0, and Equation [15] by 


It is interesting to note that 7 is now a function of x. This 
indicates that the method of separation of variables is not 
applicable to the solution of the problem. To avoid the 
mathematical complications arising from this situation, we 
shall assume that so far as the over-all effect is concerned, 
the propellant elements distributed over the entire burning 
surface respond to the same axial velocity variation at the 
mean station ¥. In actual situation each propellant element 
responds to the local velocity variation. These local velocity 
variations are of different amplitudes and are in general out 
of phase. Therefore certain cancellations of the exciting 
contributions of different propellant elements exist. The 
total amount of excitation in the actual system will be less 
than that in the assumed system if £ is selected as the posi- 
tion where the amplitude of the axial velocity variation is 
the largest. In other words, the actual system will always be 
more stable than the most unstable one of the idealized sys- 
tems with arbitrary ¥ from 0 to A. Thus we write 


J3'(n) 
= —tady| 1 4+ — bi(r)- 


au, tan &r 


Equation [A8] is solved simultaneously with Equations 
[13] and [16], with the approximate procedure as explained 
previously. The critical values of 6%, corresponding to neutral 
oscillations of frequency w is given as 


G + [G2 + (|AH|? — + F2)]'2 
where 
9 ee 2 2 
2 wig no” 
w ra) 2 win no? 
A=A,+i%A; ‘ 
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The complex quantity — = & + 7é, is still determined by the 
axial boundary condition, Equations [29] and [30]; and the 
magnitude of A for the fundamental spiral mode of oscil- 
lation is considerably less than unity. 

The magnitude of H depends on the position ¥ where the 
generation of the burned gas is assumed to be concentrated 
in the evaluation of the variation of the rate of burned gas 
generation. If ¥ = \ we have H = 1, and if ¥ = Owehave 
H = 0. For arbitrary values of 


sin? + sinh? + sinh? 
cos? sinh? cos? + sinh? 


Since &;/£, is of the order of 2 or 3 or even larger, |H| will 
remain very close to unity and become very small when £ 
approaches zero. The magnitude of AH is thus given by 
A| which is considerably less than unity. 

Now Equation [A9] shows that if | AH | <| the right-hand 
side is negative no matter whether the upper or the lower sign 
before the bracket is taken. Since bi, of the solid propellants 
are positive, we conclude that the fundamental spiral mode 
would be always stable if the velocity mechanism were the 
only exciting agent. Physically this means that the velocity 
mechanism is not able to supply sufficient amount of energy 
to the disturbance to compensate the energy dissipation due 
to the outflow through the nozzle associated with the velocity 
oscillation. The disturbance will lose its energy and die out 
eventually. 

Even if we consider the extreme case where | A//| > |, 
the minimum value of bu, compatible with any unstable oscil- 


lations 1S given by 
l + A A, (1 
G ; r no? 


(bite)min 


A-10 
AH| |AH| -1 


Sample calculations show that (bue)min remains considerably 
bigger than unity while the magnitude of ba, for solid pro- 
pellants is in general considerably less than unity, even for 
those propellants that exhibit instability in tubular grain. It 
is therefore inferred that the velocity mechanism is not the 
essential mechanism in exciting instability, at least for the 
cases that are considered in the present investigation. 

Finally it should be emphasized that this velocity mecha- 
nism can be neglected without seriously affecting the 
qualitative results only in the analysis of the stability of 
small disturbances, that is, the analysis of the behavior of the 
initial phase of the development of the large oscillations in 
solid rockets. Once the unstable small disturbances have 
grown up to considerable magnitudes, the velocity mechanism 
will become important and even may be the predominating 
mechanism, as is in the case of unstable spiral modes at large 
amplitudes. 


Are You Missing Something? 

Do you know that an annual subscription to the 
Journal of the British Interplanetary Society can be 
obtained for $1.00? By a special exchange arrange- 
ment with the BIS, these Journals can be ordered by 
ARS members from the New York office. 

Some 1953 subscriptions are still available, as 
well as a few back volumes. A list of immediately 
available issues will be sent upon request. 

Address: Secretary 
American Rocket Society 
29 West 39 Street 
New York 18, N. Y. 
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Excitation of Oscillations by Chemical 
Reaction’ 


JACK LORELL?’ and HENRY WISE’ 


Jet Propulsion Laboratory, California Institute of Tech- 
nology, Pasadena, Calif. 


EVERAL theoretical treatments (1-5)‘ have recently 
dealt with the problem of combustion instability in the 
rocket motor. These analyses have shown that inherently 
unstable (oscillatory) systems may result from the interaction 
of the hydraulic, mechanical, and hydrodynamic components 
of the rocket motor. The purpose of this paper is to suggest 
that, in addition to the factors named, the chemical kinetics 
of the combustion process may also be a source of excitation 
for pressure oscillations. 

In an open reaction system, such as a rocket motor, in 
which there exist a continuous supply and removal of reac- 
tants, the concentration changes which occur are governed by 
the rate processes of mixing, phase change, and chemical 
reaction. For the present discussion, the system oscillations 
caused by variations in physical parameters are separated 
from the fluctuations which may originate within the chemical 
system as a result of the interplay of the reaction kinetics. 
Such time-dependent kinetic variations in the concentration 
of a particular molecular or atomic species can couple with the 
hydrodynamics of a system such as a rocket motor. In gen- 
eral, chemical] systems do not fluctuate about their equilibrium 
because chemical energy is normally transformed irreversibly 
into heat, not into some form of potential energy as experi- 
enced in a mechanical vibration. 

In considering only the chemical kinetics of a system made 
up of a series of consecutive first-order reactions, the possi- 
bility of true time-dependent oscillations in the concentration 
of a given reactant is ruled out (6). However, it has been 
shown that the differential rate equations applicable to an 
idealized open reaction system of two autocatalytic reactions 
can give rise to undamped oscillations in the concentration of 
intermediate species (7). Such a system of equations has been 
employed by Walsh (8) and Frank-Kamenetsky (9) in the 
treatment of the periodicity of cool flames operative in the 
mechanism of two-stage ignition of hydrocarbons (10). Other 
periodic chemical reactions have recently been reviewed (11). 
The analysis has been extended to a chain of autocatalytic 
reactions in which the precursor is introduced by a first-order 
process and the other substances are removed by first-order 
processes (12). Such a system exhibits almost undamped 
oscillations under certain conditions. 

Since the study of the transient behavior of nonlinear sys- 
tems (such as are encountered in many chemical reactions) 
involves a complex mathematical treatment, the analysis sug- 
gested herein deals with isothermal, closed (with respect to 
the concentration of reactants) systems involving kinetics of 
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first and second order. It is obvious that, in a reaction pro- 
ceeding along the sequence, reactants — intermediates — prod- 
ucts, the rate of this over-all reaction is governed by the slow- 
est step in the mechanism. In such a system characterized 
by the absence of a feedback loop, the concentration of each 
species does not display oscillations but rises or falls mono- 
tonically to the equilibrium value. On the other hand, in the 
presence of a chain or an autocatalytic mechanism as ex- 
hibited by some chemical reactions of interest in combustion 
(13), periodic variations in the concentration of some of the 
species cannot be ruled out. 

The oscillatory behavior of such a system may be analyzed 
by linearizing the reaction-rate equations about the equilib- 
rium point, resulting in a system of linear differential equations 
with constant coefficients. If, for such a system of linear dif- 
ferential equations in the displacements from equilibrium of 
the various chemical species, the characteristic complex fre- 
quencies are all real and negative, the system is nonoscillatory. 
However, a resulting complex frequency with nonzero imagi- 
nary part indicates oscillatory behavior. An oscillating system 
is of particular interest since it contains a natural frequency 
which, upon excitation by an external disturbance, may result 
in an observable resonance. 

As an example of an oscillation induced by the kinetics of a 
reacting system, a reaction involving two consecutive auto- 
catalytic reactions may be considered. Such a system is of 
particular interest since it has been applied in the analysis of 
ignition of hydrocarbons (10). This system in which the for- 
mation of the intermediates A and B (whose respective con- 
centrations are given by r; and rz) exhibits autocatalysis may 
be represented by the rate equations 


These equations lead to a periodic time dependence of 1; and 
r. about their equilibrium values 


7, = ks/ke 
= ki/ke 


For small displacements 2, and 22 from their equilibrium 


values, one obtains from Equations [1] and [2] 


where 

= — (ks/k2) and x2 = re — (ki /ke) 
This set of simultaneous, first-order, differential equations 
may be solved directly, resulting in an expression for 2; and 2; 
which has the characteristics of a harmonic oscillator. The 
nature of the oscillation is also apparent from an examination 


of the characteristic equation of the differential expressions 
[3] and [4] 


= 
which is a quadratic equation having for its discriminant 
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Since the discriminant is negative, the solution of Equation 
[5] contains purely imaginary roots. Consequently this sys- 
tem exhibits undamped oscillations with a frequency » equal 
to 


[7] 


Further analysis shows that x; and 22 are out of phase by 
90 deg. It should also be noted that the rate constants are 
exponential functions of temperature. Therefore the fre- 
quency of the oscillation rises exponentially with temperature, 
according to Equation [7]. The order of magnitude of this 
frequency at 1000 K is from 10 to 10° sec for a reaction with 
an activation energy of 50 to 60 keal and a pre-exponential 
factor of 101%, 
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Production of Chemicals by Gas- 
Dynamical Methods’ 


GORDON J. MACLEOD? 


Daniel and Florence Guggenheim Jet Propulsion Center 
California Institute of Technology 
Pasadena, Calif. 


EVERAL proposals for the production of chemicals by the 
use of gas-dynamical methods have been advanced in recent 
years. Briefly, these methods involve (a) the utilization of 
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the possibility of ‘freezing’ equilibrium gas compositions by 
very rapid chilling (as by expansion through a de Laval 
nozzle), or (b) the production of chemicals under nonequilib- 
rium conditions, for example, by rapid heating and cooling 
through proper exposure to compression and rarefaction 
shocks. 

It is the purpose of this letter to call attention to a study 
entitled “Some Considerations in the Application of a Gas 
Turbine Cycle to the Manufacture of Nitric Oxide,” which 
was carried out by the author in 1951 and 1952 at the Cali- 
fornia Institute of Technology, in partial fulfillment of re- 
quirements for the degree of mechanical engineer. The ab- 
stract of the thesis (146 pages) is given below in order to in- 
dicate the scope of the work: 

“Problems associated with the application of a gas turbine 
cycle to the manufacture of nitric oxide have been investi- 
gated. The feasibility of quenching the nitric oxide decom- 
position reaction with a de Laval nozzle is demonstrated. 
Thermochemical studies show that yields of nitric oxide ap- 
proaching one per cent are attainable in a gas turbine cycle. 
It may be possible to overcome the serious problem of blade 
cooling by using either transpiration cooling or film cooling. 
Transpiration cooling appears to be the more promising 
method. The coolant requirements amount to approximately 
one per cent of the main stream mass flow per cooled turbine- 
blade row. The effect of transpiration and film cooling upon 
cycle performance is negligible. Although all of the problems 
associated with the application of a gas turbine cycle to the 
manufacture of nitric oxide appear to be surmountable, it has 
not been demonstrated conclusively that the turbine cycle 
will be competitive with conventional methods for nitrogen 


fixation.” 


Ethylene Oxide as a Monopropellant’ 
WILLIAM CLAY ROBISON? 


Daniel and Florence Guggenheim Jet Propulsion Center 
California Institute of Technology 
Pasadena, Calif. 


HERMODYNAMIC ealculations to determine the 

theoretical performance of ethylene oxide as a mono- 
propellant have been carried out for various possible de- 
composition reactions. The performance calculations were 
carried out, by using standard evaluation procedures,’ for 
two possible exothermic decomposition reactions, one of 
which leads to the formation of CO and CH4,, whereas the 
other leads to the production of CO, C,and H». For the proc- 
ess leading to carbon formation, two limiting cases were con- 
sidered: (a) no slippage between the carbon particles and 
the gases during expansion with thermodynamic equilibrium 
being maintained at all times; and (b) complete deposition of 
carbon in the combustion chamber. The results of calcula- 
tions, for a chamber to exit pressure ratio of 20.42:1 and in- 
jection of the liquid monopropellant at its normal boiling 
point of 10.7 C, are summarized in Table 1, where T, and T, 
denote the adiabatic flame temperature and the nozzle ex- 
haust temperature, respectively, c* is the characteristic gas 
velocity, C; represents the nozzle thrust coefficient, and J,, 
denotes the specific impulse. 
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TABLE 1 SUMMARY OF PERFORMANCE CALCULATIONS FOR 
C.H,O MONOPROPELLANT 


Decomposition c* Cy le 

reaction (°K) (°K) (m/sec) (—) (sec) 

CH, 1430 908 1230 1.295 162.2 
C:H0O0 = CO + 


2H, + C (carbon 
in equilibrium 


with gas) 638 304 891 1.400 127.6 
C.HO = CO + 

2H. + C (car- 

bon deposited in 

rocket chamber) 638 270 1025 1.385 105.4 


The factors leading to carbon formation in combustion are 
not clearly understood, although a considerable amount of 
work has been done on this subject in recent years. An ex- 
tensive literature survey on the kinetics of ethylene oxide 
decomposition has been carried out. The results of this 
study are summarized in the thesis from which this letter is 
abstracted. 


Experimental Measurement of the Heats 
of Dissociation of Hydrazine-Water and 
Hydrazine-Methyl Alcohol Systems' 


WILLIAM CLAY ROBISON? 


Daniel and Florence Guggenheim Jet Propulsion Center 
California Institute of Technology 
Pasadena, Calif. 


LTMAN and Adelman? have studied the equilibrium be- 

tween hydrazine and water in the vapor phase by following 
the pressure change in a constant-volume apparatus. They 
found the heat of dissociation of hydrazine hydrate to be 
13.97 Keal/mole. Their data on hydrazine hydrate have 
been checked within the limits of our experimental error 
Furthermore, the heat of dissociation of NeH,-CH;OH has 
been found to be 8.6 + 0.3 Keal/mole, using essentially the 
same experimental technique as Altman and Adelman. 
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(Continued from page 94) 


Collateral Reading 

4. “The Dynamics and Thermodynamics of Compressible 
Fluid Flow,” vol. 1, by A. H. Shapiro, The Ronald Press, Nev 
York, 1953, pp. 147-152 (sections on supersonic inlets and di- 
fuser efficiency). 

5. “Aerodynamics of Propulsion,” by C. Kuchemann and J. 
Weber, McGraw-Hill Book Co., Inc., New York, 1953, chap. 7 
(section on the ramjet engine). 
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Jet Propulsion News 


ALFRED J. ZAEHRINGER, Thiokol Chemical Corporation, Associate Editor 


JOSEPH C. HOFFMAN, General Electric Corporation, Contributor 


Rockets and Guided Missiles 


DDITIONAL details of the Nike guided missile system 
have been made available. The Nike is a ground-to- 
air missile which was developed jointly by Army Ordnance, 
Douglas Aircraft, Bell Telephone Laboratories, Western 
Electric Company, and Aerojet-General Corporation. Nike 
is about 20 ft long and 12 in. in diam. A solid propellant 
booster is used for launching while an acid-aniline rocket 
engine provides sustaining thrust for the missile. Range is 
stated to be from 10 to 20 miles with 30,000 to 50,000 ft 
altitude. Guidance is by a beam rider system. It was an- 
nounced that the Nike system will be used for the air defense 
of Washington, D. C. Other U. 8. cities are scheduled to 
be protected by Nike batteries in the near future. (See Figs. 
1 to 5.) 


LOKI is a new member of America’s guided missile family 
to be announced. Loki is a solid-propellant barrage-type 
antiaircraft missile based on the Taifun rocket developed in 
Germany during World War II. 


THE Air Force has released preliminary flight information 
on the Rascal, a Mach 1.5, air-to-surface missile developed by 
Bell Aircraft. Completely rocket-propelled, the Rascal is 
designed to be launched from the bomb bay of a B-36 or a 
B-50 bomber. The launching gear is said to be similar to 
that used for the X-1A, Bell’s rocket-powered aircraft. To 


keep the mother ship well out of the range of antiaircraft 
missile defenses, the Rascal is launched about 100 miles 
from its surface target. 
from a preset programming 


After the drop, initial guidance comes 
device within the missile. 


Courtesy U. 8. Army 


FIG. 2 NIKE 


Courtesy U.S. Army 
FIG. 1 NIKE 
This remotely controlled U. 8S. Army Ordnance launcher poises a 
partially erected Nike missile against the New Mexico sky. When 
the launcher has reached an 85-deg angle, the missile will be in posi- 
tion to fire 


A booster charge has started the Army’s antiaircraft Nike on its 
supersonic quest for a flying target. The lower portidn is expended 
after an initial thrust, and the guided missile proper continues toward 
its target, guided by an electronic ‘‘brain” 


Epitor’s Note: The information reported in this Section has been selected from approved news releases originating with the 


Department of Defense, private manufacturers, universities, etc., and from published news accounts in journals and newspapers. The 
reports are considered generally reliable, although no attempt has been made to verify them in detail. 
Marcu-Aprii 1954 


/ 13 


R 
by 
) 
‘ats 
o be 
have 
: 


Courtesy U. S. Army 
FIG. 3 


Left: Nike nears the target plane. A smoke pot under wing 
enabled the photographer to follow the action 
Right: Nike flashes beneath the target’s wing 


Midcourse and terminal guidance is provided by radar from 
the mother plane. The Rascal uses its rocket engine to ac- 
celerate in level or climbing flight, reaching its maximum 
speed when the propellants are exhausted. Terminal portion 
of the flight is said to be either a long glide or a quick push- 
over into a steep glide. 


TWO other British rocket engines have been disclosed by 
the Royal Aircraft Establishment. One, a 900-lb thrust 
unit, with an over-all diameter of 18 in., was used to power 
the transonic research model aircraft produced by Vickers. 
The propellants, hydrogen peroxide and methanol, were fed 
into the uncooled combustion chamber by high pressure gas 
at 500 psi. Although single-shot performance (70 sec) of 
the rocket motor was considered good, the British found that 
it was not suitable for repeated use. The RAE has had 
marked success with the Beta 2 rocket engine. One of the 
Beta’s main features is its turbopump propellant feed system. 
Driven by steam generated from hydrogen peroxide, which 
is used as the oxidizer, the turbine is of simple construction 
and said to be extremely reliable. 


ENGINEERS of the Cook Research Laboratories, Chi- 
cago, Ill., are testing parachutes at near supersonic speeds 
on the rocket track at Edwards AFB, Calif. The sled (Figs. 
6, 7) is powered by a North American Aviation rocket motor 
and is capable of being accelerated from standstill to 1500 
mph in 4.5 sec. As the sled nears the end of the 10,000-ft 
long precision-aligned track, a brake scoops up water from 
a trough between the tracks and throws it out into the air- 
stream, thus utilizing the energy of the moving water to 
stop thesled. Cook engineers are also conducting near super- 
sonic tests on parachutes with its ‘Skokie Missile” (Fig. 
8) which is dropped from aircraft at high altitudes. The 
pods at the tail of the missile house cameras which record the 
action of the parachute during drop conditions. The use of 
a landing spike (Fig. 9) facilitates recovery after the test. 


Aircraft 


ON DECEMBER 12, 1953, Major Yeager of the USAF 
flew the Bell X-1A at a speed of 1600 mph. The flight took 
place over Edwards AFB, Calif. Powered by an RMI rocket 
engine, the X-1A differs from the X-1 in having a turbine- 
driven propellant pumping system. 


THE U. 8S. Navy has awarded North American Aviation 
a contract for the production of FJ-4 carrier-based jet fight- 


FIG. 4 


Left: Nike explodes at microsecond of intercept na 


Right: Firey fragments of exploded Nike missile rip target 


FIG. 5 


Left: Wing shattered and afire, target plane starts final plunge 
Right: Motor tears from wing as target plummets to earth 


ers. Called an improved version of the FJ-3, the power 
plant is to be a Curtiss-Wright J-65 jet engine which will 
place the speed of the craft “above 650 mph.” 


ANOTHER Navy craft, the F3H Demon, has been placed 
in production. Manufactured by McDonnell Aircraft, St. 
Louis, Mo., the Demon is a single-jet, all-weather carrier- 
based fighter model powered by a Westinghouse J-40 turbo- 
jet engine. Later models are scheduled to be powered by 
an Allison J-71 power plant. The Demon combines inter- 
ceptor speed and maneuverability with the payload of an 
attack bomber. Specifications: wing span, 35 ft 4 in.; 
length, 59 ft.; height, 14 ft. Top speed was not disclosed 
by the Navy. Armament is to consist of 20-mm cannon 
and air-to-air rockets. 


THE de Havilland Corporation stated that jet-powered 
Comets will be placed on the trans-Pacific and trans-Atlantic 
passenger runs. The Comet 2, seating 44 to 48 passengers, 
operating at nearly 500 mph at ranges of over 2000 miles, 
will be delivered to Canadian Pacific Airlines during 1954 for 
the Vancouver-Sydney run. The Comet 3, cruising at better 


than 500 mph, at ranges of 2500 miles, will carry 58 passengers 
in the trans-Atlantic run for British Overseas Airways Cor- 
poration. 
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FIG. 6 ROCKET-PROPELLED TEST SLED 


Jet Propulsion Engines 


A RAMJET engine that needs no boost to reach its operat- 
ing speed is under development at Rensselaer Polytechnic 
Institute. Mechanical inlet valves and explosive ignition 
are said to be the basis for the RPI design. The information 
available indicates that the operating cycle of the ramjet 
will resemble a pulsejet engine. 


RYAN Aeronautical Company has received production 
orders for components on the General Electric J-57 jet en- 
gine and for afterburners on the newest model of the Wright 
turbojet. 


JAPAN soon expects to complete its first postwar turbojet 
engine, the Omiya Fuji JO-1, produced by Omiya Fuji 
Industries. The JO-1 is essentially an improved Ne-20 
(a Japanese version of the German BMW 003) which pow- 
ered the ‘Kikka.” The engine is of the axial-flow type and 
features a single-stage turbine and 8-stage compressor, 
delivering two thrusts around 2200 Ib. 


FRENCH jet-engine development was given a substantial 
boost with the announcement that the Snecma Vulcan turbo- 
jet completed qualification test at 10,000-lb sea-level thrust. 
The single-compressor axial-flow turbojet is expected to 
qualify as an 11,000-lb thrust production unit in 1955. Other 
jet engines developed by the French Snecma firm include the 
Atar 101D and 101E. The 101D produces 6600-lb thrust 
from a 1900-lb engine. With an afterburner, the 101D will 
give close to 8500-lb thrust. A later modification of the 101D, 
the 101E is rated at 7270-lb thrust and can produce 9400-Ib 
thrust when equipped with an afterburner. 


SEVERAL reverse-thrust devices have been in the news 
lately. The French firm of Snecma has announced a reverse- 
thrust device that utilizes 20-30% of full-stage forward static 
thrust. The unit is claimed to have been test proved on 
five types of jet engines and that it has been used in place of 
air brakes up to Mach 0.9. Total running time of 150 hr 
is claimed. Boeing is said to be conducting extensive. re- 
verse-thrust tests for the all-jet B-52. The 150-ton aircraft 
is thought to require braking power considerably greater 
than that offered by more conventional methods (viz., 
wheel brakes and parachutes) to land on present airstrips. 
Although the work is classified, Boeing has reported “real 
progress’ in their reverse-thrust development. The de 
Havilland Comet, Britain’s all-jet transport, will have re- 
verse-thrust devices on the Model III series, according to a 
recent press note. De Havilland, like Boeing, is said to be 
conducting an extensive reverse-thrust development program. 
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PARACHUTE TEST AT SUPERSONIC SPEED 


8 SKOKIE MISSILE FOR TESTING PARACHUTE _— FROM 
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— Metals for Jet Engines 


A LARGE portion of metals research during the past 
year has been concerned with developing alloys for use at 
very high temperatures in corrosive atmospheres. Also 
pressing is the need for low-weight alloys of high strength in 
modern aircraft where weight must be kept at a minimum. 
Several developments are of interest in the jet engine field. 

New aluminum powder-metallurgy products, first intro- 
duced in Europe, have become available in the United States 
for limited experimental use. One such new alloy has a ten- 
sile strength at 900 F equivalent to the best existing high- 
temperature aluminum alloys at 600 F. It is hoped that this 
high-tensile strength feature can be yet extended. 

Production increases in jet engines have used large amounts 
of nickel; however, the amount of nickel per plane has been 
decreasing despite the change from cobalt to nickel alloys. 
Replacements for cobalt-type “superalloys” used in jet en- 
gines have not yet been found. 

One factor which has prohibited the widespread use of 
molybdenum for aircraft gas-turbine blades is the high oxi- 
dation rate. Meanwhile, such problems as production of 
large pieces, low-temperature ductility, fabrication, and good 
load-carrying features at elevated temperatures are reported 
as being solved or very near solution. Molybdenum has 
been roll-cladded with nickel alloys to provide flat-surface 
oxidation resistance at 1900 F. However, protection of edges 
and corners is still a problem in turbine blade applications. 
The production of ductile welds still remains largely unsolved. 

The family of stainless steel alloys has been continually 
increasing. As a substitute for the familiar 18-8 chromium- 
nickel stainless, a new stainless low-nickel austenitic steel 
with 16% chromium, 16% manganese, and 1% nickel was 
developed. As an extender, stainless-clad steel has come 
into wider use. 


ve Pilot plant production of a new titanium-tin alloy is ex- 
pected to have wide use in jet engines and airframes. The 
alloy consists of a titanium base with 2.5% tin and 5% alu- 
minum. Having exceptional strength characteristics up to 
1000 F, the alloy can be welded. For protection against 
corrosion, zirconium is added. Current production of tita- 
nium is about 2000 tons per year, while the Department of 
Defense is aiming at a production of 25,000 tons per year. 
Therefore, the present emphasis is on increasing production. 
Progress is being made in welding and production tech- 
niques. Alpha-phase alloys have given yield strengths of 
110,000 psi, while some alpha beta-phase alloys of 180,000 
psi yield strengths have been produced. 

Production of zirconium, with high-temperature and cor- 
rosion resistance, has been begun by several private firms. 
Current price of zirconium metal sponge is about $15 per 
pound. 

Much attention also has been given to the use of cermets 
for blade materials in turbojets. Titanium carbide base 
cermets have shown good properties at 1800 F. Main ave- 
nues of research have been to improve creep and oxidation 
qualities. The low ductility has remained a serious problem 


_in jet engine applications. 


Instrumentation and Equipment 


GENERAL Electric Company has announced development 
of its FC5 flight control system to be used in production mod- 
els of supersonic aircraft coming down the production line 
in 1955. Functioning as a fully automatic relief and ma- 
neuvering system, the FC5 consists of three basic units: a 
rate gyro, an amplifier, and an actuator. The rate gyro de- 
tects and accurately measures spurious attitude rates and 
sends signals to the amplifier which discriminates the signals 
and amplifies them to a level suitable for control of the actua- 
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tor. The actuator in turn moves the aircraft’s control sur- 
faces to oppose the unwanted attitude rates. 


THE extreme heat generated by jet engines and guided 
missiles calls for a new approach to temperature measurement. 
A promising new thermometer contains paint which changes 
color with variations in temperature. Scientists at the 
Naval Research Laboratory have been able to establish a 
series of compounds that will record temperatures from 50 
to278C. It is expected that future compounds will indicate 
temperatures in the higher ranges. 


THE temperature distribution in a rocket flame can be de- 
termined by a new gage developed at the General Electric 
Research Laboratory. Based on the use of a sodium-vapor 
lamp as a standard for comparison, the gage uses an inter- 
ferometer which magnifies the radiated spectrum coming 
from the lamp for a more detailed analysis. The new gage 
ean be used to survey the temperature structure in a large 
section of a flame, in contrast to earlier methods which 
permitted temperature measurements at one point only. 
Study of the rocket exhaust helps determine the efficiency of 
the rocket motor and tells how much energy is being con- 
verted into thrust. 


THE Hycon Manufacturing Company has announced pro- 
duction of its Submicrosecond Camera which enables expo- 
sure times of one ten-millionth of a second. There are no 
moving parts to the shutter which operates on the Kerr 
effect. An electric potential causes realignment of nitro- 
benzene molecules placed between two polarized screens, 
thus allowing light to pass onto the film. The camera is 
being used by Aberdeen Proving Grounds for ballistic re- 
search work. 


Reynolds Electrical 
and 


rr 


Engineering Co., Inc. 


Electrical and Construction = 


HOUSTON 

ALBUQUERQUE 

SANTA FE 
LAS VEGAS 


AN ORGANIZATION OF ELECTRICAL 7 
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yna-gage 


an electronic instrument for 
the measurement of static and 


dynamic pressures 


Versatile, positive, dependable — 
under the most adverse conditions 


m The water-cooled capacitance pick-up permits the measurement of pres- 
sures at extremely high combustion temperatures. Remote location of the 
pressure indicator is provided by cable connection, with lengths up to 100 feet 
performing satisfactorily. The pressure indicator can be installed with the dia- 
phragm flush with the wall of the pressure chamber, eliminating surge effects. 


Write or telephone for information 


PHOTOCON researcu PRobuCcTS 


Phone SYcamore 2-4131 


421 N. Foothill Blvd., Pasadena 8, Calif. e 
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standard instrument of 
the rocket industry . . . for 
indicating combustion pressures 


DYNA-GAGE DG-101 


Output: + 10 volts at 
25,000 ohms. 


Power requirements: 
PS-102. 
Size: 


836" x 1114" x 1114” 


CAPACITANCE PICKUP 
11-30. 


QOLEO PRESSURE GAS. 
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; eee which 
MPB ball hearing 
do you need? 


Most complete ever 


offered on miniature ball bearings 


© Complete specifications - 140 different types and sizes 

© Bearings from 1/10” to 3/8” o.d. shown in actual sizes 
e Speed-load charts with conversion factor 
e Lubrication government specifications, commercial sources 
e Recommended shaft, housing fits; shaft, shoulder data 
e Radial and axial clearance graphs 
e Typical methods of using miniature ball bearings 


For the designer of precision mechanisms 
this new 20 page MPB catalog offers practical 
solutions of problems involving miniaturiza- 

: tion. MPB has compiled for you the most com- 
plete and detailed information ever offered 
on this subject. Request the new MPB catalog 
53-54 on your letterhead . . . it may help 
develop a new product idea for you. 


| Miniature P,-ccision Bearings 
1corporated Kee New Hampshire 


INSTRUMENTATION 
SPECIALIST 


to be in charge of electronic instrumentation 
in connection with research and development 
of liquid propellants. Permanent salaried po- 
sition with expanding group in the research 
department of a major chemical company lo- 
cated in the Northeast. Requirements include: 
2 years of college plus 3 years of experience 
or equivalent. Experience with electronic cir- 
cuitry, high frequency transient pressure meas- 


urement, and some knowledge of liquid pro- 
pellants or chemistry desired. Send complete 
résumé to Box 717A, % Journal of the Amer- 
ican Rocket Society. 


Upper Atmosphere Research 


A PROJECT is being conducted by the USAF Air Research 
and Development Command to shed some light on the tre- 
mendous potential—upward of 100,000 volts—that exists 
between our earth and the upper stratosphere, a problem per- 
tinent to long-range communications today and one that may 
affect future extraterrestrial flight. Balloons are carrying 
aerial electrometers to altitudes of 100,000 ft above Hollo- 
man AFB, New Mexico. This upper stratosphere project 
is a part of a continuing study being made on the terrestrial 
electrical field existing between the earth and the ionosphere. 
For over 50 years, science has been trying to discover the 
source of the 1800-amp current constantly flowing toward 
the earth. 


THE “Jet Stream” is another upper atmosphere phenome- 
non being studied by the Air Force. The aim of ‘Project 
Jet Stream” will be to study the air current that ranges over 
the northern hemisphere at altitudes between 30,000 and 
40,000 ft, and sometimes reaches speeds of 300 mph. A 
B-47 and a B-29 will be used to measure winds, tempera- 
tures, humidity, and turbulence. Air Force scientists agree 
that the jet stream will be of considerable aid to both mili- 
tary and commercial planes that can take advantage of this 
super tailwind which travels from west to east. 


A NEW parachute jump record is claimed for the Navy 
pilot who recently bailed out of his all-jet F9F at 33,000 ft 
while cruising at 450 knots. Forced to abandon his craft 
because of an explosion his parachute opened after he had 
fallen 25,000 ft; he landed safely. 


FAST service in the overhaul and testing of fuel control 
systems for jet engines manufactured by Pratt & Whitney 
(J-48 and J-57) is now available to West Coast manufactur- 
ers. Pacific Airmotive Corporation will operate the jet 
testing oe at its Burbank, Calif., plant. 


Facilities ake 


THE El Segundo Division of Douglas Aircraft announced 
completion of a new $1.5 million hangar and flight test build- 
ing at Los Angeles International Airport to house final test 
installation work on production models of the bat-winged 
F4D Skyray. 


A NEW rocket facility to be constructed by the Navy De- 
partment and operated by Reaction Motors at Rockaway, 
N. J., was announced. The construction program will in- 
clude modern research and development laboratories, 4 
fabricating facility for experimental units, project engineering 
facilities, a cafeteria and supporting facilities, comprising 
about 150,000 sq ft. It is expected that the $3.5 million fa- 
cility will be completed in 1954. 


MACH number 13 has been attained in a wind tunnel 
designed by Cornell Aeronautical Laboratory for the Air 
Force’s Arnold Engineering Development Center. With a 
Mach 20 capability, the tunnel is of the hypersonic impulse 
type. The driving force of the tunnel comes from a detona- 
tion of a helium-hydrogen-oxygen gas mixture. The lu- 
minous gases may reach 7000 F during the 1-millisee run at 
a 10,000 mph test. 


BOEING Airplane has announced that it is operating the 
only privately sponsored wind tunnel capable of speeds in 
the transonic zone. This tunnel, originally designed for 
subsonic air speeds, was modernized by a completely new test 
section, a new fan, a new cooling system, and reinforced walls. 
Horsepower has been increased from 18,000 to 54,000. The 
new test section is 8 ft high and 12 ft wide. For tests at 
high speed it can accommodate models with wing spans as 


great as 9 ft. At lower speeds, larger models can be tested. 
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Regulus Lands After Tests... Although tactical missiles are 

intended for one-way flight only, considerable savings are real- 

ized in tests by being able to recover the particular unit. The 

landing here is monitored by a mother control plane and a photo- 

graphic jet as the Regulus opens its tail parachute to slow down 
the run at Edwards AFB, Muroc, Calif. 


Ejection Seat... This newest ejection seat is currently being 
installed in the F-94C Starfire jet interceptor. At left, the seat 
is in the down position and ready for ejection. By squeezing a 
lever on the right side of the seat activates the ejection propellant. 
At right, the seat is shown as it would appear in the process of 
being ejected. The guide rails steer the seat well above the 
top of the cockpit sill and give it a high trajectory after ejection 


Newest Jet Trainer... The latest jet trainer to be placed in 
operation is this North American TF-86, a two-place version of 
the Air Force F-86 Sabre Jet. The “hot” trainer is similar to the 
F-86 but has a tandem cockpit, dual controls, and is about 5 ft 
longer than the Korea-famed version. The speed is rated at 
over 650 mph, a service ceiling of 45,000 ft, and a combat radius _ ,Regulus Gets Rocket Boost ... Rocket boosters enable the 
of 600 miles. The trainer is to be used for advanced pilot train-s Regulus to get quickly into the air from a short ramp or launch- 
ing in high-speed flying, gunnery, and dive bombing é ing platform on a submarine or surface ship 


Aeroelastic Effects in Supersonic Flight ... These film strips, shot by automatic rocket-borne camera, illustrate the reaction of high- 

speed flight on different wing designs. The photographs were made in tests carried out by the Lockheed Aircraft Corporation. Strip 1 

shows a thin, straight wing in supersonic flight. Strip 2 shows a straight wing made deliberately flexible. The wing in strip 3 was a swept 

wing derived by sweeping back the straight wing in strip 2 and increasing the stiffness about 20 per cent. The reaction to high-speed 
flight was violent before failing 
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American Rocket Society News 


JOHN E. SCOTT, JR.,! GEORGE B. MATTHEWS,’ and ARTHUR A. KOVITZ? 


Introduction 


HE 8th Annual Convention of the 

American Rocket Society, held at the 
Hotel McAlpin in New York City from 
December 2 to December 4, 1953, again 
reflected the rapid growth of the Society 
with the presentation of a record number 
of twenty-one papers. The broadening 
interests of the Society were expressed in 
the wide diversification of the topics pre- 
sented. Fundamental research was repre- 
sented by the theoretical aspects of igni- 
tion and the aerodynamics of combustion; 
development activities were represented 
by such topics as the testing of large-thrust 
rocket propulsion systems; and the Space 
Flight Symposium provided a fitting cli- 
max to the Annual Convention. In addi- 
tion, the subject of the Honors Night Din- 
ner, Long Range Applications of Rocket 
Power, provided an opportunity for specu- 
lation as to the ultimate position of rocket 
propulsion in the power plant spectrum of 
the future. A large attendance represent- 
ing a nationwide interest assured the most 
effective exchange of technical informa- 
tion on a nonsecurity basis, both through 
formal and informal discussions. A good 
many of the papers were well worth listen- 
ing to, and some of them provoked spirited 
and informative discussion from the floor. 
Some other papers, however, did not seem 
to awaken the interest of the audience. 

It is to be hoped that in future years, 
more open discussion can be stimulated on 
controversial technical subjects, because, 
the interchange of ideas in this way can go 
a long way toward the solution of many of 
the serious problems in the field of rockets 
and jet propulsion. 


Session I 


Experimental Testing Techniques of Ex- 
tremely High Thrust Rocket Propulsion 
Systems, by R. F. Gompertz, Edwards 
Air Force Base. It was pointed out that 
the increasing requirements for extensive 
experimental testing of long-range guided 
missile rocket engines have severely taxed 
the capacity of available test facilities in 
this country. The requirements for sys- 
tem testing were outlined in which it was 
stressed that rocket test facilities be de- 
signed with the weapon system in mind; 
i.e., the test facility should be designed to 
reflect the ultimate mission of the rocket 
propulsion system. 

Several solutions to the problem of 
optimum test installations design which 


1 Research Assistant. Stud. Mem. ARS. 
2 Research Assistant. Mem. ARS. 
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Forrestal Research Center, Princeton University, Princeton, N. J. 


are currently under consideration were 
outlined by the author. An analysis was 
given comparing the simultaneous utiliza- 
tion of one single test cell for more than 
one large test article versus utilization of 
multiple noninterfering preparation areas 
supporting a single test stand. Test readi- 
ness and compartmentation of instrumen- 
tation, propellant lines and other test con- 
nections were major factors in this analy- 
sis. It was concluded that multiple-use 
test stands or preparation areas are the 
possible solution. At the conclusion of the 
paper, a color movie was presented show- 
ing the test facilities at Edwards AFB. 

A prepared discussion of Mr. Gompertz’ 
paper by R. F. Rose of the Jet Propulsion 
Laboratory, CalTech, was read. It ex- 
pressed the opinion that simultaneous 
occupancy of the test stand by more than 
one project should be rigorously avoided, 
since the possibility of mutual damage is 
too high a price to pay for the time saved. 
Development of proper handling gear and 
preparatory areas should eliminate the 
necessity for multiple occupancy. Mr. 
Rose felt that the paper was too super- 
ficial an analysis of the problem, and that 
more could have been said about test stand 
design. 


Multidirection Vibration Tester, by A. 
Bohr, Reaction Motors, Inc. Present-day 
military considerations, particularly in 
missile equipment installations, have led 
to the need for vibration testers capable 
of producing simultaneous oscillations in 
several directions over a wide range of 
amplitudes and frequencies. He described 
the construction of such a test system from 
commercially available parts in which 
electrically excited vibration shakers were 
used to vary the stem position in four-way 
valves and thus oscillate a load to a maxi- 
mum of approximately 50 cycles per second 
in both horizontal and vertical directions 
by means of hydraulic cylinders. During 
the discussion following the paper, it was 
pointed out that qualification test speci- 
fications for aircraft instruments require 
systems capable of producing frequencies 
as high as 500 cycles per second. 


Structural Design Considerations for a 
High Altitude Sounding Rocket of the 
Viking Type, by R.C. Lea, Glenn L. Martin 
Co. The need always exists for improve- 
ment in mass ratio by reduction of rocket 
empty weight. Among many problems 
are the evaluation of aerodynamic heating 
effects at presently attainable Mach num- 
bers above 4.0, and the development of 
quantitative design criteria for shock and 
impact loads. In addition, structural feed- 
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back of operating control frequencies to 
the control system by heavy oscillating 
masses appears to be one of the more diffi- 
cult problems of structural design. 


Session II-A 


Introduction to the Analysis of Supersonic 
Ramjet Power Plants, by G. Sears and B. 
March, Marquardt Aircraft Co. The 
effects of the various ramjet power plant 
components, altitude, Mach number, ete. 
upon ramjet engine performance were 
analyzed. The problems and limitations 
which determine the region of application 
of ramjet engines were also discussed. The 
equations necessary for the calculation of 
thrust and fuel consumption were derived. 
Simple charts were presented for the rapid 
estimation of ramjet engine performance. 


The 350,000 Pound Rocket Test Stand at 
Lake Denmark, N. J., by B. N. Abramson, 
D.S. Brandwein, and H. C. Menes, Naval 
Air Rocket Test Station. The test stand, 
whose basic structure is a steel cantilever, 
is constructed of steel and concrete. The 
firing attitude range in this stand is from 
vertical to horizontal with controls of 
operations being effected from the under- 
ground level in a blockhouse 250 feet away. 
The test stand itself is located on a steep 
hill such that the vertical distance from 
the firing mount to the ground is 60 feet. 
This provides protection from the effects 
of the exhaust jet. In addition, water for 
cooling the exhaust jet is provided. In- 
struments are installed to measure the 
usual rocket engine parameters, and 
enough channels are provided to allow for 
any additional measurements which might 
be needed for a complete engine test. 


The Problem of Cooling a Rocket Flame 
Deflector, by T. F. Reinhardt, Naval Air 
Rocket Test Station. The problem was 
that of cooling a rocket jet deflector plate 
which was to be exposed to the exhaust of 
a rocket having a possible maximum thrust 
of 100,000 Ib and located only 51 in. 
away. Several possible solutions to this 
cooling problem were considered, and four 
of these were discussed by the author. 
The methods considered were as follows: 
1 Utilization of the heat capacity and 
thermal conductivity of the plate to keep 
the surface temperature below the melting 
point. 2 Provision of a cooling water 
source beneath the deflector plate to ab- 
sorb the heat. 3 Flooding of the plate 
with water on its upper surface. 4 Cool- 
ing of the exhaust jet by water injection 
through high pressure nozzles. 

Each of these methods was explained in 
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jetail, and a comparison of the advantages 
and disadvantages of each was made. 
Experimental test results were presented 
to illustrate the major factors involved in 
ach method. 


Session II-B 


Some Aspects of Design and Fabrication 
of Liquid Propellant Engines, by H. 
Davies, Reaction Motors, Ine. Such 
topics as the choice of materials for par- 
ticular applications and the problems of 
valing against high temperature, high 
pressure gases were stressed in this paper. 
Strength-to-weight ratio, high and low- 
temperature strength properties, machina- 
ility, and compatibility with various 
propellants are involved in the selection of 
materials. The presentation concluded 
yith a comparison of various types of 
static and sliding seals now in use and a 
short discussion of the major factors which 
fect serviceability and ease of mainte- 
nance. A comment prepared by C. E. 
Hawk of the Aerojet-General Corporation 
emphasized the implication of the paper 
that rocket engines have reached a stage of 
development so advanced that designers 
concern themselves less with the problems 
of explosions or burnouts than with the 
more sophisticated questions of ease of 
production, durability, and field servicing. 
Mr. Hawk also noted the increasing im- 
portance of control systems and their in- 
tegration into over-all rocket motor designs. 


Rocket Preliminary Design and Proposal 
Preparation, by S. Lehrer, Reaction 
Motors, Inc. The author outlined the 
organization of his company’s proposal 
effort, with particular details on the inter- 
ution between various departments. He 
ilso discussed the mechanics of proposal 
writing and the structure of a special 
proposal report. 


Stability of Ethylene Oxide, by E. M. 
Wilson, Aerojet-General Corp. This 
paper was a supplement to several existing 
reports on the stability of ethylene oxide 
under normal operational hazards. It 
reported the results of detonation tests 
under conditions of sharp restriction in 
flow passages. Detonation did not take 
place in any of the pressurization tests. 
Similar stability and conspicuous absence 
of detonation was observed in the ejection 
tests. In an additional series of tests, 
ethylene oxide at temperature up to 150° 
was splashed against a stainless-steel 
block heated to temperature from 400 to 
1000° F. Again no detonations were 
present. However, experiments with ni- 
tromethane produced immediate detona- 
tions under similar pressurization and 
splashing conditions. An interesting dis- 
cussion followed the presentation in which 
such factors as time of contact with heated 
surfaces, contamination of such surfaces, 
ind incompatibility (when heated) of 
metals other than stainless steels were 
mentioned as important considerations in 
determining the stability of this com- 
pounc 


Session III-A 


The Isothermal Compressibilities of Some 
Rocket Propellant Liquids and the Ratios 


of Specific Heats, by G. Kretschmar, 
USNOTS, Inyokern. The author de- 
scribed a method for the determination of 
the isothermal compressibility in which a 
thin glass capsule containing the liquid 
was pressurized inside a thermostated 
glass vessel (or piezometer) that is filled 
with mercury up to a fiducial point. 
With a knowledge of the adiabatic com- 
pressibility, obta‘ned by the use of an 
ultrasonic interferometer, the author ob- 
tained the specific heat ratio simply by 
dividing the isothermal into the adiabatic 
compressibility. 


One-Dimensional Steady Adiabatic Flow 
in a Constant Area Channel with Mass 
Addition at Constant Enthalpy and Negli- 
gible Kinetic Energy, by E. W. Price, 
USNOTS, Inyokern. Five equations, 
including continuity, momentum, energy, 
heat capacity, and state are written in 
terms of the six variables, pressure, density, 
velocity, temperature, enthalpy, and mass 
flow. When one of these variables is 
assumed known, then the other five may 
be determined. In this paper, the inde- 
pendent variable was chosen as the mass 
flow. Finally, application of the results to 
a solid propellant rocket was discussed 
briefly. 


Investigation of Boiling Heat Transfer 
and Burnout to JP-4, by C. M. Beighley and 
L. Dean, Bell Aircraft Corp. A _ trans- 
parent annular test section surrounding 
an electrically heated stainless-steel tube 
yas used to obtain heat transfer data in 
the nonboiling, nucleate boiling, and film 
boiling regions, heat flux values being in- 
creased until the test section burned out 
or the limit of the power generator was 
reached. The nonboiling heat transfer 
data were correlated by the Seider-Tate 
equation with a constant of 0.0213 instead 
of 0.023. It was observed that a sand- 
blasted tube caused approximately 50 per 
cent increase in the heat transfer coefficient 
over the regular and polished tubes in the 
nonboiling region. The nucleate boiling 
range was investigated quite thoroughly 
up to the point where film boiling begins 
and the wall normally fails. In the film 
boiling region, the authors found that it 
was possible to obtain a 200 per cent in- 
crease in heat flux with less than a 20 per 
cent increase in wall temperature. Fi- 
nally, the formation of coke on the heater 
tube surface was discussed. It was con- 
cluded that coke deposition did not occur 
below a given heater tube surface tem- 
perature, that temperature depending 
upon the constituents in the fuel. 


Session III-B 


Study of the Combustion of Fuel Droplets 
Descending Through an Oxidizing Atmos- 
phere, by D. Charvonia, Purdue Univer- 
sity, and H. Wood, Virginia Polytechnic 
Institute. The studies were made with 
three organic fuels, namely, triethyl- 
amine, allylamine, and cyclohexene, de- 
scending through decomposition vapors of 
white fuming nitric acid. Interruption of 
a light beam by the falling droplet starts 
operation of a Fastax camera photograph- 
ing the descent of the droplet. It was 
found from these photographic records 
that ignition is initiated in the vapor 
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region in the wake at a considerable dis- 
tance behind the droplet, producing a 
flame front which propagates at a speed 
greater than droplet velocity, and hence 
overtaking the droplet. Imgnition delays 
were measured for a wide range of condi- 
tions for each of the three fuels men- 
tioned. Increasing vapor temperature de- 
creased the ignition lag of two of the fuels, 
leaving the third relatively unaffected; a 
fuel temperature increase showed a tend- 
ency to decrease the ignition lag; and fuel 
droplet size was found to have a marked 
effect on ignition delay, appearing to give 
a minimum ignition delay at some definite 
drop size. Attempts to explain these varia- 
tions from a theoretical approach have 
thus far proved unsuccessful. 


Ignition and Combustion in a Laminar 
Mixing Zone, by F. E. Marble and T. C. 
Adamson, California Institute of Tech- 
nology. This paper considered the theo- 
retical approach to the problem of igni- 
tion and combustion in a laminar mixing 
zone between two parallel streams, one of 
which consists of a cool combustible mix- 
ture and the other being composed of hot 
combustion products. After discussion of 
the physical and chemical relations which 
apply to the problem, the initial portion 
of the mixing zone was investigated where 
the heat evolved through chemical reac- 
tion is small, and the problem is solved by 
consideration of a perturbation to the 
solution for mixing without combustion. 
This gives the analytical picture of the 
early portion of the mixing zone. 

The complete course of development of 
a laminar flame front was then examined 
by means of an approximate treatment 
which effectively extended the solution 
throughout the remainder of the flame 
zone. This analysis indicated that a new 
class of combustion problems is open to 
investigation through application of the 
usual boundary layer concepts. 

In a prepared comment, J. A. Bierlein 
pointed out that the mathematical de- 
scription is in harmony with experience. 
He also pointed out the benefits which 
might be forthcoming from an extension 
of the analysis to consider fractional and 
second-order reactions in an attempt to 
obtain a mass of data depicting flame 
stand-off distance over a wide range of 
values for all the disposable parameters. 
It might then be possible to predict the 
complete blow-off or quenching behavior 
of a system on the basis of a few experi- 
mental results, without prior knowledge of 
the chemical kinetics or the transport 
properties of the system. The discussion 
from the floor which followed indicated a 
widespread interest in this paper. 


Erosive Burning of Some Composite 
Solid Propellants, by L. Green, Jr., Aero- 
jet-General Corp. The paper reported ex- 
perimental data on the erosive burning 
characteristics of several composite solid 
propellants. It was found that the data 
were best correlated in terms of a reduced 
mass velocity. In discussing the existence 
of an apparent “threshold velocity for ero- 
sion,” the author pointed out that this is 
probably due to the conditions at the fore 
end of the grain. In order to determine the 
best method of correlation and to resolve 
the question as to the existence of the 
threshold velocity, it was pointed out that 
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experiments of increased precision are 
necessary. 

In a prepared comment by R. N. Wim- 
press, read by the vice-chairman, it was 
pointed out that the existence of a thres- 
hold velocity has been definitely shown for 
JPN and ballistite and must be considered 
in any design using these propellants. It 
was also mentioned that the chamber 
pressure and the gas velocity must be con- 
sidered as being two separate effects on the 
“erosion ratio” e. 


Session IV-A 


Oxidant Pumps, by W. Mizen, Bendix 
Aviation Corp. The results of develop- 
ment work on a positive displacement 
plunger-type pump for 90 per cent H,O, 
were described. A single plunger test 
pump was built using Type 440-C steel, 
bright hardened to 55 R, for the plunger 
housing and Type 416 steel chrome- 
plated for the plunger. All other mechani- 
cal components were made of 18-8 stain- 
less steel with neoprene as the O-ring seal 
material. This pump performed well, but 
failure after 9 hr of operation was caused 
by severe pitting of the Type 440-C dis- 
charge valve seat. To extend the life of 
the pump, materials such as glass, syn- 
thetic sapphire, and dry molded alumina 
ceramic were tried for the plunger housing 
and plunger materials. The test results 
which were reported evoked active dis- 
cussion on the relative merits of all kinds 
of materials for such service. 

Two prepared comments on this paper 
were read. The first comment, by D. 
Chatfield of RMI, complimented Mr. 
Mizen on providing a good insight into the 
use of some rather unusual materials in 
notching or moving parts. With respect 
to the use of ceramics or high strength 
glass as bearings, he felt that this applica- 
tion would be limited to moderate speeds 
and thrust loads. He also verified that the 
use of steels of the 400 series is usually 
avoided because they cannot be success- 
fully passivated. 

The second comment was prepared by 
David K. Hart, USNOTS, Pasadena, 
Calif. In discarding conventional corro- 
sion-resistant materials, Mr. Hart felt 
that too much haste was shown. He 
pointed out that pitting occurred only in 
areas which were not critical. If the fail- 
ure was due to pitting of the discharge 
valve seat, it seems reasonable to suppose 
that the installation of a Teflon valve seat 
would have solved the problem without 
resorting to the longer term investigation 
of ceramics. 


An Analytical Procedure for Mixed Acid, 
by J. Clark, Naval Air Rocket Test Sta- 
tion. The problem of determining the 
composition of mixed nitric acid, a liquid 
that is in common use as an oxidizer for 
rockets, was described. Mixed acid is 
approximately 83 per cent HNO, 15 per 
cent H2SO,, plus small amounts of water 
and NO:. The method in common use is 
based on three titrations, a procedure de- 
veloped by Scott and modified at the 
author’s laboratory. This method is 
entirely satisfactory except for the H2SO, 
titration. The author’s line of attack re- 
solved itself into a conductimetric end- 
point method for the H2SO, after most of 
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the interfering HNO; had been removed. 
In conclusion, it was stated that the 
method is probably reliable to 0.1 per 
cent H.SO,, which is a considerable im- 
provement over the current procedure. 


Experience with the Application of Hydro- 
gen Peroxide for Production of Power, by 
H. Walter, Worthington Pump and 
Machinery Corp. The paper constituted 
a recounting of personal experiences in 
Germany with hydrogen peroxide in appli- 
cations for the production of power. This 
historical account began in 1933 with a 
proposal for a submarine engine using 
hydrogen peroxide, and extended through 
World War II in which peroxide rocket 
engines were used to propel aircraft and 
guided missiles. The major emphasis of 
the paper was on the reliability and rela- 
tive hazards involved in the handling, 
storage, and power applications of per- 
oxide, and a number of interesting experi- 
ences were related to illustrate these prop- 
erties. After a period of research and de- 
velopment on the sensitivity of hydrogen 
peroxide and mixtures of this oxidizer and 
various organic fuels, a number of per- 
oxide-drive devices were tested, accepted, 
and introduced into service without a 
single accident directly attributable to 
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Low-Frequency Combustion Stability of 
Liquid Rocket Motor With Different 
Nozzles, by Sin-I Cheng, Princeton Uni- 
versity. The paper reported the results 
of a theoretical investigation of the effects 
of the nozzle flow on the stability of low- 
frequency oscillations in liquid propellant 
rockets. The fundamental equations were 
developed on the basis of assumptions out- 
lined in previous analyses of low-frequency 
instability. The resulting characteristic 
equation was solved to determine the 
boundaries of stable operation. These new 
stability boundaries thus represent a gen- 
eral solution for all rockets with a given 
combustion chamber and feed system, but 
with different Laval nozzles. It was 
generally concluded that a decrease in con- 
traction ratio with constant length of con- 
verging section or an increase in length 
with constant contraction ratio would 
have a stabilizing effect on low frequency 
oscillations. 


Low-Speed Combustion Aerodynamics, by 
R. A. Gross and R. Esch, Harvard Uni- 
versity. The authors introduced the 
problem of the two-dimensional low-speed 
steady flow field containing a finite length 
discontinuity which represents a flame in 
unbounded space with uniform parallel 
flow at infinity. By making the assump- 
tions that the density ratio and pressure 
ratio across the flame are constant, the 
heat of reaction is constant and the flow 
field is everywhere irrotational, the flow 
field can then be described by classical 
potential theory where the flame discon- 
tinuity is represented by a source sheet. 
Two approximate solutions to flow fields 
were presented. The first was a straight 
line flame of finite length inclined at an 
angle to the free stream, and the second 
was the V-flame shape. The experimental 
data showed good agreement only for 


small values of the ratio of flame speed to 
the uniform velocity at infinity. 

The general discussion which followed 
the presentation of the paper dealt mainly 
with the validity of the basic assumptions 
and their possible verification experi- 
mentally. On the basis of the experi- 
mental data obtained by the author, it was 
pointed out that the source sheet is a 
powerful tool for use in solving problems 
of this general type. 


Instantaneous Rocket Flame Temperature 
Measurement, by J. H. Hett and J. B. 
Gilstein, New York University. The 
paper described some measurements of 
instantaneous combustion chamber and 
exhaust temperatures in a 1500-lb thrust 
alcohol-oxygen rocket. The authors gave 
an account of their work while at the GE 
Malta Test Station from January 1950 to 
August 1951. The temperature instru- 
ment used was an instantaneous two-path 
optical pyrometer. Since the continuous 
background of the spectrum of a hydro- 
carbon flame is believed to be in thermal 
equilibrium with the gas, this instrument 
uses for temperature measurement a 
region of this spectrum that is free from 
line and band structure. Planck’s radia- 
tion law and the usual negative exponen- 
tial function for absorption are used for 
the purpose of analysis. Spectrographic 
data showed that the radiation was not 
quite Planckian so that the recorded tem- 
peratures should be considered only as 
“indicated” temperatures. Temperatures 
of about 2130° K at the first Mach cone 
were shown in contrast to the Bundy and 
Strong figure of 2360° K, using a sodium- 
line reversal method. The _ indicated 
average temperature was found to be 
1900° K=+50° with an emissivity of 
0.15+0.05. A spectrographic study of the 
radiation from the chamber through a 
quartz window was undertaken to deter- 
mine the departure from gray-body radia- 
tion. Because of this departure, the au- 
thors were again obliged to speak of “indi- 
cated’’ temperatures. An average tem- 
perature of 2750° K is presented, whereas 
the calculated theoretical maximum ten- 
perature for the mixture ratio used was ap- 
proximately 3200° K. 


Epitor’s Note: The papers summarized 
in this review will appear in future issues of 
JET PROPULSION as space becomes available. 
In the meantime, preprints may be obtained 
from the ARS Secretary at 25 cents per copy. 


Space Flight Symposium 


HE final session of the ARS 1953 

Convention was a symposium on space 
flight held in the Ballroom of the Hotel 
McAlpin on Friday morning, December 4. 
This symposium took the form of a panel 
discussion followed by a formal question 
period. The members of the panel were: 
Andrew G. Haley, chairman, of Haley, 
Doty, and Schellenberg; John R. Pierce of 
Bell Laboratories, Inc., main speaker; 
J. H. DeWitt, Jr., of Stations WSM and 
WSM-TV; Krafft A. Ehricke of Bell Air- 
craft Corp.; E. E. Francisco, Jr., WSPG; 
R. W. Porter, General Electric Co.; Dar- 
rell C. Romick, Goodyear Aircraft Corp.; 
Milton W. Rosen, Naval Research Labo- 
ratory; Kurt R. Stehling, Forrestal Re 
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search Center, Princeton University; A. 
R. Tocco, Hughes Aircraft Co.; R. C. 
Truax, U. S. Navy Bureau of Aeronautics; 
and Paul F. Winternitz, New York Uni- 
versity. 

Each of the panel members presented a 
short talk on a particular phase of space 
travel, after a brief introduction by the 
chairman. Col. DeWitt discussed com- 
munications between the earth and a 
satellite vehicle, mentioning the possible 
use of such a vehicle as a television trans- 
mitting station and demonstrating the ease 
of such transmission by the relative success 
of radar communication to the moon. 

Mr. Ehricke spoke about a new system 
for supplying orbital vehicles by means of 
winged rockets, and discussed the problem 
of calculating the optimum position and 
direction of take-off from an orbit for in- 
terplanetary travel. 

Mr. Francisco commented on the role of 
research in advancing the science of space 
flight. He emphasized the need for more 
fundamental work in high altitude re- 
search and mentioned the contributions 
made in this field by the White Sands 
Proving Ground. 

Mr. Romick discussed control systems 
for space ships, stressing the importance 
of guidance. He suggested automatic 
computers feeding servo control systems 
as a possible method for flight correction 
of initial directional deviations. 

Mr. Rosen then spoke on the topic, 
“Margin for Error,” in which he pointed 
out a factor usually overlooked in satellite 
or interplanetary vehicle calculations, 
namely, the extreme effects which small 
errors or variations in specific impulse or 
propellant flows can produce in final veloc- 
ity and altitude. He applied a reasonable 
value of propellant ‘‘outage’’ to the cal- 
culations of three satellite designs; and 
showed how far short of predicted perform- 
ance they might fall. 

Mr. Stehling spoke on the hazards of 
radiation in an orbital vehicle. He demon- 
strated how secondary emissions might 
produce much more harmful effects than 
the primary cosmic or solar radiation, and 
he suggested the possibilities of using water, 
paraffin, or electromagnetic force fields 
to shield human beings from such rays. 

Mr. Tocco remarked on some promising 
new materials for rocket vehicles with re- 
gard to the problems of overheating during 
re-entry. He mentioned the use of ti- 
tanium, porous metals, and ceramics as 
examples. 

Dr. Winternitz gave a brief discussion 
on the improvement of specific impulse as 
a means of obtaining orbital velocities 
from multistage chemically fueled rockets. 
He was followed by Cmdr. Truax who 
made a plea for better economic analysis 
of satellite vehicle designs, illustrating the 
importance of such considerations by 
analyzing the relative costs of returnable 
and expendable third stages of a missile of 
the type suggested by von Braun. 

Dr. Porter then outlined the policy of the 
American Rocket Society toward space 
fight, accentuating the importance of 
interplanetary exploration as a scientific 
challenge, independent of any military or 
practical applications. 

Dr. Pierce, the main speaker, then ad- 
dressed the audience on the over-all 
problem of communication in space. His 
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analyses and conclusions were based en- 
tirely on existing equipment and tech- 
niques and indicated that, in general, 
transmission and receipt of information 
between Earth and a moving space 
vehicle or between Earth and a fixed inter- 
planetary body would pose fewer prob- 
lems than similar communication be- 
tween two points on the earth’s surface. 
He proposed 3-centimeter wave-length 
radiation and analyzed the power require- 
ments, size of equipment, signal-to-noise 
ratio, and velocity corrections involved in 
transmitting or receiving such waves on a 
moving vehicle outside the earth’s atmos- 
phere. On the topic of Earth-Moon com- 
munication, Dr. Pierce was even more 
optimistic, declaring that telephone and 
telegraph contact between these two bodies 
would be possible with power supplies of 
only a few hundred watts. In conclusion, 
he speculated that if interfering noise can 
be reduced and power-to-size ratio can be 
increased at their present rates of develop- 
ment, communication with the stars may 
not be unreasonable. 

At the close of this formal part of the 
session, an interesting question period was 
held. Written questions were submitted 
by those in attendance, addressed to in- 
dividual panel members. Such topics as 
the use of solid propellants for space 
rockets, the minimum skin temperatures 
possible in descent through the atmos- 
phere, the use of powdered metals for pro- 
pellants, video communication from Earth 
to Mars, the technical definition of 
gravity, and the use of nuclear power 
plants in rockets were treated by the 
various speakers. 


Chemical Industry’s Stake 
in Rockets 


 & EDWARD Pendray, ARS Director, 

spoke at a dinner meeting of the 
Chemical Public Relations Association on 
February 3 at the Hotel Roosevelt in 
New York. 

After delivering a talk on rocket history 
and the beginning of the ARS, Pendray 
told the group of some 60 public relations 
men from 30 major chemical companies 


MRS. ESTHER C. GODDARD, WIFE OF THE LATE DR. ROBERT H. GODDARD, 
LUNCHEON TENDERED HER BY BRIG. GEN. GEORGE G. EDDY, 


MISSILE TEST STATION, WSPG; 


what rocket development could mean 
to their industry. He contrasted the 
“early days’? when liquid oxygen used 
to be “swiped by the container” to today 
when it moves into missile test centers 
by the tank carload. He also pointed out 
that 29 of the 79 technical papers pre- 
sented in the ARS JourNat since Sep- 
tember, 1951, have dealt directly or in- 
directly with chemicals, as well as 31 of 
111 papers presented at technical meetings 
since January, 1949. 


ARS in Dallas Airpower 
Symposium 


AMONG the speakers at the recent Air- 
power Symposium sponsored by the 
Dallas Council on World Affairs was 
Frederick C. Durant III, last year’s ARS 
President and current President of the 
International Astronautical Federation. 
Durant’s talk—on space flight—followed 
speeches delivered by General Nathan F. 
Twining, USAF Chief of Staff; James H. 
Smith, Jr., Under Secretary of the Army; 
Earl D. Johnson, and Arthur Godfrey. 


‘Dream Come True”’ Says 


Mrs. Goddard at WSPG 


66" HIS rocket range would have been 

a dream come true for my husband 
if he had lived to see it,’’ remarked the 
widow of the late Dr. Robert H. Goddard, 
during her recent visit to White Sands 
Proving Ground. 

In New Mexico for the opening of an 
exhibit of some of her husband’s work at 
Roswell, Mrs. Goddard was honored 
guest and speaker at a banquet arranged 
by the New Mexico-West Texas Section 
on January 22. She spoke on the “human 
aspects” of Dr. Goddard’s work from 
1925 to 1941 in New England and New 
Mexico and also showed motion pictures 
of static and flight tests of the rocket 
pioneer’s experimental units, fired at 
Eden Valley near Roswell. 

During her visit to WSPG, Mrs. God- 
dard was also tendered a luncheon by 
Brig. Gen. George G. Eddy, commanding 
general of WSPG. 


IS SHOWN AT 
COMMANDING GENERAL OF 
WSPG. AT RIGHT IS CAPTAIN PHILIP D. QUIRK, COMMANDER OF THE U.S. NAVY ORDNANCE 


E. E. FRANCISCO AND FRANK L. KOEN, JR., PAST-PRESIDENT 


AND PRESIDENT, RESPECTIVELY, OF NEW MEXICO-WEST TEXAS SECTION ARE AT LEFT 
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RESIDENT Andrew G. Haley served 

as chairman of a technical session on 
rocket propulsion at the IAS Twenty- 
Second Annual Meeting on January 29 
in New York. 

Four papers were presented before an 
audience of about 60 people at the session, 
which was arranged by ARS. 

“Ballistics of an Evaporating Droplet” 
by C. C. Miesse of Aerojet-General Cor- 
poration, Azusa, Calif., analyzed the be- 
havior of droplets—particularly fuel drop- 
lets injected upstream in a ramjet burner— 
when droplet velocity and size, air velocity, 
and physical properties of the surrounding 
air are taken into account. Making sev- 
eral basic assumptions, Miesse pointed 
out the possibility of determining analyti- 
cally the velocity profile of an evapo- 
rating droplet. 

“Gas Torch Igniter for Rocket Motors” 
by George N. Woodruff of Reaction 
Motors, Inc., Rockaway, N. J., described 
an ignition system which promises, ac- 
cording to Woodruff, greater reliability 
in firing of rocket engines over a wide range 
of thrusts. Hydrogen and oxygen were 
the gases used principally during tests 
on which the paper is based. 

In “Toxicity and Health Hazards of 
Rocket Propellants,’’ Stephen Krop of the 
Army Chemical Center, Maryland, gave 
a detailed analysis of the characteristics 


and potential handling dangers of rocket 


fuels, oxidizers, and their decomposition - 
products. Included were ammonia, ani- 
line, hydrazine, JP-4, hydrogen peroxide, 
liquid oxygen, red fuming nitric acid and 
ethyl, methyl, and furfuryl alcohol. 

The final paper, “A Light Weight Elec- 
tric Analog Computer With no Moving 
Parts,’ was presented by Paul H. Savet of 
Arma Corporation, Mineola, N. Y. It 


New York Section to Spon- 
sor Education Program 


Or MARCH 24 in New York City, 

a group of outstanding university 
scientists will discuss the engineering prepa- 
ration advisable for students looking 
forward to participation in the field of — 
space flight. 


The panel, to be moderated by Dr. Mar- 


tin Summerfield, Editor-in-Chief of Jet 
Proputsion and Professor of Jet Propul- 
sion at Princeton University, will include 
Dr. Charles L. Mantell, Chairman of the 
Department of Chemical Engineering, 
Newark College of Engineering; Prof. 
Charles H. Townes, Executive Officer, 
Department of Physics, Columbia Uni- 
versity; Prof. William R. MacLean, De- 
partment of Electrical Engineering, Poly- 
technic Institute of Brooklyn; Prof. Frede- 
rick K. Teichmann, Chairman of the 
Department of Aeronautical Engineering, 
College of Engineering, New York Uni- 
versity; and Dr. Alfred Bornemann, Pro- 
fessor of Metallurgy, Stevens Institute of 
Technology. 

Undergraduate engineering and science 
students from colleges in the metropolitan 
area, as well as interested high-school 
science students will be invited to the 
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ANDREW G. HALEY A 


Haley Presides Over ARS-IAS Session 


THE ARS-IAS 


JOINT MEETING 


At. 
desctibed a device which performs basic 
calculus operations directly on a-c signals 
and multiplies and divides statically by 
means of heat transfer, eliminating the 
need for servomechanisms. The mini- 
ature computer was developed principally 
for applications in navigation and piloting 
of aircraft and guided missiles. 


meeting which will be held at the American 
Museum of Natural History at 8:00 p.m. 
Arranging the program for the New York 
Section is C. W. Chillson, Past-President 


of the ARS. 
2 
IAF is to Award 
Loeser Medal 


HE International Astronautical Fede- 

ration announces that the subject 
of the competition for the first annual 
Guenter Loeser Memorial Medal will be 
“On the Economics of Space Flight.” 

The award, honoring Dr. Guenter 
Loeser, former IAF vice-president killed 
in 1953 ‘in a plane crash doing consulting 
work for the U. 8. Air Force, was estab- 
lished at last year’s IAF Convention. 
Entries should be submitted to the Secre- 
tary, International Astronautical Federa- 
tion, P. O. Box 37, Baden, Switzerland, 
in five copies, not later than December 31, 
1954. They should be forwarded in two 
envelopes, the outer marked ‘Loeser 


Medal Competition.”” The inner (contain- 
ing the entry) should be blank as it will be 
inscribed with a code mark in place of the 
author’s name. 


Four More Corporate 
Members 


outstanding companies in the 

guided missile and aircraft manu- 
facturing fields, one chemical company, 
and an aviation publishing company are 
recent additions to the ARS Corporate 
Member lists. 

They are Lockheed Aircraft Cor- 
poration, Burbank, Calif., which just an- 
nounced the formation of a Missile Sys- 
tems Division; Bell Aircraft Corporation, 
Buffalo, N. Y., designers and producers 
of liquid propellant rockets and rocket 
engines; Linde Air Products Co., Division 
of Union Carbide and Carbon Co., New 
York, N. Y., makers of liquid oxygen; 
and Aeronautical Digest Publishing Co., 
Washington, D. C., publishers of Aero 
Digest and the Aircraft Yearbook. 


Going to Innsbruck? 


Reservation forms are being circu- 
lated to all ARS Members for a 21-day 
tour to the 5th Annual Convention of 
the International Astronautical Fed- 
eration at Innsbruck, Austria. 

The tour will leave by plane from 
New York on July 25 and will include, 
tentatively, four days in Paris, three 
days in Switzerland, six days in Inns- 
bruck for the IAF Convention, August 
1-7 (with side trips optional to such 
places as Bolzano, Italy, Oberam- 
mergau, Garmisch-Partenkirchen), four 
days in Germany (including a trip up 
the Rhine), and two days in Belgium. 
| Return to New York will be on August 
15. 

Cost of the trip to ARS members 
and their guests will be $750 all-inclu- 


sive (except for extras such as drinks 
and gifts). Seats will be assigned as 
reservations are received by the travel 
Details are included in the 


| agency. 
form. 
Any member who has not received a 
| reservation form can obtain one by 
writing to the national office. 


ARS Fellow Honored 


N FEBRUARY 8, Dr. Theodore von 
Ka4rma4n was honored at a reception 
at the Italian Embassy in Washington, 
D. C. Ambassador Alberto Tarchiani, 
on behalf of the Italian Government, pre- 
sented Dr. von K4rm4n with the Order “‘Al 
Merito della Repubblica” in recognition of 
his scientific achievement and contribu- 
tions to the field of aero- and thermody- 
namics. In 1949, ARS honored Dr. von 
Karman by awarding him its Fellow mem- 
bership 
Mrs. Andrew G. Haley (in ARS Presi- 
dent Haley’s absence) and IAF President 
F. C. Durant III attended the reception. 
Others also attending included Asst. 
Secretary of Defense Donald A. Quarles; 
USAF Chief of Staff General Nathan F. 
Twining; and USAF Asst. Chief of Staff 
Lt. Gen. L. L. Craigie. 
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> Section Doings 


Detroit. Results of the election of 
officers for this recently organized group 
are as follows: Laurence M. Ball, Presi- 
dent; Fred Klemach, Vice-President; 
David Buell, Secretary; Dr. Richard B. 
Morrison, Treasurer. On the Board of 
Directors are Dr. Arthur B. Ash, Charles 
W. Williams, Dr. Leland G, Cole, and Dr. 
Arthur A. Locke. 

Florida. On January 26 more than 
350 persons listened to a talk on space 
flight by Wernher von Braun of Redstone 
Arsenal. This was the second public 
mecting of the newly formed group which 
is centered at Patrick Air Force Base. 
The first, held on November 20, 1953, had 
ARS Director Dr. Richard W. Porter as 
principal speaker. Officers of the Sec- 
tio: are Keith K. McDaniels, President; 
Major Prentice B. Peabody, Vice-Presi- 
dent; and 1st Lt. Theodore A. DeBacker, 
Secretary-Treasurer. 

National Capital. Officers recently 
elected are: James R. Patton, Jr., Presi- 
dent; Esterly C. Page, Vice-President; 
Ethna White, Secretary; Richard K. Neu- 
main, Treasurer. Directors are Harry 
J. Archer, J. D. Gilchrist, Andrew G. 
Haley, Charles F. Marsh, Dr. S. F. Singer, 
and Richard B. Snodgrass. 

The Section gave a dinner in honor of 
Dr. Theodore von Kaérmdén on Friday, 
January 29. There were approximately 
60 members in attendance. Andrew G. 
Haley, ARS President, discussed his pro- 
gram for the coming year. Dr. von Ka4r- 
man gave a very interesting talk on the 
future of rockets and included some in- 
formative comments on foreign develop- 
ments. 

On Friday, February 12, a two-hour 
film, ‘History of Rocket Development,” 
was shown to approximately 100 members 
and guests at the second meeting of this 
new section. The film, released through 
the British Ministry of Supply, traced the 
development of rockets in Germany from 
the very early work of Opel, Valier, Hey- 
landt, and Zucker to the Rheinbote, 
Rheintochter, 


Wasserfall and A.4,. Of 


JAMES J. WYLD 


Announcement was made at the 
ARS 8th Annual Honors Night Dinner 
on Dec. 3, 1953, of the death of James 
J. Wyld, pioneer developer of rocket 
power, past president of ARS, and one 
of the founders of Reaction Motors, 
Inc., Rockaway, N. J. 

Mr. Wyld, 41 years’ old, died of a 
heart ailment on Dec. 1 at his home in 
Pompton Lakes, N. J. A graduate 
of Princeton University where he had 
achieved a brilliant scholastic career, 
he helped to form RMI in 1941 after 
having done research work on high 
speed aerodynamics for NACA for a 
time. His work at RMI established 
him as one of the leading figures in the 
rocket field and much of his research 
and development work has seen frui- 
tion in liquid propellant engines which 
power today’s aircraft and guided 
missiles. 

A member of the Board of Directors 
of RMI at his death, Mr. Wyld had 
served as a consultant to the Atomic 
Energy Commission in recent years. 
He leaves his wife, Mrs. Helen Wyld, a 
daughter, Marie, 4, and a son, Robert, 
2. 


particular interest were combat shots of 
the Me.163 rocket interceptor attacking a 
formation of Flying Fortresses. The film 
closed with the take-off of the space rocket 
from Fritz Lang’s Frau Im Mond. 

New Mexico-West Texas. The first 
issue of Missile Away, published in 
December, 1953, was so successful that 
it is planned to print it quarterly. This 
year’s officers are Frank L. Koen, Jr., 
President; Edward L. Brown, Vice- 
President; and G. Harry Stine, Secretary- 
Treasurer. The Board of Directors will 
consist of E. E. Francisco, Jr., Chairman, 
Clyde Tombaugh, Dr. Russell K. Sher- 
burne, Charles W. Mansur, John D. Pa- 
trick, R. Gilbert Moore, Fred Koether, 
and John S. Piech. Honorary directors 
are Brig. Gen. George G. Eddy, com- 


DR. RICHARD W. PORTER, ARS VICE-PRESIDENT, DISCUSSES FLORIDA SECTION 


ACTIVITIES PRIOR TO A MEETING HE ADDRESSED RECENTLY. 
ARE MAJOR PRENTICE B. PEABODY, VICE-PRESIDENT; 
DR. PORTER; 


DEBACKER, SECRETARY-TREASURER; 
PRESIDENT OF THE SECTION. 
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LEFT TO RIGHT 
IsT LT. THEODORE A. 
AND KEITH K. MCDANIELS, 


manding general of White Sands Proving 
Ground, and Capt. Philip D. Quirk, 
commander of the U. S. Navy Ordnance 
Missile Test Station, WSPG. 

New York. New officers are Michael 
J. Samek, President; Thomas Reinhardt, 
Vice-President; W. Dorwin Teague, Jr., 
Secretary; and Robert J. Thompson, Jr., 
Treasurer. 

On January 29 William Wright of 
Reaction Motors, Inc., talked to a group 
of about 40 members and guests at the 
Engineering Societies Building on “Patents 
and the Rocket Engineer.” A film from 
the Glenn L. Martin Company on the 
Viking #1 and #2 firings was shown. 

Northern California. C. C. Ross of 
Aerojet-General Corporation, Azusa, spoke 
on January 14 on “Combustion Insta- 
bility in Rocket Motors” at a dinner meet- 
ing in San Francisco. New officers are 
William J. Barr, President; M. A. Pino, 
Vice-President; John S. Amneus, Secre- 
tary; and Carl W. Messinger, Jr., Treas- 
urer. Directors are A. J. Eggers, Marjorie 
Evans, S. 8S. Sorem, Austin Bryant, and 
A. Quarterman. 

Southern California. Dr. Leo Stool- 
man, a member of the propulsion research 
group at Hughes Aircraft Company, Cul- 
ver City, spoke on ‘Ramjets” at the 
December 17 meeting. New officers are 
H. S. Seifert, President; C. M. Me- 
Closkey, Vice-President; W. D. Crossland, 
Secretary-Treasurer. The Board of Di- 
rectors includes 8S. K. Hoffman, D. I. 
Baker, T. F. Dixon, R. D. Geckler, R. F. 
Gompertz, B. T. Morris, and W. L. Rogers. 


ARS SECTION PRESIDENTS 


Alabama GEORGE HENDERSON 
Redstone Arsenal 

Arizona ANTHONY R. TOCCO 
Hughes Aircraft Co. 

Chicago KENNETH H. JACOBS 


American Machine & 
Foundry Co. 
Cleveland-Akron 
DARRELL C. ROMICK 
Goodyear Aircraft Corp. 


Detroit LAWRENCE M. BALL 
Chrysler Corp. 
Florida K. K. McDANIEL 
Boeing Airplane Co. 
Indiana J. M. BOTJE 
Purdue University 
Maryland IVAN E. TUHY 
Glenn L, Martin Co. 

National Capital 


JAMES R. PATTON, JR. 
Office of Naval Research 
New Merico- 

West Texas FRANK L. KOEN, JR. 
White Sands Proving Ground 
MICHAEL J. SAMEK 
American Electro-Metal Corp. 

Niagara Frontier 
WILLIS SPRATTLING 
Bell Aircraft Corp. 


New York 


Northeastern 

New York GEORGE E. MOORE 
General Electric Co. 

Northern 
California WILLIAM J. BARR 
Detroit Controls Corp. 

Southern 
California H. S. SEIFERT 


Calif. Inst. of Tech. 
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dirs ss S. SEIFERT, California Institute of Technology, Associate Editor 


Introduction to Aeronautical Dynamics, 
by M. Rauscher, John Wiley & Sons, 
Inc., New York, N. Y., 1953, 664 pp. 
$12. 

Reviewed by H. J. Stewart 
California Institute of Technology 
Jet Propulsion Laboratory 


This book presents essentially the con- 
tent of a group of courses which Prof. 
Rauscher taught for many years in the 
Aeronautical Engineering Department of 
the Massachusetts Institute of Technology. 
With this background the material obvi- 
ously has the advantage of having been 
used many times; thus the methods of 
presentation and the choice of examples 
are well suited to the subject. 

The material covered in the book con- 
sists of various topics in solid dynamics, 
such as motion of particles and of rigid 
bodies and vibrations, and topics in fluid 
mechanics, such as one-dimensional com- 
pressible-flow theory, and two-dimensional 
perfect-fluid theory and the elements of 
viscous-fluid theory, which are required 
as a background for the discussion of the 
dynamics of aircraft. The subject ma- 
terial is carried far enough to discuss two- 
dimensional airfoil theory, the Prandtl 
monoplane-wing theory and vibrations of 
systems with more than one degree of free- 
dom; however, the treatment stops short 
of the final goal—the dynamics of airplanes 
in flight. 

Since the subject matter covered is in- 
tended for presentation to undergrad- 
uate classes, Professor Rauscher has wisely 
limited the mathematical methods to the 
most elementary level required. In the 
treatment of two-dimensional airfoil the- 
ory he introduces and uses the theory of 
complex variables and conformal trans- 
formation in order to discuss Joukowsky 
airfoils and von Kérmén-Treftz airfoils. 
He also introduces the elements of Fourier 
analysis in his treatment of ‘vibrations. 
With these two exceptions the mathemati- 
eal level is only that which should have 
been acquired by any undergraduate 
junior student. 


The Complete Book of Outer Space, J. 
Logan, Editor, Maco Magazine Corp., 
New York, N. Y., 1953, 144 pp. $0.75. 

Reviewed by A. J. ZAEHRINGER 

Thiokol Chemical Corporation, Elkton, 

Md. 


Despite the pretentious title of this 
magazine-style book, it represents a 
simple and concise picture of rocketry. 
Willy Ley presents a conventional pic- 
ture of rocket history and in a later sec- 
tion discusses the possibilities of life be- 
yond our solar system. Dr. von Braun 
has two sections dealing with the space 
station and urging standardization and 
immediate starts for an orbiting rocket. 
A very optimistic section on space medi- 
cine is discussed by Dr. Heinz Haber which 


essentially is a condensation of his recent 
book, “Man in Space.” A section on 
space suits is given by Dr. Menzel of the 
Harvard Observatory in which he visuali- 
zes the space suit as a sort of manned 
automobile tire. Robert Haviland of the 
GE Project Hermes gives an interesting 
picture of multistage rocket development 
in the Bumper program. One of the most 
interesting and longest sections is the his- 
tory of rocket motors by James Wyld. 
Developments from the Paulet engine to 
the author’s experiences with early RMI 
rocket motors are presented. Other sec- 
tions of the book are: Legal Aspects of 
Space Travel, by Oscar Schacter of the 
UN; Exploitation of the Moon, by Hugo 
Gernsbach; Interstellar Flight, by Dr. 
Leslie Shephard, a director of the British 
Interplanetary Society; and two sections 
by the editor, The Spaceship in Science 
Fiction, and the Flying Saucer Myth. 
Also presented are ‘“‘charts’”’ of Moon and 
Mars flights, space-ship timetables, tem- 
peratures of the planets, and a space- 
travel dictionary. 

The publishers are to be commended for 
presenting an interesting book of semi- 
technical caliber at low cost. Although 
it is paper bound, the printing and illus- 
trations (almost 200) are excellent. How- 
ever, many of the illustrations of science- 
fiction type could have been weeded out 
together with many irrelevant illustra- 
tions, viz., a drawing of the F-99 Bomare 
in the section on high-altitude research. 
Nevertheless, the general technical ac- 
curacy of the material is excellent, and the 
authors reflect a highly optimistic outlook 
toward manned space flight. Although 
somewhat scrambled in method of pres- 
entation, this book can be expected to 
give many persons a fairly accurate pic- 
ture of space flight. 


Flying Saucers, by D. H. Menzel, Har- 
vard University Press, Cambridge, 
Mass., 1953, 319 pp. $4.75 

Reviewed by J. E. 
California Institute of Technology 
Jet Propulsion Laboratory 


Based on his own experience with opti- 
cal displays in the atmosphere, Professor 
Menzel presents an interesting and plau- 
sible explanation of the flying-saucer phe- 
nomena. History is apparently studded 
with recordings of public consternation at 
strange visitors from the skies, and this 
book includes a fascinating summary of 
these historical events. 

Flying Saucers is more than a historical 
summary. An adequate description of the 
physical concepts which explain optical 
displays in the skies is also given in the 
book. In addition, the author outlines 


some good fundamental rules for evaluat- 
ing critically observations of any kind. 
Professor Menzel unfortunately makes the 
same error he criticizes in other treat- 
ments of the subject. 


Throughout the 


book he forces the conclusion that the 
flying-saucer phenomena arise from natu- 
ral causes and are certainly not the re- 
sult of interplanetary vehicles. It is 
understandable that the author is more 
willing to believe in optical tricks than in 
space ships, but his analysis of flying 
saucers presents no new facts to ‘“‘prove” 
either conclusion; rather it falls into 
the category of hypothesis, as do other 
writings on this subject. Otherwise it is 
an excellent book, one toward which ergi- 
neers and scientists will feel sympathetic. 

The illusions of the atmosphere as <e- 
scribed in this book are enough to m:ke 
one cautious in accepting any explanation 
of flying saucers on the basis of the ob- 
servations that have been reported to 
date. 

The examples which are used to illus- 
trate the theory are unusual for college 
texts and are especially suitable for modern 
students interested in aeronautics or in jet 
propulsion. For example, the author con- 
siders the motion of a rocket in a plane- 
tary gravitational field as an example 
in particle dynamics. He also considers 
the theory of a ramjet, in a low-speed zero- 
drag approximation, as an example in one- 
dimensional fluid flow. Again, in the 
section on dynamics of rigid bodies, he has 
an introduction of the effects of rotating 
propellers on the motion of an airplane ina 
spin. 


Mechanical Vibrations, by W. T. Thom- 
son, Prentice-Hall, Inec., New York, 
N. Y., 1953, 252 pp. $6. 

Reviewed by R. H. BopEN 
North American Aviation 


This textbook covers in detail second- 
order linear systems. Many examples are 
presented of such systems and their elec- 
trical analogs, both alone and in com- 
bination as many-degree-of-freedom sys- 
tems. Since it is presented clearly and in 
great detail, this book is a worth-while 
reference for the practicing engineer, al- 
though the material is not new. 

The objective of presenting the funda- 
mentals of vibration theory and providing 
a general background for advanced study 
in the field is achieved. The author does 
not stray from a presentation of the funda- 
mentals of second-order linear systems. 
Some discussion of first-, third-, and 
higher-order systems, some digression into 
nonlinear types, to define more clearly the 
limits of application of the fundamental 
information and also some further <is- 

“cussion of stability, which is mentioned 
only briefly, would enhance the value of 
the text without overstepping the bounds 
of an introductory course in vibration 
analysis. Another worth-while addition 
would be a more extensive discussion of 
transient phenomena in simple systems 
and their relationship to steady-state per- 
formance. In general, the book appears 
limited but quite satisfactory as the basis 
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of an introductory course in vibration 
‘analysis, particularly for the mechanical 
engineer. 


Information Theory, by S. Goldman, 
Prentice-Hall, Inc., New York, N. Y., 
1953, 385 pp. $9. 


Reviewed by R. M. Srewarr anv E. 
RECHTIN 
California Institute of Technology 
Jet Propulsion Laboratory 


This book is one of the first published 
attempts to rewrite the work of Shannon 
and Wiener in another author’s language 
for textbook purposes. Because of this 
unique status, if for no other reason, the 
book will probably be very popular and 
certainly should be examined by all 
specialists in the field. These reviewers, 
however, were disappointed in certain 
aspects of the book. In particular, the 
rewriting of Shannon’s work did not seem 
as successful as the original work, and that 
of Wiener is considerably less satisfactory 
than that of Bode and Shannon (Proc. 
IR#, vol. 38, 424, April 1950). The 
author expresses unnecessary pessimism 
when he concludes that ‘‘the actual de- 
sign of an optimum filter in any practical 
ease is invariably laborious.” 

Goldman’s language differs from Shan- 
non’s in several ways. Various unusual 
words (quadiva, liniva, quadratic content, 
linear content, binits) are used for old 
concepts in lieu of well-established words, 
and some old words (random noise, er- 
godic) are given new meanings. Goldman’s 
definition of ‘information received” is 
different from Shannon’s in a way which 
may lead to considerable confusion. In 
the last chapter, the author delves into 
speculative theory of his own, usually a 
source of misunderstanding in a book writ- 
ten for instructional purposes. 

The author is to be commended, how- 
ever, for his effort to clear the air in a 
rapidly growing field; but these reviewers 
feel that the effort was not particularly) 
successful. 


Meteorological Instruments, third edition, 
revised, by W. E. K. Middleton and 
A. F. Spilhaus, University of Toronto 
Press, Toronto, Canada, 1953, 286 pp. 
$11.50. 


Reviewed by B. HELFAND 
North American Instruments, Ine. 
Altadena, Calif. 


Instrumentation for the meteorologist 
includes an unusual assortment of de- 
vices and techniques designed to enable 
him to specify the weather. The devices 
range in complexity from a simple mercury 
and glass thermometer to a search radar, 
and it is with these instruments that he 
must make measurements for use in 
weather analysis, forecasting, and climato- 
logical work. 

The book, now in its third edition, de- 
scribes the state of the art as practiced on 
this continent and in most of the English- 
speaking countries and is based on Mr. 
Middleton’s lecture course given for a 
number of years to graduate students in 
meteorology at the University of Toronto. 
Dr. Spilhaus, Dean of the Institute of 
Technology, University of Minnesota, 
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appears as co-author for the first time in 
this edition. He is well known in this 
country for his contributions in the field 
of meteorological instruments. 

The book is written in an easy-to-read 
conversational style and is well illustrated. 
The content is largely descriptive, typical 
instruments and their usage being covered 
in considerable detail. Analytical treat- 
ment is provided where necessary to give 
an understanding of static and dynamic 
performance. Many references are pro- 
vided to lead the reader to original sources 
of information and commercial sources of 
particular instrumentation. War-time 
developments in the application of radio 
and radar to upper-air-measurement prob- 
lems constitute much of the new material 
added in this edition. 

The introductory chapter deals with the 
special properties and requirements for 
meteorological instruments. Five sub- 
sequent chapters cover the measurement 
of atmospheric pressure, temperature, 
humidity, precipitation, and wind and 
ground stations. One chapter is devoted 
to sunshine recorders, one to instruments 
for investigating clouds, and two to the 
investigation of the upper air. The final 
chapter deals with special assemblies and 
systems including portable, mobile, and 
automatic weather stations. 

This book is the only comprehensive 
source of information on meteorological 
instruments written in English; it is 
highly recommended as a text for me- 


teorology students at both the under- 
graduate and graduate levels and also as a 
reference for practicing meteorologists 
and those in other fields with an interest 
in temperature, pressure, and flow meas- 
urements. 


Principles of Automatic Controls, by |’. F. 
Nixon, Prentice-Hall, Inc., New York, 
N. Y., 1953, 409 pp. $9.35. 

Reviewed by R. J. Parks 
California Institute of Technology 
Jet Propulsion Laboratory 

This book was written as an outgrowth 
of a company training course given by the 
author at the Glenn L. Martin Compe ny. 

In the author’s words, it is intended «3 

reference, self-study, or undergradu: te- 

level text. A very minimum of beck- 
ground training is assumed. As state: in 
the preface, a first-year course in physivs is 
assumed and a course in electrical circ tits 
would be helpful but not essential. 
Consistent with this premise, the bok 
starts out in a fairly elementary man ver, 
discussing linear differential equations, the 

Laplace transform, block-diagram nota- 

tion, and the theory and application of the 

Nyquist criteria. These topics are lis- 

cussed from the applied standpoint, wiih a 

minimum of mathematical rigor. A con- 

siderable number of simple examples are 
discussed and, in addition, problems are 
included at the end of each chapter. 

The book then discusses other methods 
of servo analysis such as the inverse open- 
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loop, transfer-function method, the log- 
decibel method, and the root-locus tech- 
nique. These subjets are well discussed 
from the standpoint of one unfamiliar 
with the field. These techniques and 
that of the Nyquist criteria are discussed 
completely enough to point out the various 
pitfalls which sometimes occur if these pro- 
cedures are applied blindly without a com- 
plete understanding of their development. 

Several chapters are devoted to giving 
the reader a feeling for the problems en- 
countered in the design of feedback sys- 
tem: together with typical solutions to 
these problems. Sections are included on 
nun-erical integration and the techniques 
of automatic computing, both analog 
and digital, as well as their application to 
servo design. 

The important subject of noise and its 
effe.ts on closed-loop systems is discussed 
onl, briefly and qualitatively. Six appen- 
dixes include information for reference 
incliding the Runge-Kutta integrator 
tecl nique and a rather complete table of 
Laplace transforms. 

I: is the opinion of the reviewer that the 
book is well written and accomplishes its 
stated goals. It is not a high-level tech- 
nical book, and many readers will want to 
skin. over the first few chapters, but it 
would make excellent introductory reading 
for persons whose experience in the field of 
controls is limited. 


Weight Strength Analysis of Aircraft 
Structures, by F. R. Shanley, McGraw- 
Hill Book Co., Inc., New York, N. Y., 
1952, 394 pp. $8.50. 

Reviewed by M. L. WrILu1ams 
California Institute of Technology 
Professor Shanley’s book represents the 
first attempt to gather together in a com- 


| 
prehensive treatment a rational method for 


evaluating the interaction of weight and 
strength in aircraft structures. Analytical 
expressions are developed for the weights 
as a function of the design parameters 
and the strength level, or structural index. 
As opposed to solely statistical-trend 
data, the analytical formulation permits 
a greater understanding of the weight 
breakdown and contributes to more effi- 
cient design. In particular it is shown 
that the minimum structural weight can 
be separated into two parts, one depending 
on the loads, material properties, and load 
transmission path, and the other upon the 
size and geometrical arrangement. 

The difficulties involved in the direct 
application of the material to aircraft 
design stem from the proper evaluation or 
prediction of certain constants in the 
analytical expressions which in many cases 
require additional statistical data. 

The book is separated into three parts: 
(I) Principles of Optimum Structural De- 
sign; (II) Structural Weight Equations, 
covering shell structures and high-aspect 
ratio, thick- and thin-walled wings with 
and without sweepback; and (III) Ma- 
terial Properties and Behaviors, a some- 
what different, although related, subject 
dealing with elastic and inelastic effects 
particularly in column behavior. 

I: is believed that one of the major 
coniributions of the volume is to present 
the problem and recommend a formulation 
upon which to base future refinements and 
adv neces, 
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Jet Propulsion Engines 


Jet Propulsion for Converti-planes, by 
David C. Prince, American Helicopter, vol. 
32, Oct. 1953, pp. 6-9, 13-16. 

Supersonic Ram-jet Performance, by 
J. Lukasiewicz, Aircraft Engineering, vol. 
25, Oct. 1953, pp. 298-306. 

Small Turbines—The Heart of Modern 
Aircraft Accessory Power Systems, by 
L. W. Bieuer and R. L. McManus, Soc. 
Automotive Engrs., Prepr. no. 159, Sept. 
1953, 9 pp., 8 figs. 

The Performance of Some Typical 
Turbo-jet Engine Exhaust Systems with 
Particular Reference to the Effects of 
Swirl, by P. F. Ashwood and P. J. 
Fletcher, Gt. Brit. Aero. Res. Counc. 
Curr. Pap. no. 130 (formerly ARC Tech. 
Rep. No. 14736, National Gas Turbine 
Est. Rep. no. R. 112), 1953, 13 pp. 

Compressor Testing at Ansty. (Arm- 
strong Siddeley Motors, Ltd.), Aero- 
plane, vol. 85, Nov. 13, 1953, pp. 658-663. 

Effect of a Turbojet Engine on the 
Dynamic Stability of an Aircraft, by Hans 
C. Vetter, J. Aero. Sct., vol. 20, Nov. 
1953, pp. 797-798. 

Turboprops, What Europe Has Learned, 
by James Hay Stevens, American Avia- 
tion, vol. 17, Nov. 13, 1953, pp. 28-29. 

A Method of Calculating the Perform- 
ance of Turbojets and Turboprops, Useful 
Expressions for Research in Propulsion 
and in the Mechanics of Flight, by 
Gaspar Santangelo, L.’ Aerotecnica, vol. 33, 
June 1953, pp. 204-219, 230 (in Italian). 

On Some Possible Development Trends 
in Aircraft Engines, by Lorenzo Poggi, 
L’ Aerotecnica, vol. 33, June 1953, pp. 
199-204 (in Italian). 

Control for Hypothetical Ramjet Must 
Manage 300,000 hp, Fit in 4 cu. Ft., by 
J. C. Wise, SAE Journal, vol. 61, Nov. 
1953, pp. 85-86. 

Some Stall and Surge Phenomena in 
Axial-Flow Compressors, by Merle C. 
Huppert and William A. Benser, J. Aero. 
Sci., vol. 20, Dec. 1953, pp. 835-845. 

Fairchild Monocoque Turbojet, by 
Randolph Hawthorne, Aviation Age, vol. 
20, Dec. 1953, pp. 30-31. 

Design for Turbojet Simplicity, by 
A. L. Ervin, Aviation Age, vol. 20, Dec. 
1953, pp. 28-29. 

A Slipstream View of Aircaft Propulsion, 
1903-1953, by Frank W. Caldwell, Aero. 
Engng. Rev., vol. 12, Dec. 1953, pp. 108- 
117, 131-132. 

The Compound Engine, by C. Casci, 
L’ Aerotecnica, vol. 33, Aug. 1953, pp. 
263-269 (in Italian). 

Pressure Ratios for Aero-Engines, by 
J. M. Stephenson, Aircraft Engng., vol. 25, 
Dec. 1953, p. 371. 

A Fuel-Air Turbine Starter, Aeroplane, 
vol. 85, Dec. 11, 1953, pp. 796-798. 

Jet Ignition Developments Moving 
Fast, by Joseph S. Murphy, American 
Aviation, vol. 17, Jan. 4, 1954, pp. 44-45. 


The Case for the By-Pass Jet Engine, 
by James Hay Stevens, American Avia- 
tion, vol. 17, Jan. 4, 1953, pp. 28-30. 

Use of Electric Analogs for Calculation 
of Temperature Distribution of Cooled 
Turbine Blades, by Herman H. Eller- 
brock, Jr., Eugene F. Schum, and Alfred 
J. Nachigall, NACA TN 3060, Dec. 1953, 
116 pp. 

Some Factors Pertaining to the Use of 
Air Breathing Propulsion for the Accelera- 
tion of High Altitude Sounding and Other 
Long Range Ballistic Missiles, by J. W. 
Luecht, Soc. Automotive Engr. Prep., No. 
145, Sept. 1953, 5 pp., 9 figs. 


Rocket Propulsion Engines 


The Effect of Variation of Propellant 
Density on Rocket Performance, by 
Jack Lorell and Albert Hibbs, Calif. Inst. 
Tech. Jet Prop. Lab. Rep. no. 20-66, 
July 1952, 33 pp. 

Some Limiting Factors of Chemical 
Rocket Motors, by W. N. Neat, J. Brit. 
Interplan. Soc., vol. 12, Nov. 1953, pp. 
249-274. 

Propulsion Efficiency of a Rocket, by 
Enio Mattioli, L’Aerotecnica, vol. 33, 
Aug. 1953, pp. 270-274 (in Italian). 

Flow in Injection Orifices, Application 
to Rocket Engines, by R. Kling and R. 
Leboenf, Recherche Aéron., no. 35, Sept.- 
Oct. 1953, pp. 35-41 (in French). 


Heat Transfer and Fluid 
Flow 


The Coriolis Effect, by Paul E. Wylie, 
J. Roy. Aeron. Soc., vol. 57, Oct. 1953, pp. 
655-658. 

The Fluid Mechanics of Detonation 
Waves, by L. G. Dawson, J. Roy. Aeron. 
Soc., vol. 57, Oct. 1953, pp. 618-626. 

An Analysis of the Air Flow Through the 
Nozzle Blades of a Single Stage Turbine, 
by I. H. Johnston, Gi. Brit. Aero. Res. 
Counc. Curr. Pap. no. 131 (formerly ARC 


Tech. Rep. no. 14300, National Gas 
Turbine Est. Memo. no. M. 108), 1953, 
8 pp., 10 figs. 


The Design and Testing of Supersonic 
Nozzles, by R. Harrop, P. I. F. Bright, 
J. Salmon, and M. J. Caiger, Gt. Brit. 
Aero. Res. Counc. Reps. and Memo. 
no. 2712 (formerly ARC Tech. Reps. nos. 
12114, 13383; Royal Aircraft Est. Reps. 
Aero. 2293 2293a), 1953, 40 pp. 

Turbulence in Supersonic Flow, by 
Leslie S. G. Kovasznay, J. Aero. Sci., 
vol. 20, Oct. 1953, pp. 657-674, 682. 

Recent Developments in Convective 
Heat Transfer with Special Reference to 
High-Temperature Combustion Cham- 
bers, by Martin Summerfield, Heat 
Transfer, A Symposium, Univ. of Mich., 
1952, pp. 151-171, 25 refs. 

Evaporation and Spreading of Isooctane 
Sprays in High-Velocity Air Streams, by 


Donald W. Bahr, NACA RM E53114, 
Nov. 1953, 35 pp. 

Air Cooling Methods for Gas Turhine 
Combustion Systems, by F. J. Bagiey, 
Gt. Brit. Aero. Res. Counc. Curr. Pap. no. 
133 (formerly ARC Tech. Rep. 14:48, 
National Gas Turbine Est. Rep. R. 151), 
1953, 34 pp. 

Cavitation Tests on Hydrofoils in Cas- 
cade, by Fukusaburo Numachi, 7'r:ns. 
ASME, vol. 75, Oct. 1953, pp. 1257-1269, 

On the Propagation and Structure of the 
Blast Wave, Part I, by Akira Sakurai, 
J. Phys. Soc. of Japan, vol. 8, Sept.-Oct. 
1953, pp. 662-669. 

Calculation of Transpiration Covled 
Gas-Turbine Blades, by J. N. B. Livin- 
good and E. R. G. Eckert, Trans. AS\/E, 
vol. 75, Oct. 1953, pp. 1271-1278. 

On the Evaluation of the Accuracy of the 
Coefficient of Discharge in the Basic 
Flow-Measurement Equation, by A. L. 
Jorissen, Trans. ASME, vol. 75, Oct. 
1953, pp. 1323-1326. 

A Rapid Method for the Determination 
of the Thermal Conductivity of Organic 
Liquids for Use in Heat Transfer Calcu- 
lations, by S. Baxter, H. A. Vodden, and 
Davies, J. Appl. Chemistry, vol. 3, part 10, 
Oct. 1953, pp. 477-480. 

Local Heat-Transfer Coefficients on 
Surface of an Elliptical Cylinder, Axis 
Ratio 1:3, in a High-Speed Air Stream, 
by R. M. Drake, Jr., R. A. Seban, and 
D. L. Doughty, and 8S. Levy, T'rans. 
ASME, vol. 75, Oct. 1953, pp. 1291-1302. 

New Aspects of Natural-Convection, 
Preliminary Study of Effect of Frictional 
Heating, by Simon Ostrach, Trans. 
ASME, vol. 75, Oct. 1953, pp. 1287- 
1290. 

Flow Through Nozzles and Related 
Problems of Cylindrical and Spherical 
Waves, by Yu Why Chen, Comm. Pure 
Appl. Math., vol. 6, May 1953, pp. 179- 
229. 


On the Solution of the Differential 
Equation Occurring in the Problem of 
Heat Convection in Laminar Flow Through 
a Tube, by Milton Abramowitz, J. Math. 
and Physics, vol. 32, July—Oct. 1953, pp. 
184-187. 

Note on the Prandtl Number for Dis- 
sociated Air, by C. Frederick Hansen, J. 
Aero. Sci., vol. 20, Nov. 1953, pp. 789-790. 

Measurement of the Profile Drag of 
Compressor and Turbine Cascades and 
the Effect of Wakes in Exciting Vibra- 
tion, by J. M. Stephenson, J. Roy. Aeron. 
Soc., vol. 57, Nov. 1953, pp. 722-725. 

Evidences of an Inherent Error in 
Measurement of Total-Heat Pressure at 
Supersonic Speeds, by James S. Murphy, 
Aero. Engng. Rev., vol. 12, Nov. 1953, pp. 

7-51. 

Analogy Between Mass and Heat 
Transfer With Turbulent Flow, Edmund 
E. Callaghan, NACA TN 3045, Oct. 1958, 
19 pp. 


Eprtor’s Note: This collection of references is not intended to be comprehensive, but is rather a selection of the most significant and 


stimulating papers which have come to the attention of the contributors. 
literature is unavailable because of security restrictions. 


The readers will understand that a considerable body of 
We invite contributions to this department of references which have not come 


to our attention, as well as comment on how the department may better serve its function of providing leads to the jet propulsion app!ica- 


tions of many diverse fields of knowledge. 


130 


JET PROPULSION 


Sk 
IV 
by 
Bi 
14 
Re 
21 
4 
Mi 
Ka 
19% 
] 
Fil 
| No 
ber 
pp. 
\ 
Pla 
In 
263: 
Kn 
Tur 
Mo 
Fl 
in X¢ 
Equi’ 
jet 
Ord. 
App 
Color 
Flame 
Hydro 
M. F. 
TN 30 
Pras 
Tempe 
Millar 
Soc. / 
Nov. 1 
Dev 
tive Di 
Loefile 
P. repr. 
-§ Stud 
tion in 
Robbin 
Chemisi 
Re No 
Burn’ 
10. Fi, 


LON 


The Mechanics of Film Cooling, by 
Eldon L. Knuth, Calif. Inst. Tech. Jet 
Prop. Lab. Memo. 20-85, Sept. 1953, 108 


preliminary Investigation of a New 
Type of Superonic Inlet, by Antonio 
Ferri and Louis M. Nucci, NACA Rep. 
1104 (formerly T'N 2286), 1952, 19 pp. 

The Measurement of Heat Transfer and 
Skin Friction at Supersonic Speeds, 
Pt. III. Measurements of Overall Heat 
Transfer and of the Associated Boundary 
Layers on a Flat Plate at M: = 2.43, by 
R. J. Monaghan and J. R. Cooke, Gt. 
Brit. Aero. Res. Counc. Curr. Pap. no. 
139 (formerly ARC Tech. Rep. 15348, 
Roval Atreraft Est. Tech. Note No. Aero. 
2120), 1953, 38 pp. 

The Measurement of Heat Transfer and 
Skin Friction at Supersonic Speeds, Part 
IV. Tests on a Flat Plate at M = 2.82, 
by 8. J. Monaghan and J. R. Cooke, Gt. 
Brit. Aero. Res. Counc. Curr. Pap. no. 
140 (formerly ARC Tech. Rep. 15400, 
Ro, al Aircraft Est. Tech. Note No. Aero. 
2171), 1953, 27 pp. 

Analogy Among Heat, Mass, and 
Momentum Transfer, by Richard F. 
Kaser, Ind. Engng. Chem., vol. 45, Dec. 
1953, pp. 2634-2636. 

Heat Transfer in Forced Convection 
Film Boiling, by LeRoy A. re 4 
Norman R. LeRoy, and James A. Ro 
bers, Ind. Engng. Chem., vol. 45, Dec. 1953, 
pp. 2639-2646. 

Velocity Distribution Between Parallel 
Plates, by W. G. Schlinger and B. H. Sage. 
Ind. Engng. Chem., vol. 45, Dec. 1953, pp. 
2636-2639. 

Sweat Cooling, a Review of Present 
Knowledge and Its Application to the Gas 
Turbine, by P. Grootenhuis, and N. P. W. 
Moore, Iron and Steel Institute. Sym- 
posium on High Temperature Materials 
for Gas Turbine, 1952, pp. 281-289, 47 ref. 

Surface Temperature and Pressure Dis- 
tributions on a Circular Cylinder in Super- 
sonic Cross Flow, by L. W.. Walter and 
H. H. Lange, Naval Ord. Lab. NAVORD 
Rep. 2854 (Aeroballistic Res. Rep. no. 
180), June 1953, 10 pp. 

Flow Properties of Strong Shock Waves 
in Xenon Gas as Determined for Thermal 
Equilibrium Conditions, by Alexander P. 
Sabol, NACA TN 3091, Dec. 1953, 29 pp. 

On the Stability of Isotropic Turbulence, 
by Robert Betchov, Maryland University. 
Institute for Fluid Dynamics and Applied 
Mathematics, Tech. Rep. BT-15, Sept. 
1953, 26 pp. 

Jets, by M. Z. Krzywoblocki, Naval 
Ord. Test Sta. Tech. Memo. no. 1576, Sept. 
1953, 83 pp. 


Combustion 

Application of an Electro-Optical Two- 
Color Pyrometer to Measurement of 
Flame Temperature for Liquid Oxygen- 
Hydrocarbon Propellant Combination, by 
M. F. Heidmann and R. J. Priem, NACA 
TN 3033, Oct. 1953, 39 pp. 

Practical Application of Engine Flame 
Temperature Measurements, by G. H. 
Millar, O. A. Uyehara, and P. S. Myers, 
Soc. Automotive Engrs. Prepr. no. 196, 
Nov. 1953, 9 pp. 24 figs. 

Development of an Improved Automo- 
tive Diesel Combustion System, by Bruno 
Loeffler, Society of Automotive Engrs. 
Prepr. no. 188, Nov. 1953, 16 pp. 47 figs. 

Studies on Carbon Formation and Reac- 
tion in Acetylene, by Jerry D. Frazee and 
Robbin C. Anderson, Texas Univ. Dept. 
Chemistry Combustion Kinetics Proj. Tech. 
Rep. no. 8, Aug. 1953, 17 pp. 

Burning Velocity Determinations. Part 
10. Flame Propagation Along Horizontal 
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Though miniature in size and light in 
weight, Bourns Gage Pressure Potentiom- 
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ball bearing supported linkage system trans- 
___ mits the movement of the Bourdon tube to the 
sliding contact. Instruments are not damaged 
by shock, vibration or nominal over-pressure. 


Bourns designs and manufactures other 
precision potentiometers for measuring position, 
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The Chance-Vought Cutlass was 
an advanced design in 1947. It is 
a leader today. The same is true 
of Airborne’s LINEATOR, which 
was developed for the Cutlass 
flight control system. Today, seven 
years later, there is not another 
“tee” type linear actuator like it. 

The same basic model is used 
in the latest of the Cutlass series, 
and in the McDonnell Banshee. 
Modifications of the LINEATOR 
are specified equipment in the 
McDonnell Demon; its Air Force 
companion, F-101; and the Mar- 
tin P5M patrol bomber. 

Conforming to MIL-A-8064 
(USAF), the LINEATOR is most 
adaptable where light weight and 
short length, for a given stroke, 
are desirable features. A ball bear- 
ing jack screw enables it to handle 
1500 Ib. maximum operating load 
in either tension or compression. 

Airborne has set the pace in the 
actuator field with advanced de- 
signs like the LINEAToR. As air- 
craft configurations change and 
speeds increase, count on Air- 
borne for more of the same. For 
information on the LINEATOR and 
other actuators, see our literature 
in the 1.A.S. Catalog. 
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Tubes, by M. F. Hoare and J. W. Lin- 
nett, Trans. Faraday Soc., vol. 49, pt. 9, 
Sept. 1953, pp. 1038-1049. 

Investigation of Flame Temperatures in 
a Single-Cylinder Spark-Ignition Engine, 
by J. H. Potter, and E. B. Dillaway, 
Trans. ASME, vol. 75, Oct. 1953, pp. 
1311-1321. 

Combustion of a Low-Volatility Fuelina 
Turbojet Combustion Chamber-Effects of 
Fuel Vaporization, by V. V. Holmes, A. J. 
Pahnke, O. A. Uyeharo, and P. 8. Myers, 
Trans. ASME, vol. 75, Oct. 1953, pp. 
1303-1310. 

Temperature and Gas-Analysis Sur- 
veys in the Combustion Zone of a Gas- 
Fired Gas-Turbine Combustor, by K. L. 
Rieke, Trans. ASME, vol. 75, Oct. 1953, 
pp. 1233-1239. 

The Slow Combustion of Cyclopropane. 
II. Analytical Results and Mechanism, 
by A. C. McEwan and C. F. H. Tipper, 
Proc. Roy. Soc., vol. 220A, Nov. 10, 1953, 
pp. 266-277. 

Effect of Diffusion Processes and Tem- 
perature on Smoking Tendencies of 
Laminar Diffusion Flames, by Rose L. 
Schalla, NACA RM E52122, Dec. 1953, 
23 pp. 
Formation and Combustion of Smoke In 
Bunsen Flames, by Thomas P. Clark, 
Ind. Engng. Chem., vol. 45, Dec. 1953, pp. 
2785-2789. 

Thermodynamic Charts of Combustion 
Products, by Ernesto Macioce, L’ Aero- 
tecnica, vol. 33, Aug. 1953, pp. 288-290 
(in Italian). 

Method of Electronic Statistical Meas- 
urement For the Study of Combustion 
Phenomena In Reciprocating Engines, by 
R. Vichnievsky, and Weissmann, Re- 
cherche Aéron., no. 35, Sept. —Oct. 1953, pp. 
423-435 (in French), 36 references. 

A Note on Organ-Pipe Resonance in 
Ducted Burners, by J. C. Truman and 
N. J. Lipstein, J. Aero. Sci., vol. 20, Dec. 
1953, pp. 846-847. 

Effects of Additives on Pressure Limits 
of Flame Propagation of Propane-Air 
Mixtures, by Frank E. Belles and Doro- 
they M. Simon, NACA RM E53129, Dec. 
1953, 36 pp. 


Fuels, Propellants, and 
Materials 


Stability of Hydrazine in Liquid Am- 
monia, by John A. Krynitsky and Homer 
W. Carhart, Naval Res. Lab. Rep. no. 4219, 
Sept. 1953, 16 pp. 

Hydrazine, by Lawrence P. Lessing, 
Sci. Amer., vol. 189, July 1953, pp. 30-33. 

Hydrogen Peroxide, Part One, by W. C. 
Schumb, C. N. Satterfield, and R. L. 
Wentworth, Mass. Inst. Tech. Depts. 
Chemistry and Chemical Engng., Div. 
Indus. Cooperation Rep. no. 42, Sept. 1953, 
247 pp. 

Hydrocarbon and Nonhydrocarbon De- 
rivatives of Cyclopropane, by Vernon A. 
Slabey, Paul H. Wise, and Louis C. Gib- 
bons, "NACA Rep. 1112, 1953, 18 pp. 

Formation, Stability and Crystal Struc- 
ture of the Solid Aluminum Suboxides, 
Al,O and AlO, by Michael Hoch and Her- 
rick L. Johnston, Ohio State Univ. Res. 
Found. Tech. Rep. 280-10, Aug. 1953, 6 


p. 

° Effective Thermal Conductivities of 
Magnesium Oxide, Stainless Steel, and 
Uranium Oxide Powders in Various Gases, 
by C. S. Eian and R. G. Deissler, NACA 
RM E53G03, Oct. 1953, 18 pp. 

Hydrogen Peroxide. Part 4, by W. C. 
Schumb, C. N. Satterfield and R. L. 
Wentworth, Mass. Inst. Tech. Dept., 
Chem. and Chemical Engng., Div. Indus. 
Cooperation, Rep. no. 45, Nov. 1953, 194 


pp. 


Developing Aviation Fuels and Lubri- 
cants, by C. B. Davies, J. Roy. Aeron, 
Soc., vol. 57, Nov. 1953, pp. 700-708. 

Velocity of Sound in Liquid Propellants, 
by Robert S. Wicks, Irving L. Odgers, 
and Rollow 8S. Pickford, Jr., Calif. /nst, 
Tech., Jet Prop. Lab., Memo. 20 84, 
July 1953, 11 pp. 

Properties, Production and Uses of 
Hydrazine, by James E. Troyan, /nd. 
Engng. Chem., vol. 45, Dec. 1953, pp. 
2608-2612. 

Volumetric Behavior of Nitric Acid, by 
H. H. Beamer, W. H. Corcoran, and B. H. 
Sage, Ind. Engng. Chem., vol. 45, ec. 
1953, pp. 2699-2704. 

Hydrogen Peroxide. Part 2, by W C. 
Schumb, C. N. Satterfield, and R. L, 
Wentworth, Mass. Inst. Tech., Di pts. 
Chem. and Chemical Engng., Rep. 43, | dec. 
1953, 232 pp. 

The Anelasticity of Molybdenum. F nal 
Report Covering the Period May 1, 1°50 
July 14, 1953 and 16th Quarterly Re: ort 
Covering the Period March 1-July 14, 
1953, Battelle Mem. Inst., 43 pp. 

The Improvement of the Ductility of 
Molybdenum Weldments, Summary Re- 
port Covering the Period Aug. 1, 1951 
July 14, 1953 and 16th Quarterly Report 
Covering the Period March 1—July 14, 
1953, Battelle Mem. Inst., 13 pp. 


Physical-Chemical Topics 

Minimum Spark-Ignition Energies 0: 12 
Pure Fuels at Atmospheric and Reduced 
Pressure, by Allen J. Metzler, NACA ?}.M 
E53H31, Oct. 1953, 28 pp. 

The Oxidation of Carbon Monoxide in 
the Presence of Ozone, by David Garvin, 
Princeton Univ. James Forrestal 
Cent., Chemical Kinetics Proj., Tech. lep. 
no. 5, Sept. 1953, 22 pp. 


Over 85% of the torque wrenches 
used in industry are 


TORQUE WRENCHES 


Read by or Feel. 


@ Permanently 

@ Practically Indestructible 
@ Faster—Easier to use 

@ Automatic Release 


@ All Capacities 


in inch ounces... inch 
pounds . . . foot pounds 
(All Sizes from 0-6000 

ft. Ibs.) 


manufacturer, 
design ond 
production man 
should hove 
this valuable 
dota. Sent upon 
request. 
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The Reaction Between Nitric Oxide and 
Ammonia On Supported Catalysts, by 
John R. Kelso, Abderdeen Proving Ground, 
Ballistic Res. Lab. Memo. Rep. no. 699, 
July 1953, 15 pp. 

Thermodynamic Properties of Boron 
and Aluminum Compounds. Progress Rep. 
no. 2, Pennsylvania State Univ., Dept. 
Chemistry, Jan. 1-Sept. 30, 1953, 50 pp. 

The Viscosity of Vapor Mixtures of 
Hydrogen Peroxide and Water, by C. N. 
Satierfield, R. L. Wentworth, and 8S. T. 
Demetriades, Mass. Inst. Tech., Dept. 
Chemical Eng., Div. Indus. Cooperation, 
Rep. no. 39, Aug. 1953, 36 pp. 

Kinetics of the Initial Reaction Between 
an Aldehyde and Hydrogen Peroxide In 
Aqueous Solution, by C. N. Satterfield, 
and L. C. Case, Mass. Inst. Tech., Dept. 
Chem. Eng., Div. Indus. Cooperation, Rep. 
no. 40, Aug. 1953, 19 pp. 

A Mollier Chart for Moisture-Saturated 
Air, H. Reichert, Aircraft Engng., vol. 25, 
Oct. 1953, pp. 321-323. 

Isotherms and Thermodynamic Func- 
tions of Ethane at Temperatures between 
0°C and 150°C and Pressures up to 200 
atni., by A. Michaels, W. vanStraaten, 
an’ J. Dawson, Maryland Univ. High 
Pr ssure Lab., Tech. Rep. no. 1, Oct. 1953, 
1 pp., 11 tab. 

The Thermal Decomposition of Organic 
Nitrates. III. The Effect of Additives 
on the Thermal Decomposition of Ethyl 
Nitrate, by Joseph B. Levy, Naval Ord. 
Lab., NAVORD Rep. 2897, Aug. 1953, 20 

The Solubility of Oxygen in Nitric Acid 
Mixtures, by Glenn D. Robertson, Jr., 
David M. Mason, and William H. Cor- 
coran, Calif. Inst. Tech. Jet Prop. Lab. 
Prog. Rep. no. 20-176, Aug. 1953, 8 pp. 

The Heats of Mixing of Some Aromatic 
Hydrocarbons and Phtholates, by P. 
Meares, Trans. Faraday Soc., vel. 49, 
Oct. 1953, pp. 1133-1139. 

The Oxidation of Hydrocarbons, II. 
Examination of the Products of the Vapour 
Phase Oxidation of n-Heptane and of 
Methylcyclohexane, by F. H. Garner, R. 
Long, and R. G. Temple, Trans. Fara- 
day Soc., vol. 49, Oct. 1953, pp. 1193-1197. 

Energy Transfer in Hydrocarbon Solu- 
tions, by E. J. Bowen and B. Brockle- 
hurst, T’rans. Faraday Soc., vol. 49, Oct. 
1953, pp. 1131-1133. 

Ionization in Solutions of Nitrogen Di- 
oxide in Nitric Acid from Optical Ab- 
sorbance Measurements, by Scott Lynn, 
David M. Mason, and W. H. Corcoran, 
Calif. Inst. Tech., Jet Prop. Lab. Prog. 
Rep. 20-175, Sept. 1953, 11 pp. 

An Approximate Theory of Inelastic 
Collisions, by James F. Hornig and J. O. 
Hirschfelder, Wisconsin Univ., Naval Res. 
Lub., Rep. ONR-5, Oct. 1953, 31 pp. 

Calculation of the Heats of Formation of 
Gaseous Free Radicals and Ions, by 


When ordinary temperature limits are exceeded... 


CERAMIC COATING 


OF PARTS MADE OF 


STAINLESS INCONEL Oy 
HASTELLOY 


SHOULD BE CONSIDERED. 


Write for complete information and prices. 


LANGDON RD. & PENN. R.R., CINCINNATI 13, 0. Pcev™ 
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BARROWS PORCELAIN ENAMEL CO. 3P: 


iid to work on 


These career positions have been created by Lockheed’s program of 
diversified development — a program that means more scope for your 
ability and more opportunity for promotion because it is diversified. 


In addition to more 
opportunity, you receive 
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DYNAMICISTS and 


AERO- 


nuclear energy 
jet transports 
super-sonic fighters 


continuing development of the Super 
Constellation and other production models 


increased pay rates now in effect peiic® 
generous travel and moving allowances 
the chance for you and your family to 
enjoy life in Southern California 


DYNAMICS ENGINEERS 


Lockheed invites qualified aerodynamicists and aerodynamics engineers 


to apply for these positions. Coupon below is for your convenience. 


LOCKHEE 


Mr. E. W. Des Lauriers, 

Engineering Recruiting, Dept. JP-A-3 
Lockheed Aircraft Corporation 
Burbank, California 


Dear Sir: Please send me an application form and illustrated 
brochure describing life and work at Lockheed in California. 


AIRCRAFT CORPORATION 
BURBANK, CALIFORNIA 


my name 


my street address 


my city and state 
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subminiature 
accelerometer 


iL detect and measure 
static or dynamic 
acceleration in ranges from 
+2 gto +100 g, the Model 
F accelerometer is available 
for applications requiring 
miniature components. 


Weighing less than one 
ounce, this accelerometer 
can be attached to light 
machine parts without 
altering their dynamic 
performance, or it may be 
inserted into small spaces 
previously inaccessible for 
vibration measurements. At 
no sacrifice in performance, 
extremely small size and 
weight are emphasized in 


the Model F accelerometer. 


bine Als. Ge 


3 


Please request 
Bulletin No. 4.0. 


LABORATORIES 
los Angeles 64, Calif. 


J. L. Franklin, J. Chem. Phys., vol. 21, 
Nov. 1953, pp. 2029-2033. 

Prediction of Liquid-Vapor Equilibria 
from an Empirical Equation of State for 
Mixtures, by George Galloway Duplaise 
Town, Wisconsin Univ., Naval Res. Lab., 
Rep. no. ONR-4, Aug. 1953, 14 pp. 

The Intermolecular Potentials for Some 
Simple Non-Polar Molecules, by Ed- 
ward A. Mason and William E. Rice, 
Wisconsin Univ., Navul Res. Lab., Rep. no. 
ONR-6, Nov. 1953, 34 pp. 

Theory of the Formation of Bubbles, by 
Shumsuke Takagi, J. Applied Phys., vol. 
24, Dec. 1953, pp. 1453-1461. 


Instrumentation and 
Experimental Techniques 


Measurement of the Size Distribution 
of Spray Particles, by L. K. Wheeler and 
E. S. Trickett, Electronic Engineering, 
vol. 25, no. 308, Oct. 1953, pp. 402-406. 

An Investigation into the Possibilities of 
Adapting the Hot Wire Anemometer to 
Liquid Phase Measurement, by David T. 
Harrje, Princeton Univ., Dept. Aeron. 
Engng., Rep. no. 237, Oct. 1953,’40 pp. 

Analysis of Transient Aerodynamic 
Heating Problems with GEDA Equip- 
ment, Goodyear Aircraft Corp., Rep. No. 
GER-5434, June 1953, 17 pp. 

The Use of Radar to Measure the True 
Air Speed of an Aircraft, by J. Warner, 
and K. J. Heffernan, J. Roy. Aeron. Soc., 
vol. 57, Oct. 1953, pp. 653-655. 

A Thermistor Calorimeter for Heats of 
Wetting, Entropies from Heats of Wetting 
and Absorption Data, by A. C. Zettle- 
moyer, G. J. Young, J. J Chessick, and 
F. H. Healey, J. Phys. Chem., vol. 57, 
Oct. 1953, pp. 649-652. 

Burning Time and Ignition Temperature 
Apparatus for Metal Powders, by Holger 
C. Andersen and Lawrence H. Belz, Rev. 
Sci. Inst., vol. 24, Oct. 1953, pp. 1004. 

The Polarographic Determination of Tin 
in Propellants, by C. Rubaudo, Picatinny 
Arsenal, Tech. Rep. 1974, Oct. 1953, 6 pp. 

Pressure Response in Supersonic Wind- 
Tunnel Pressure Instrumentation, by 
Arnold L. Ducoffe, J. Applied Phys., vol. 
24, Nov. 1953, pp. 1343-1354. 

Wire Resistance Strain Gage Practice, 
by Chr. Rohrbach, Arch. Tech. Messen. 
no. 213, Oct. 1953, pp. 237-240, 25 refer- 
ences (in German). 

Control System Requirements for Ram- 
jet Missiles, by James C. Wise, Aviation 
Age, vol. 2, Nov. 1953, pp. 64-70. 

A New Space Rate Sensing Instrument, 
by Josseph Lyman, Aero. Engng. Rev., vol. 


| 12, Nov. 1953, pp. 24-30. ‘ 


General Theory and Operational Char- 
acteristics of the Gyratron Angular Rate 
Tachometer, by Roland E. Barnaby, John 
B. Chatterton, and Fred H. Gerring, Aero. 
Engng. Rev., vol. 12, Nov. 1953, pp. 31-36, 
106 


Monitor for Automatic PiJots, by John 
J. Hess, Jr., Electronics, vol. 26, Dec. 1953, 
pp. 174-177. 

Guided-Weapons Telemetering, the 
Two Principal Systems Developed by the 
Ministry of Supply (Gt. Brit.), Flight, 
vol. 64, Oct. 23, 1953, pp. 556-558. 

Radio-Frequency Mass Spectrometer 
for Upper Air Research, by John W. 
Townsend, Jr., Naval Research Lab., Rep. 
No. 3928, Jan. 1952, 8 pp. 

High Frequency Pressure Indicators 
for Aerodynamic Problems, by Y. T. Li, 
NACA TN 3042, Nov. 1953, 52 pp. 

Very Accurate Measurements of Fringe 
Shifts in an Optical Interferometer Study 
of Gas Flow, by Ernst H. Winkler, Rev. 
Sci. Instrum., vol. 24, Nov. 1953, pp. 
1067-1068. 


_ TEMPERATURE 


New, Improved Probe Accurately 
Measures Stagnation Temperature of 
a Rapidly Moving Airstream. 


torn 


MODEL 49126-A ts 


‘ 


Advanced aerodynamic design of the 
Giannini Adiabatic Temperature 
Probe permits rapid, precise temper- 
ature measurements in moving air by 
stagnating samples of Ambient Air 
in a thermo-isolated chamber. Com- 
plete Adiabatic Rise of the stagnated 
air is captured and indicated by the 
Giannini “Ultra-Fast Response” (0.5 
seconds at speeds greater than 100 
mph) Resistance Element, Thermo- 
couple, or Resistance Bulb. Recovery 
factor of better than 0.99 makes the 
probe ideal for True Air Speed in- 
strumentation. A heating jacket to 
provide anti-ice control is optional. 
Write for complete engineering 
information. 


G. M. GIANNINI & Co., INC. 
AIRBORNE INSTRUMENT DIVISION 
PASADENA 1, CALIFORNIA 
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A Goachatnstic Field Angular Orienta- 
tion Indicator for a Rocket, by Donald M. 
Packer, Naval Res. Lab., Rep. 4247 
(Upper Atmosphere Res. Rep. 19) Nov. 
1953, 16 pp. 

Rocket-Borne Servo Tracks the Sun, 
by D.S. Stacey, G. A. Stith, R. A. Nidey, 
and W. B. Pie tenpol, Electronics, vol. 27, 
Jan. 1954, pp. 149 151. 


Terrestrial Flight, Ballis- 
ties, and Vehicle Design 


Taking Man to Higher Altitude, by A. 
M. Mayo, Soc. Auto. Engrs., Prep. no. 
144, Sept. 1953, 8 pp., 11 figs. 

Fabrication of the Orbital Vehicle, by 
K. W. Gatland, A. M. Kunesch, and A. E. 
Dixon, J. Brit. Interplan. Soc., vol. 12, 
Nov. 1953, pp. 274-285. 

A Study of the Applicability of the Un- 
steady One Dimensional Isentropic 
Theory to an Experimental Gun, by A. E. 
Seizel, Naval Ord. Lab., N AVORD Rep. 
2602 (Aeroballistic Res. Rep. 79), July 
1952, 34 pp. 

Design Trends—Guided Missiles, by 
Clark B. Millikan, Aero. Engng. Rev., 
yo!. 12, Dec. 1958, pp. 82-87, 129-130. 

Douglas X-3, ‘‘The Flying Stiletto,’’ 
American Aviation, vol. 17, Dec. 17, 1953, 
p. 17. 

Design of Military Aircraft, by J. H. 
Kindelberger, Aero. Engng. Rev., vol. 12, 
Dec. 1953, pp. 42-49, 125-126. 

A Guide to High-Altitude Rocket De- 
sign, by Richard C. Lea, Aviation Week, 
vol. 59, Dec. 28, 1953, pp. 28-31. 


Space Flight 


Solar Space, Part I. A New Map of 
Solar Space, by Wayne Proell, J. Space 
Flight, vol. 5, Oct. 1953, pp. 1-9. 

Solar Space, Part II. The Probability 
of Planets in Solar Space, by Wayne 
Proell, J. Space Flight, vol. 5, Nov. 1953, 
pp. 1-7. 


(strophysics, Aerophysics, 
and Atomic Physics 


The Propagation of Weak Spherical 
Shocks in Stars, by G. B. Whitham, 
Comm. Pure and Appl. Math., vol. 6, 
Aug. 1953, pp. 397-414. 

Radio Astronomy, by Hans Herbert 
Klinger, J. Franklin Inst., vol. 256, Oct. 
1953, pp. 353-366. 

Air Force Cambridge Research Center. 
Geophysics Research Directorate, Geo- 
physical Research Paper No. 19 (AFCRC 
Technical Report 59-9), Dec. 1952, 530 pp., 
International Symposium on Atmospheric 
Turbulence in the Boundary Layer, 
Massachusetts Institute of Technology, 
June 4-8, 1951. 

Radio Echo Studies of Meteor Ioniza- 
tion, by T. R. Kaiser, Advances in Physics, 
vol. 2, Oct. 1953, pp. 495-544. 

V-2 Rocket in Upper Atmosphere Re- 
search, by Charles F. Green, Aero 
Digest, vol. 67, Nov. 1953, pp. 20-26. 

Measurement of Variations in At- 
mospheric Refractive Index with an 
Airborne Microwave Refractometer, by 
Howard E. Bussey and George Birnbaum, 
4 Res. Nat. Bur. Stand., vol. 51, Oct. 1953, 

171-178. 

The Use a Clouds for Locating the Jet 
Stream, by Vincent J. Schaefer, Aero- 
plane, vol. 85, Oct. 30, 1953, pp. 599-602. 

Astronautics, by Antonio Eula, L’ Aero- 
tecnica, vol. 33, June 1953, pp. 231-244 
in Italian). 
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ENGINEERS 


cae! for Lockheed’s expanding Missile Systems Division 


ae Recently formed from other Lockheed Engineering organizations to 
prepare for the era of automatic flight, the Missile Systems Division 


deals exclusively with missiles and their component systems. 
Its expansion has created “ground-floor” openings for: _ a 


Aerodynamicists 
the 
Thermodynamicists 
Dynamicists 
Design Engineers “A” and 
Electronic Engineers 


(micro-wave techniques, electronic 
components, circuits or design flight 
instrumentation) 


In addition to outstanding career opportunities, the Lockheed 
Missile Systems Division offers you high salaries commensurate with —_ 
your experience, generous travel and moving allowances, and a 
better life for you and your family in Southern California. 


Address inquiries to L. R. Osgood, Dept. JP-M-3 Lockheed Missile 
Systems Division, 7701 Woodley Avenue, Van Nuys, California. 


Electro-Mechanical Engineers 
(circuits or servomechanisms) 


Applied Mathematicians 
and Physicists 


Mechanical Engineers 


(precision mechanisms or packaging 
of electronic equipment) 


LOCKHEED 


VAN NUYS, CALIFORNIA | DIVISION | 
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World's Largest Stretch-Wrap Forming Machine 


The Model 60 is Hufford’s largest stretch-wrap 
forming machine to date, exceeding the capacity of 
all existing stretch-wrap forming machines in the 
world today. 


With its 703,720 pounds of “pull” it is now form- 
ing stabilizer fairings, panels and fuselage skin parts 
for the new FJ-2 FURY. New dies will utilize its 
huge capacity on sheets measuring 6 feet wide by 


ullord_ 


MODEL 60 
CAPACITY: 350 TONS 


to work for North American Aviation Inc., Columbus, Ohio 


31% feet long. Its tremendous strength enables form- 
ation of sheet thicknesses as great as % inch! 


The Model 60 opens another chapter in stretch- 
wrap forming history, almost doubling the capacity 
of the largest Hufford press heretofore built. It is 
Hufford’s answer to unusual requirements of the 
present in preparation for the commonplace of the 
future! 
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You are looking at 1/1000 of a second in the history of 
aviation. It occurred at a fraction past 4:31 p.m. on 
January 20, 1949. 


This was the Zero moment which marked the official 
launching of the first successful pilotless bomber to be 
approved by the U.S. Air Force—the Martin B-61 
Matador. 


The picture is historic for a very significant reason: it re- 
cords the tradition-shattering payoff of an entirely new 
development in the aircraft industry, known as Martin 
ystems Engineering. This is a science and a method of 


history | 


developing spaceborne systems as total solutions of Oper- 
ations problems. 


The Martin Matador is far more than the thing you 
glimpse here. Behind it is an integrated network of fa- 
cilities designed to give this important new weapon sim- 
plicity of operation and extreme mobility. These 
components add up to the total solution of one of the 
most formidable security problems of our time. 


They also add up to one of today’s most important devel- 
opments: the full story of Martin Systems Engineering. 


You will hear more about Martin! 


THE GLENN L. MARTIN COMPANY 


BALTIMORE MARYLAND 


Mort. 
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NITROGEN 
TETROXIDE 


offers outstanding advantages 


to designers of f 


- wi SPECIFIC IMPULSE: Nitrogen Tetroxide exceeds many other well- 
known oxidizers in pounds of thrust developed per pound 
of fuel consumed per second. 

EASY T0 HANDLE: Nitrogen Tetroxide may be shipped, piped and 
stored in ordinary carbon steel equipment. It possesses 
high chemical stability, high density, low freezing point, 
and a reasonably low vapor pressure. 


Nitrogen Tetroxide is available at low cost in 125-pound I.C.C. approved 
steel cylinders and 1-ton containers. 


ones _ Address your inquiry to the Product Development Department 


abi tits 


40 RECTOR STREET, NEW YORK 6, N. Y. 


Technical service and development on Nitrogen Tetroxide—formerly handled by the Product Development Department, 
Solvay Process Division—are now handled by Nitrogen Division, Allied Chemical & Dye Corporation. 
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The 66W Series of STRAYOPOWER Hydraulic 
Pumps embrace a range of models from 2 
te 10 gpm at 1500 rpm with continuous pres- 


continuous speed of 3750 rpy 

speeds to 4500 rpm. 
4) 


wres to 3000 psi. Designed for maximum a 


STRATOPOWER 


HYDRAULIC PUMPS 


To the layman, the pinpointing 
of a destination for a guided mis- 
sile is still pure magic. But to the 
engineer, the control and guid- 
ance of rockets and guided mis- 
siles in flight simply means 
another application for 
STRATOPOWER quality and 
advanced design. 


The unerring performance of 
STRATOPOWER Pumps has 
been demonstrated times with- 
out number down through the 
years. At sea level and at 
heights still to be achieved these 
perfect examples of precision 
engineering provide the fluid 


Get the full story on STRATOPOWER 
constant and variable delivery Pumps 
for your hydraulic circuits. 


ta 
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ent, 
on. STARBUCK AVENUE 


Marcu-Aprib 1954 


N 


WATERTOWN: N.Y. 


power that positively answers 
the question of weight vs. horse- 
power as well as the equally 
important requirement of long- 
lived dependability. 


In the designing of any high 
pressure hydraulic circuit there 
are definite advantages in 
STRATOPOWER Pumps. 
Whether yours is a problem for 
constant or variable delivery, 
high or low temperature opera- 
tion, capacities from .25 to 30 
gpm, STRATOPOWER will pro- 
vide the pump to 3000 psi that 
will resolve that problem NOW! 


BOUND Controlled All the Way by= 


WATERTOWN DIVISION 
THE NEW YORK AIR BRAKE COMPANY 


730 Starbuck Avenue « Watertown, N. Y. 


Kindly send me information on STRATOPOWER 


Hydraulic Pumps 


(C0 Constant delivery [] Variable delivery 


Name 


Company. 


Address 
City. 


Zone___ State. 
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RANOF 


Consolidated Engineering 


CORPORATION 
300 North Sierra Madre Villa, Pasadena 15, California 


FOR 


Sales and Service through €E€ Instruments, Inc., a subsidiary with offices in: Pasadena, Chicago, Dallas, Detroit, New York, Philadelphia, Washington, D. a 
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ANALYTICAL INSTRUMENTS 


in many research and development projects 
today, full scale structures or scale models are 
tested and re-tested hundreds of times with only 
slight changes in the test conditions. Where this 
occurs, as it does in aerodynamic wind tunnels 
or rocket and reaction motor static-test stands, a 
staggering amount of test data may be accumu- 
lated in analog form. Before this data can be 
evaluated by the test engineer, it must be con- 
verted into digital—i.e., numerical—form. CEC’s 
high speed, high accuracy SADIC digital data- 
processing systems are un- 
equalled for this purpose. 
Assembled from “building 
bleck” components to meet 
the precise requirements of 
the test project, they auto- 
matically convert measure- 
ments from thermocouples 
and resistance pickups into 
numerical form on paper 
tape, punch cards or tabula- 
tion machines immediately 
available on completion of the test. Where large 
quantities of such test data have to be processed, 
automatic digital data conversion with a Con- 
solidated SADIC system pays for itself many 
times over by eliminating the costly, laborious, 
time-wasting necessity for extracting thousands 
of discrete points by hand from an enormous 
accumulation of analog plots. 


Rik W3H 


SCIENCE AND INDUSTRY 


Jet PRopuLs!0N 


Nort 

Batt 
Cali, 
Nort! 
West 


Pac ITO 
Geor 


React 
Lond 
REyNc 
ING ( 
RYAN 
Batte: 
Calif. 
Srarua 
W este 
STURTE 


Warver’ 
H mb 


— 


Marcy 


AV 
a 
y 
Bo 
re 
Co 
e 2 Do 
H J 
es 
K 
; Exe 
Far 
G 
Gi 
For 
G. 
GEN 
De 
Gta) 
Gyr 
Co 
Hur 
Ch, 
Kau 
Gec 
Kou! 
Sch 
Lavo 
The 
Lock 
M1 
Hal 
Hei 
% 
4 
- 
SUMME 
ELECC 
9 tre 
TENNE 
th 0.8. 
Trre: 
EFLI 
= 


+t+tt +t 


Index to Advertisers 


AgroJET-GENERAL Corp. 


outside back cover 


AIRnORNE ACCESSORIES CoRP......... 132 
Gray and Rogers, Philadelphia, Pa. 

Auuicp CuemicaL & Dye Corp. 
Alhert Sidney Noble "Advertising, New York, 


Doyle, Kitchen & McCormick, New York, N. 


Buller & Combes, Providence, R. I. 
Barrows PorcELAIN ENAMEL Co...... 133 
Pe ry Brown, Inc., Cincinnati, Ohio : 
Benorx AviATION CorP. 
ScintiLLA MAGNETO Div..........- 116 
The Shaw Co., Los Angeles, Calif. a : om 
M. Carty Co., Los Angeles, Calif. 
CoNSOLIDATED ENGINEERING CorP.... 140 
Hirson & Jorgensen, Inc., Los Angeles, Calif. 
Dov Las AIRCRAFT Co, 7 
J. Walter Thompson Co., Los Angeles, Calif. ase 
ELE: 112 an 
Kcy-Christopher, Hollywood, Calif. 
ExceLco DEVELOPMENTS, INC....... 142 
ENGINE AND AIRPLANE 
ORP. 
Gaynor & Co., Inc., New ‘ork, N.Y. ++ 
Forp INSTRUMENT 70 +4 it 
G. M. Basford Co. New ‘York, 
Deutsch & Shea, New York, N.Y. H i = ange 
GianninI, G. M., & Co. INc........... 134 
Western Advert. ‘Agency, Inc., Los Angeles, Calif, ttt 
GyROMECHANISMS, Inc... . .inside back cover ss: 
Corydon M. Johnson Co. Inc., Bethpage, N. Y. ee | Complete—all calibrating circuits built-in. 
Hurrorp MAcHINE Works, INc.. . 136 N 
ude D. Graham, Lee Aupiles, Calif.” | eeds only standard DC oscilloscope. 
Kavpp, C. B., & Sons. Developed by Fairchild in its program for transistorizing Missile 
George Homer Martin Associates. Guid 
vidance Systems and other intricate electronic circuits. 
Schaefer and Favre, New York, N.Y. 
144 Rapidly Plots Static and Dynamic Characteristics of ALL Transistors 
The Picard Advertising Co., New York, N. Y. ... point contact and junction. Designed on basic principles, to meet 
LocKHEED AIRCRAFT CorRP............ 133 
MIssILe Systems Div.............. 135 uture transistor needs. 
Hal Stebbins Inc., Los Angeles, Calif. 
citi die Ce.......... 127 Complete Families of Curves obtainable in 10 incremental steps for 
Heintz & Co.,Inc., Los Angeles, Calif. 
rInc., , each of 5 ranges. Sweeping technique shows up anomalies. 
Martin, THE GLENN L., Co.......... 137 9 ping P 
Vansant, Dugdale & Co., Baltimore, Md. 
= MINIATURE PreEcISION BEARINGS INC.. 118 
Packard & Kraft, Worcester, Mass. 
N.Y. Arr Brake Co. 139 Presents on the Scope: 
Humbert & Jones, New York, N.Y. 
NortH AMERICAN AVIATION.......... 76 Alpha vs. Emitter Current 
Batten, Barton, Durstin & Osborn, Los Angeles, x 
Calif. Collector, Emitter and 


NortHrop Aircrart, Inc. .......... 
West-Marquis, Inc., Los Angeles, Calif. 
PHorocon RESEARCH PRopvucTs...... 117 

George Burtt, Hollywood, Calif. 
Reaction Morors, Inc.. .inside front cover 
London Advert. Agency, Newark, N. J. 


Transfer Characteristics 


Collector Characteristics 
in Grounded Emitter Connection 


{EYNOLDS ELectTrRIcAL & ENGINEER- 
Ryan AERONAUTICAL CO........... 73 ¢ 
Batten, Barton, Durstin & Osborn, Los Angeles, Al Sc } LD 
Calif. 
STATHAM LABORATORIES............. 134 7 
Western Advert. Agency, Los Angeles, Calif. Cuided, (ttiles Dvition 
STURTEVANT, P.A.,Co............... 132 
Ross Llewellyn, Inc., Chicago, Ill. Wyandanch, L. 1., N.Y 
Summers Gyroscope Co............. Other Divisions: Ai ivisi m, Md. i ivisi i 
TELECOMPUTING CoRP............... 120 |° ° 
. Hal Stebbins, Inc., Los Angeles, Calif. x FAIRCHILD ENGINE & AIRPLANE CORPORATION : 
ENNEY ENGINEERING, INC........... 72 e 
TirerLex,INc.............. .. 143 |* Please send me your Detailed Transistor Analyzer Technical Bulletin. ° 
Join Falkner Arndt, Philadelphia, ‘Pa. 
| 
Wavertown Div., N. Y. AIR Baars Co. 139 |; Name and Company. 
Humbert & Jones, New York, N. 
Address 
City Zone State 
N | Marcu-Aprin 1954 141 


4, 


EXCELCO DEVELOPMENTS 


: INCORPORATED 
MILL STREET - BOX 230 - 
Me 


For Skill And Precision 
In The Development 
And Manufacture Of... 


ROCKET MOTORS 


and 
GUIDED MISSILE COMPONENTS 


Fabricators of 


COMPLETE ROCKET MOTORS 
SPHERES FOR PRESSURE TANKS 
NOZZLES OF ALL TYPES 

INNER & OUTER THROAT SECTIONS 
AIRCRAFT SEATS & BULKHEADS 
ELECTRONIC CHASSIS & DETAIL ASSY 
LOX BOILERS 
TORUS TANKS 
SPECIAL MACHINING 


vot Aer 


TESTING AFTER FABRICATION 


DEVELOPMENTS 


j 
ix 
| 
— 
: HELL-ARC WELDING FORMING OF STOCK TO PRECISION DIMENSIONS 


teflex Connector 


Problems 
New Modifiable 


tr 


We have the experience to solve most complex connector 
problems involving extreme altitudes, temperatures and 
pressures with space and weight limitations. And within standard 


INSTRUMENTATION Connector —07. For moisture 
and corrosion resistance, temperature ranges of 
—65°F. to +400°F. Made of Teflon, plug and 
receptacle mated weigh only 74 oz. Length 2". 
Insulation properties will permit 3500 volts at sea 
level, 1200 volts at 50,000 feet altitude. Can be 
made with 2 or 3 pins, current 7 amperes. 


CUSTOM WIRING SYSTEMS—For accessory, 

strumentation and radio shielded applications. Can 
be furnished with Titeflex or Standard AN Con- 
nectors. Can be sheathed with one or more layers 
of various metal braids, fiber glass or nylon, and 
jacketed with silicone or various other compounds. 
Titeflex will be glad to design, develop and pro- 
duce complete wiring systems to your specifications. 


= 


design requirements, we can develop special-duty connectors 


as a part of complete wiring systems. 


MOISTURE-PROOF and resistant to synthetic lubri- 
cants. For extreme temperature changes in ranges 
of —65°F. to +400°F., high altitudes up to 65,000 
feet. Resists salt spray, corrosion, vibration. This 
Titeflex Connector is radio shielded, has positive 
retention of pins and sockets. 5" in length. Mates 
with connectors that conform to MIL-C-5015. 


SPECIAL —O7 CONNECTOR. Designed to solve 
your connector problems in instrumentation with a 
real saving in space and weight since this connector 
has no protuberance beyond the flange. Can be 
designed as an integral part of your wiring or 
instrument components. Available in 1, 2 or 3 pin 
arrangements—current 7 amperes, size | "in length. 
Receptacle and plug weigh only 11 grams. 


All TITEFLEX Connectors can be furnished with thermocouple pins and sockets. 


TITEFLEX, INC. 
578 Frelinghuysen Avenue, Newark 5, N. J. 


WRITE TODAY for specific information—or send 
us your specifications. Whatever your requirements, | 
we can usually provide the right answer. Our Engi- 
neering Staff will be glad to | 
discuss your problem without | Please send me your catalog on the Titeflex Connector. 
obligation. | 


Have your representative call 0. 
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FAMILY OF PRODUCTS a 
HELP YOURS 


FREQUENCY METERS 

Three frequency meters accurate to 0.001% cover 
ranges from 10 to 2000 MC. Model LA-5 covers the 
10 to 100 MC range, LA-6, 100 to 500 MC and LA-61 
500 to 2000 MC. 


PULSE GENERATOR—MODEL LA-592D 


A double-pulse generator with wide 
range control, excellent pulse shape. 
Eliminates necessity for many instru- 
ments usually required ir in an electron- 


ics laboratory. 


OSCILLOSCOPE—MODEL LA-239C 

The new and improved Lavoie oscilloscope offers wider 
frequency range, greater sensitivity and faster rise time. 
The Lavoie Camera Adapter may be added quickly and 
without modification. 


PRECISION CRYSTAL OVEN 75° or 85° 
Maintains crystal temperature to 
within 0.025° C at normal room tem- 
perature—to within .15° C over outside 
range from -40° F to 150° F. Available 
for HC 6 or HC 183 crystal units. 
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MORGANVILLE, NEW JERSEY 
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